Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,478,245 C→A 8.8% L310M (CTG→ATG)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,478,2450CA 8.8% 184.4 / 3.0 67L310M (CTG→ATG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base A (4/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGACCGTACCGGG  >  CP009273/1478102‑1478393
                                                                                                                                               |                                                                                                                                                     
cGATTGTTGAAGAAGCCGATGACACTGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAAc                                                                                                                                                  >  1:639717/1‑149 (MQ=255)
 gATTGTTGAAGAAGCCGATGACACCTAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACa                                                                                                                                                 >  1:621825/1‑149 (MQ=255)
        gaagaaGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGGGCGGGAAATCCGCGATACCCCCGATGCGATGAACAAGCTGaa                                                                                                                                          >  2:665059/1‑149 (MQ=255)
        gaagaaGCCGATGACACCGAGAGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGaa                                                                                                                                          >  1:878100/1‑149 (MQ=255)
         aagaagCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAAc                                                                                                                                         <  1:805194/149‑1 (MQ=255)
             agCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTAc                                                                                                                                     <  2:1157129/149‑1 (MQ=255)
              gCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATGACGCCGATGCGATCAACAAGCTGAACTTAca                                                                                                                                    >  1:565490/1‑148 (MQ=255)
               ccGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGc                                                                                                                                   <  1:569342/149‑1 (MQ=255)
                 gATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGCGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCa                                                                                                                                 >  1:337986/1‑149 (MQ=255)
                         cgtgcgcGATCAGTGGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATTAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGaaa                                                                                                                         <  1:750115/146‑1 (MQ=255)
                         cGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGGACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGaaa                                                                                                                         >  1:1036749/1‑149 (MQ=255)
                          gAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGACCAAGCTGAACTTACGCCATACCGAAAt                                                                                                                        >  2:498393/1‑149 (MQ=255)
                                    cgTTGCAGGCGA‑TTTTGACGCCTTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCttt                                                                                                             <  2:164406/148‑1 (MQ=255)
                                       tGCAGGCCATTTTTGACGCCGTAGACGAACCGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTAt                                                                                                           <  1:15362/149‑1 (MQ=255)
                                            gCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGcgcg                                                                                                      <  2:1108564/149‑1 (MQ=255)
                                               aTTTTGTACGCCGTAGACGACCTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCt                                                                                                   <  1:513102/149‑1 (MQ=255)
                                                tttttGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCAATACAGCAATCTTGCCGCTTTATGCGCGGAtt                                                                                                  >  1:168904/1‑149 (MQ=255)
                                                 ttttGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                                                                                                 <  2:399833/149‑1 (MQ=255)
                                                 ttttGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGATGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                                                                                                 <  1:1228464/149‑1 (MQ=255)
                                                 ttttGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGATGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                                                                                                 <  1:318184/149‑1 (MQ=255)
                                                 ttttGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGATGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                                                                                                 <  1:318708/149‑1 (MQ=255)
                                                 ttttGACGCCGGAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTTAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCttt                                                                                                 <  2:240742/149‑1 (MQ=255)
                                                 gttgtactcgggaGATGAACTGAGTCAAGGGAGGCGTTGCGACATTCGGATCTTTAGGACCGGCGAGAGGGGAATCCGCGATACCGCCGATGCGCTGCACAAGCGGAACGATCACACTACCGAAATCTTCCCGCTTTATGCGAGGCtct                                                                                                 >  2:423766/13‑147 (MQ=255)
                                                    tgccggcGTAGACGAATTGAGTCAGGAAAGCCATGGCGACATTCTGACCTTTATGAGCTGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGACCAAGCTGAACTCAATCTTTACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGa                                                                                              <  2:519448/143‑1 (MQ=255)
                                                      aCGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAAc                                                                                            >  2:435553/1‑149 (MQ=255)
                                                         ccGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGc                                                                                         <  2:472375/149‑1 (MQ=255)
                                                         ccGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGc                                                                                         <  1:776551/149‑1 (MQ=255)
                                                            tcgACGCACTGAGTCAGGAACGCCATGGCGAGATTCCGATCCTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGATCAAGCTGAACTTACGCCATCCCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGaa                                                                                      <  1:418491/147‑1 (MQ=255)
                                                              gACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACa                                                                                    >  1:1118682/1‑149 (MQ=255)
                                                               aCGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGGACaa                                                                                   >  1:677592/1‑149 (MQ=255)
                                                               aCGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACaa                                                                                   >  1:812328/1‑149 (MQ=255)
                                                                     tGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAACCTACTATTTCGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAgg                                                                             >  1:1214647/1‑149 (MQ=255)
                                                                     ggggTCAGGAAAGCCATGGCGGCGTTCGGATCTTTATGGTCGGCGAGCGGGGAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCGTACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACGGCGAACAAAATAgg                                                                             <  2:403713/146‑1 (MQ=255)
                                                                     cgAGTCAGGAACGCCATGGCGACATTCTCATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAgg                                                                             <  2:403427/148‑1 (MQ=255)
                                                                      gAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTGTACGCGGCTTTCGAACAGCGAACAAAATAggg                                                                            <  1:607929/149‑1 (MQ=255)
                                                                        gTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTa                                                                          <  2:1218912/149‑1 (MQ=255)
                                                                         tCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGGCGATGCGCTGAACAAGCTGAACTTAGGATTTATTTTAAGCTTTCCGCTTTATGTGCCGTTTTTGAACAGCGAAATAAATAGGGTAt                                                                         >  1:35161/1‑149 (MQ=255)
                                                                         tCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGTTAt                                                                         >  1:876488/1‑149 (MQ=255)
                                                                         tCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTAt                                                                         >  2:165240/1‑149 (MQ=255)
                                                                           aGGAAAGCCATGGAGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCGGAACTTACGCCATACCGAACTCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTc                                                                       <  1:519557/149‑1 (MQ=255)
                                                                               aggccagtgCGACATTCTCATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGAGGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGt                                                                   <  1:365735/141‑1 (MQ=255)
                                                                                aggcATGGCGACATTCTGATCTTTATGAGCGGCGCGCGGGAAATCCGCGATACCGCCGTTGCGCTGAACCAGCTGACCCGACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGATATTCCAGTc                                                                  <  2:493326/146‑1 (MQ=255)
                                                                                 gCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCCGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACGTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCg                                                                 <  1:245732/149‑1 (MQ=255)
                                                                                  ccATGGCGACATTCTGATCTTTATGAGCGGAGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGc                                                                >  1:642197/1‑149 (MQ=255)
                                                                                    aTGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGcac                                                              >  1:583697/1‑149 (MQ=255)
                                                                                        cGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGTTTTCGAACAGCGAACAAAATCGGGTATTCCAGTCGCACAGCg                                                          >  2:493388/1‑149 (MQ=255)
                                                                                         gACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCgg                                                         >  1:634566/1‑149 (MQ=255)
                                                                                                tGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGg‑cgc                                                 >  2:627763/1‑149 (MQ=255)
                                                                                                  aTCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACACAATAGGGTGATCCCGTCGCACAACGGACGG‑CGCAt                                               >  1:1028805/1‑149 (MQ=255)
                                                                                                      ttATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATtgtg                                           <  2:440209/149‑1 (MQ=255)
                                                                                                               gcgAGCGGGAAATCCGCGATACCGCCGCTGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGAcc                                  <  1:67322/149‑1 (MQ=255)
                                                                                                                    cGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGt                             <  1:498393/149‑1 (MQ=255)
                                                                                                                    cGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCACTGTGCTGGCGACCAACgg                             >  2:86836/1‑148 (MQ=255)
                                                                                                                    cGGGAAATCCGCGATACCGCCGATGCGATGGAGGATCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGt                             >  2:881857/1‑149 (MQ=255)
                                                                                                                      ggAAATCCGCGATCCCGCCGATGCGCTGAACAAGCTGAACTCACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTcg                           <  2:172673/149‑1 (MQ=255)
                                                                                                                      ggAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAACCTCGCCGCTTTATGCGCGGCTTTCGAAAAGCGAACACAATATGGCACTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTcg                           <  1:1166939/149‑1 (MQ=255)
                                                                                                                      ggAAATCCGCGATACCGCCGATGCGCTGAACAAGACGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTcg                           >  2:267004/1‑149 (MQ=255)
                                                                                                                      ggAAATCCGCGATACCGACGATGCGATGAAACAGCTGAACTTACGCCAAAGTGAAATCTTGCCGCTTTAGGCGCGGCTTGCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTcg                           >  2:247435/1‑149 (MQ=255)
                                                                                                                      cgAACTCCGCGATCCCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTcg                           <  2:408027/148‑1 (MQ=255)
                                                                                                                         aaTCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCCATGCGCTGACGACCAACGtgacgg                        >  1:334324/1‑144 (MQ=255)
                                                                                                                           tCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAACTAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATGGTGCTGGCGACCAACGTCGCGGaa                      >  2:776271/1‑149 (MQ=255)
                                                                                                                            ccGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCCGaaa                     >  1:667992/1‑149 (MQ=255)
                                                                                                                              gcgATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACg                   >  2:134790/1‑149 (MQ=255)
                                                                                                                                     cgccgATGCGCTGAACGAGCTGAACTTACGCCATACCGGAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAACATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGa            <  2:150649/149‑1 (MQ=255)
                                                                                                                                      gccgATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGAc           <  1:189282/149‑1 (MQ=255)
                                                                                                                                      gccgATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGAc           >  2:764124/1‑149 (MQ=255)
                                                                                                                                          aTGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGCGCCGCCGCCCAACGCCGCGGAAACGCCGCtccccgca       >  2:600520/1‑142 (MQ=255)
                                                                                                                                             cgcTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGACCGTACCg    <  1:528444/149‑1 (MQ=255)
                                                                                                                                             cgcTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGACCGTACCg    <  2:230067/149‑1 (MQ=255)
                                                                                                                                               cTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGAAGGGTGGAATTGT‑CTGTCAAAAAAAGTCTCTTATACGTCGCTGACCGTACCggg  >  1:1207954/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGATTGTTGAAGAAGCCGATGACACCGAGCGCGATCAGTTGCAGGCGATTTTTGACGCCGTAGACGAACTGAGTCAGGAAAGCCATGGCGACATTCTGATCTTTATGAGCGGCGAGCGGGAAATCCGCGATACCGCCGATGCGCTGAACAAGCTGAACTTACGCCATACCGAAATCTTGCCGCTTTATGCGCGGCTTTCGAACAGCGAACAAAATAGGGTATTCCAGTCGCACAGCGGACGG‑CGCATTGTGCTGGCGACCAACGTCGCGGAAACGTCGCTGACCGTACCGGG  >  CP009273/1478102‑1478393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑