Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,037,805 A→G 7.7% intergenic (+64/+72) dppA‑IV → / ← PP_5284 dipeptide transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,037,8050AG 7.7% 372.5 / 14.3 130intergenic (+64/+72)dppA‑IV/PP_5284dipeptide transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (52/64);  new base G (6/4);  total (58/72)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 131 total aligned reads displayed.

CGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAGCAC  >  NC_002947/6037656‑6037948
                                                                                                                                                      |                                                                                                                                                
cgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGa                                                                                                                                                  <  1:2554756/150‑1 (MQ=255)
  aCCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACg                                                                                                                                                <  1:979349/150‑1 (MQ=255)
     ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAGGGGTATCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                             >  2:631387/1‑150 (MQ=255)
     ttGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                             >  1:871150/1‑150 (MQ=255)
     ttGCTGCGCCAAGGCGTCGAGGGGTACCAGATGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                             <  2:2368943/150‑1 (MQ=255)
     ttGCTGCGCCAAGGCGTCGAGGGGTACCAGATGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                             <  2:1083308/150‑1 (MQ=255)
      tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  1:5033832/150‑1 (MQ=255)
      tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  1:5033994/150‑1 (MQ=255)
      tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  1:1400079/150‑1 (MQ=255)
      tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  2:195430/150‑1 (MQ=255)
      tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  1:847076/150‑1 (MQ=255)
        ctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAccc                                                                                                                                         <  1:1900840/150‑1 (MQ=255)
            gcCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                                     >  1:1817923/1‑150 (MQ=255)
            gcCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                                     >  1:1817938/1‑150 (MQ=255)
             ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTa                                                                                                                                    <  2:2915134/150‑1 (MQ=25)
             ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAAACCGCGCCTACAGGGGTAGCGCCAGCCGGGCCGGTGGCGCGTccccgggg                                                                                                                                    >  1:1002605/1‑146 (MQ=255)
              cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAg                                                                                                                                   <  2:2108178/150‑1 (MQ=255)
                aGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGa                                                                                                                                 <  1:1358241/150‑1 (MQ=255)
                    gTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCgg                                                                                                                             <  2:1692603/150‑1 (MQ=255)
                     tCGTGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCggg                                                                                                                            <  1:3486690/150‑1 (MQ=21)
                      cGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGt                                                                                                                           <  1:5496579/150‑1 (MQ=255)
                      cGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGt                                                                                                                           <  2:5325511/150‑1 (MQ=255)
                       gAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGtt                                                                                                                          <  1:3081668/150‑1 (MQ=255)
                          gggTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTCCCGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTAc                                                                                                                       <  1:2370596/150‑1 (MQ=21)
                              aCCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCcgc                                                                                                                   >  1:5149721/1‑150 (MQ=255)
                                cAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTCGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                 >  2:1209569/1‑150 (MQ=255)
                                      aGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCCTGTAGGAGCGGTTTTACCCGCGAAAggg                                                                                                            <  2:1314308/150‑1 (MQ=255)
                                        cccGTTGGGGGGGCTGGATTGCAGCACGGTGAGGGTGACCACGTACGCCTGTACCGGCCTCTGCGGGGGTAAACCCGCTCCTACCGGGGTAGCGCCAGCCGGGATGTGGACGCTACCCCTGTCGGCGCGGTTTTACCCGCGAAAGGGccc                                                                                                         <  2:142784/150‑1 (MQ=255)
                                        cccGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           <  1:4432655/116‑1 (MQ=255)
                                        cccGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTAc                                                                                                                                           >  2:4432655/1‑116 (MQ=255)
                                        cccGTGCGGGCGGCTGGATTACAGCACGGTGACGGTGACCACGTACGCCTTTACCGGCCTCTTCGCGGGTAACCCCGCCCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGccc                                                                                                         <  2:142814/150‑1 (MQ=37)
                                         ccGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAATCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCg                                                                                                        >  2:4497423/1‑150 (MQ=255)
                                                gggggcTGGATTTCAGCATGGTGTCGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGGGGTTAAACCCGCTCCTACAGGGTTAGCGCCAGCCTGGCAGTTGACGCAACCCCTGTAGGAGGGGTTTTACCCGCGAAAGGGCCCGTACTGAc                                                                                                 <  1:428696/147‑1 (MQ=255)
                                                      tGGATTTCAGCACGGGGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTAGCGCCCGCCTGGCAGTTGACGCTACCCCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCCGTACTGACCCtata                                                                                           <  2:4283046/150‑1 (MQ=38)
                                                      tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCGCCAACAGGGGTAGCGCCCGCCGGGCAGTTGAAGCGACCCCGGTAGGAGCCGGGTTAACCGCGCAAGGGCCCGTACTGACCCaata                                                                                           >  1:860108/1‑150 (MQ=37)
                                                        gATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAAACCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAACGGCCCGTACAGACCCTATAAg                                                                                         >  1:3635690/1‑150 (MQ=255)
                                                               gCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCCCTCCGACAGGGGTAGCGCCAGCCTGGCCGTTGACCCTACCCCTGTAGGAGCGGGTTTAACCGCGAAAGGGCCCGGACTGACCCGATAAGGCTCTCa                                                                                  >  2:2341669/1‑150 (MQ=38)
                                                                aaaGGTCACGCTCACAAACAAAGCCTCAACCCCCCTCTTCGCGGTTAAACCCGCACAAACAGGGGTAGCGCAACCCTGCCATTTCACCCAACCCCTTTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAg                                                                                 <  2:509213/149‑1 (MQ=255)
                                                                      gacggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCATTTGACGCTACCCCTGTTGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTcc                                                                           <  1:5529806/150‑1 (MQ=255)
                                                                        cggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCt                                                                         <  2:4525282/150‑1 (MQ=255)
                                                                          gtgacCAAGTAAGCCTGTACGGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCa                                                                       <  2:2444027/150‑1 (MQ=255)
                                                                            gacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGa                                                                     >  2:5496056/1‑150 (MQ=255)
                                                                            gacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCAAGAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGGAGGAGCGGGTTTACCCGCGCAAGGGCCCGTACTGACCCTATAAGGCTCGCCGTACTCCCTCAGa                                                                     >  2:3064378/1‑150 (MQ=39)
                                                                                   tAAGCCTGTACCGGCCTCTTCGCGGTTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCATTTGACGCTACCCCTTTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAAt                                                              <  1:451553/150‑1 (MQ=255)
                                                                                   tAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACTGGGTTAGCGCCAGCCTGGCTTTTGACGCTACCCGTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAAt                                                              <  1:451548/150‑1 (MQ=38)
                                                                                    aaGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCAACAGGGGTAGCGCCAGCCTGGCCGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTCCTGACCCTATAAGGCTCTCAGTACGCCCTCAGACCCCACTc                                                             >  2:467029/1‑150 (MQ=38)
                                                                                      gCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCa                                                           >  2:2760648/1‑150 (MQ=255)
                                                                                       cgtGTACCGGCCTCTTCGCGGGTAAAGCCGCTCCTCCCGGGGTAGCGCCAGCCTGGCTGTGGACGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAc                                                          <  2:5204605/148‑1 (MQ=38)
                                                                                       ccTGTACCGGCCTCTTCGCGGTTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAc                                                          <  2:1498729/150‑1 (MQ=255)
                                                                                       ccTGTACCGGCCTCGTCGGGGTTAAACCCGCTCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAc                                                          <  1:3936368/150‑1 (MQ=38)
                                                                                        cTGTACGGTCCTCTTGGCGGTTAAACCCGTTCCTACTGGGGTAGCGGCAGCCTGGAATTGGACGCTACCCCTGTTGGAGCGGTTTTACCCGCTAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAcc                                                         <  2:1434805/150‑1 (MQ=255)
                                                                                                   actgcgcgGGAAAACCCGCTCCTCACGGGTTAGCGCCAGCGTGCCAGTGGACGCCCCCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGGCCACAATCCACCCATGCTCGAcc                                              <  1:1135519/146‑1 (MQ=255)
                                                                                                    tttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCa                                             >  2:183109/2‑150 (MQ=18)
                                                                                                      tCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATc                                           <  2:1817923/150‑1 (MQ=255)
                                                                                                      tCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATc                                           <  2:1817938/150‑1 (MQ=255)
                                                                                                           ggTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCCCCCATGCTCGACCATCGACaa                                      >  2:4047259/1‑150 (MQ=255)
                                                                                                            gtcaaCCCGCTCCTACCGGGGTAGCGCCAGCCTGGCTGTTGCCGCGACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCTTCGACAAc                                     <  2:3883886/147‑1 (MQ=255)
                                                                                                            gtcaaCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAAc                                     <  2:4088620/147‑1 (MQ=255)
                                                                                                            gTAAACCCGCTCCTACTGGGGTAGCGCCAGCCTGGCTGTTGACGCTACCCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAAc                                     <  2:894207/150‑1 (MQ=255)
                                                                                                             tAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACg                                    >  2:1667863/1‑150 (MQ=255)
                                                                                                               aaCCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGc                                  >  1:2210904/1‑150 (MQ=255)
                                                                                                                  ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACACTCCACCCATGCTCGACCATCGACAACGGCTcc                               >  2:2587929/1‑150 (MQ=255)
                                                                                                                  ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                               >  1:2678385/1‑150 (MQ=255)
                                                                                                                  ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                               >  1:3668861/1‑150 (MQ=255)
                                                                                                                    gCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcccc                             >  2:4828787/1‑150 (MQ=255)
                                                                                                                        cTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCa                         >  1:351688/1‑150 (MQ=21)
                                                                                                                         tACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAcc                       >  2:1989520/1‑150 (MQ=18)
                                                                                                                         tACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCAGGCTCGACCATCGACAACGGCTCCCCaccacc                       >  2:4376520/1‑150 (MQ=18)
                                                                                                                         tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                        >  1:1216920/1‑150 (MQ=255)
                                                                                                                         tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                        >  1:4539212/1‑150 (MQ=255)
                                                                                                                         tACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                        >  1:3816516/1‑150 (MQ=255)
                                                                                                                         tACAGGGGCAGCCCCAGCCTGGCACTACACCCAACCCCTGAAGGACCGGGAATACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCACAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                        <  2:2841150/150‑1 (MQ=255)
                                                                                                                           cAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                               <  1:2715494/141‑1 (MQ=255)
                                                                                                                           cAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCa                      >  2:548751/1‑150 (MQ=255)
                                                                                                                            aGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCaa                     >  2:12807/1‑150 (MQ=255)
                                                                                                                             ggggTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAat                    <  2:3281416/150‑1 (MQ=255)
                                                                                                                              gggTAGCGCCAGCCTGGCTGTTGACGCTACCCCTTTAGGAGGGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatg                   <  1:5086064/150‑1 (MQ=255)
                                                                                                                                ggaGCGCCAGCCTGGCCGTTGACGCTACCCCTTTTGGTGGGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgat                 <  1:1547124/148‑1 (MQ=255)
                                                                                                                                gTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAcagg                                                                                                  >  1:5269031/1‑66 (MQ=14)
                                                                                                                                gTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAcagg                                                                                                  <  2:5269031/69‑4 (MQ=14)
                                                                                                                                gTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgat                 >  2:5181924/1‑150 (MQ=255)
                                                                                                                                 tAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatg                >  1:3508815/1‑150 (MQ=255)
                                                                                                                                   gcgcCAGCCTGGCTGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat              <  1:2760648/150‑1 (MQ=255)
                                                                                                                                   gcgcCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCCCAATCCACCCAGGCTCGACCGGCGACACCGGCGCCCCGGCGCCCatgatgag              >  2:2998599/1‑149 (MQ=38)
                                                                                                                                   gcgcCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat              >  2:749635/1‑150 (MQ=255)
                                                                                                                                   gcgcCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat              <  2:2678385/150‑1 (MQ=255)
                                                                                                                                    cgcCAGCCTGGCTGTTGCCGCTACCCGTTTTGGTGGGGTTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgatg             <  1:599664/150‑1 (MQ=255)
                                                                                                                                       cAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgatg             <  1:2734954/147‑1 (MQ=255)
                                                                                                                                       cAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgatg             >  2:2734954/1‑147 (MQ=255)
                                                                                                                                        aGCCTGGCATTTGACTCTACCCCTGTAGGAGCGGGTTTAACCGCGAAAGGGCCCGAACTGACCCTATAATGCTCTCAGTACTACCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGTCTCCCCATCACCAATGATCATGTGAt         >  2:4687775/1‑150 (MQ=255)
                                                                                                                                        aGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGAt         >  1:4838156/1‑150 (MQ=255)
                                                                                                                                         gCCTGGCAGTTGACGCTACCCCTGTTGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATc        <  2:2906064/150‑1 (MQ=255)
                                                                                                                                         gCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATc        >  2:687866/1‑150 (MQ=255)
                                                                                                                                          ccTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATcc       <  2:2046464/150‑1 (MQ=255)
                                                                                                                                           cTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCa      <  2:5556891/150‑1 (MQ=255)
                                                                                                                                            tGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAg     >  2:2347988/1‑150 (MQ=255)
                                                                                                                                            tGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGAGGTGATCCAg     >  1:5319291/1‑150 (MQ=255)
                                                                                                                                             ggCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAGc    >  1:1232081/1‑150 (MQ=255)
                                                                                                                                              gCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAGCa   >  2:1237159/1‑150 (MQ=255)
                                                                                                                                               cAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAGCAc  <  2:871150/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCAGCAC  >  NC_002947/6037656‑6037948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: