Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,504 A→T 6.2% intergenic (‑54/+102) PP_1662 ← / ← PP_1663 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,5040AT 6.2% 255.3 / 6.9 81intergenic (‑54/+102)PP_1662/PP_1663hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base T (2/3);  total (40/41)
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.43e-01

CCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAG  >  NC_002947/1859356‑1859650
                                                                                                                                                    |                                                                                                                                                  
ccTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACaa                                                                                                                                                   <  2:156425/150‑1 (MQ=255)
  ttCGTAGGAGATTCTCGAGTCGTCGCGTTCGGTTACTTTGACAATCAGCGCACACCTCCCCTACTCCACGCTGTTTGCGACAATGTATCGCACGTCTGCTGTATAGATTTTGGTAGCAGACTAAATAAGTAAAACAGGCATGAAACAAGt                                                                                                                                                 <  2:219546/150‑1 (MQ=255)
     gTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTc                                                                                                                                              >  2:249441/1‑150 (MQ=255)
        ggACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATg                                                                                                                                           <  2:1167642/150‑1 (MQ=255)
         gACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGtt                                                                                                                                                <  1:1438140/144‑1 (MQ=255)
         gACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGtt                                                                                                                                                >  2:1438140/1‑144 (MQ=255)
          aCATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGCACCCTCCACTTCTTCAAGATGTCTGCGAAAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGaa                                                                                                                                         <  1:130749/150‑1 (MQ=255)
               ctcGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCATgcctgttt                                                                                                                                    <  2:1839282/150‑8 (MQ=255)
                  gATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTTTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAg                                                                                                                                 <  2:671501/150‑1 (MQ=255)
                   aTTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAgg                                                                                                                                <  1:249441/150‑1 (MQ=255)
                     tcgtcgCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGa                                                                                                                              <  1:55521/150‑1 (MQ=255)
                       gtcgCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGAtg                                                                                                                            >  2:1793585/1‑150 (MQ=255)
                           cgTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGAtgctgc                                                                                                                        <  1:1266792/150‑1 (MQ=255)
                           cgTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGAtgctgc                                                                                                                        <  1:786768/150‑1 (MQ=255)
                                ggTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgg                                                                                                                    >  2:402135/1‑149 (MQ=255)
                                ggTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATAAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgg                                                                                                                    <  1:402135/149‑1 (MQ=255)
                                 gTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggt                                                                                                                  <  2:372230/150‑1 (MQ=255)
                                      ttGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCt                                                                                                             >  2:1896837/1‑150 (MQ=255)
                                      ttGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCGTTATAAAGTTTGTTAGCAGACTAAGTAAGTAAAACAAGCATGAAACAAGTTTCATGAAGCTCGAAGGCAGGCTGCCTGGGTGGGCt                                                                                                             >  1:1554021/1‑150 (MQ=255)
                                        ggACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAAAACGCATTAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCttt                                                                                                           <  1:1019639/150‑1 (MQ=255)
                                          aCAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGttt                                                                                                                                               <  1:363166/112‑1 (MQ=255)
                                          aCAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGttt                                                                                                                                               >  2:363166/1‑112 (MQ=255)
                                             gTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgc                                                                                                      >  2:650005/1‑150 (MQ=255)
                                              tCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACCGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGGCGGGGTGGGCTCTTTAgca                                                                                                     >  1:998017/1‑150 (MQ=255)
                                              tCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                     <  2:1306077/150‑1 (MQ=255)
                                                 tCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtgg                                                                                                                <  1:1829753/136‑1 (MQ=255)
                                                 tCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtgg                                                                                                                >  2:1829753/1‑136 (MQ=255)
                                                       aaCTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGaaa                                                                                            >  2:1788528/1‑150 (MQ=255)
                                                        aCTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGaaaa                                                                                           >  1:877074/1‑150 (MQ=255)
                                                         cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGAt                                                                                                                             >  2:862874/1‑115 (MQ=255)
                                                         cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGAt                                                                                                                             <  1:862874/115‑1 (MQ=255)
                                                         cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAg                                                                                          >  2:1169315/1‑150 (MQ=255)
                                                          tCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCATgc                                                                                                                                          <  1:101545/101‑2 (MQ=255)
                                                          tCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCATgc                                                                                                                                          >  2:101545/1‑100 (MQ=255)
                                                           ccACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGcc                                                                                        <  1:1711493/150‑1 (MQ=255)
                                                              cttcttCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCg                                                                                     >  1:1027891/1‑150 (MQ=255)
                                                              cttcttCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCg                                                                                     <  2:877074/150‑1 (MQ=255)
                                                                 cttcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCAc                                                                                  >  2:597262/1‑150 (MQ=255)
                                                                  ttcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACt                                                                                 >  1:1942394/1‑150 (MQ=255)
                                                                   tcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTg                                                                                <  2:1027891/150‑1 (MQ=255)
                                                                     aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAg                                                                                                                                        >  1:444884/1‑92 (MQ=255)
                                                                     aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAg                                                                                                                                        <  2:444884/92‑1 (MQ=255)
                                                                             tGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCCCACTGGGCGGGGCtt                                                                      >  2:1560267/1‑150 (MQ=255)
                                                                                gACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                   >  1:157164/1‑150 (MQ=255)
                                                                                gACAGTCTGTCGCAGGACTGATTTATCGAAGTTGTGAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAGAAGCCCCGCCCAGTGCGGGGCttttc                                                                   <  2:419354/150‑4 (MQ=255)
                                                                                  cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCGCGCTGTGGGCGGGGCTTCAGGt                                                                 >  2:931903/1‑150 (MQ=255)
                                                                                       tgtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgc                                                            <  1:62728/150‑1 (MQ=255)
                                                                                        gtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgca                                                           <  2:1942394/150‑1 (MQ=255)
                                                                                        gtCGCACGGCTTCTGTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgca                                                           <  1:273115/150‑1 (MQ=255)
                                                                                         tCGCACGTCTTCTTTATAGAGTTTGTTAGCAGGCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCttttt                                                                                                         <  2:853867/103‑1 (MQ=255)
                                                                                         tCGCACGTCTTCTTTATAGAGTTTGTTAGCAGGCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCttttt                                                                                                         >  1:853867/1‑103 (MQ=255)
                                                                                         tCGCACGTCTTCTTGATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcag                                                          >  2:1044045/1‑150 (MQ=255)
                                                                                                  ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACACGCCAGGTGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgcgc                                                 >  2:1441366/1‑150 (MQ=255)
                                                                                                  ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgg                                                                                                                    >  1:1711282/1‑83 (MQ=255)
                                                                                                  ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgg                                                                                                                    <  2:1711282/83‑1 (MQ=255)
                                                                                                    cttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGa                                               >  1:1687676/1‑150 (MQ=255)
                                                                                                    cttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGa                                               >  2:1029167/1‑150 (MQ=255)
                                                                                                     tttATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAGAAGCCCCGCCCAGTGCGGGGCttttctt                                                                 <  1:205734/131‑6 (MQ=255)
                                                                                                     tttATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAGAAGCCCCGCCCAGTGCGGGGCttttctt                                                                 >  2:205734/1‑126 (MQ=255)
                                                                                                        ataGAGTTTGTTAGCAGCCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAg                                           >  2:994577/1‑150 (MQ=255)
                                                                                                          agagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGtt                                         <  1:1896837/150‑1 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                       <  1:1044045/150‑1 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                       >  1:424043/1‑150 (MQ=255)
                                                                                                                 gTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggt                                                                                                                  <  2:1415458/70‑1 (MQ=255)
                                                                                                                 gTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggt                                                                                                                  >  1:1415458/1‑70 (MQ=255)
                                                                                                                  ttAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCt                                                                                                                      <  2:295806/65‑1 (MQ=255)
                                                                                                                  ttAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCt                                                                                                                      >  1:295806/1‑65 (MQ=255)
                                                                                                                             agtaagtaAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCa                      >  1:634797/1‑150 (MQ=255)
                                                                                                                                 agtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGa                  >  2:159942/1‑150 (MQ=255)
                                                                                                                                 agtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGa                  <  1:597262/150‑1 (MQ=255)
                                                                                                                                 agtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGAGCTTCAGGTATTGCAGGCAGTCCGCGAACAGCTGCCCTTCACCGGCCTGCCATCGa                  <  2:998017/150‑1 (MQ=255)
                                                                                                                                   taaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttctt                                                                 <  1:224009/101‑6 (MQ=38)
                                                                                                                                   taaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttctt                                                                 >  2:224009/1‑96 (MQ=38)
                                                                                                                                   taaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGAcc                <  2:157164/150‑1 (MQ=255)
                                                                                                                                     aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGt              <  1:1793585/150‑1 (MQ=255)
                                                                                                                                      aaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg             <  2:1687676/150‑1 (MQ=255)
                                                                                                                                         cttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                              >  1:1645619/4‑114 (MQ=255)
                                                                                                                                         cttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                              <  2:1645619/111‑1 (MQ=255)
                                                                                                                                               gAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgc    >  1:845883/1‑150 (MQ=255)
                                                                                                                                                 aaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTa                                                                >  1:885596/1‑88 (MQ=255)
                                                                                                                                                 aaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTa                                                                <  2:885596/88‑1 (MQ=255)
                                                                                                                                                 aaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcag  >  2:1700689/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAG  >  NC_002947/1859356‑1859650

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