Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,000 G→T 6.0% intergenic (+46/‑151) PP_3009 → / → PP_3010 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,0000GT 6.0% 216.9 / 2.9 67intergenic (+46/‑151)PP_3009/PP_3010hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base T (2/2);  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 = 6.16e-01

ATGACTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGT  >  NC_002947/3401856‑3402147
                                                                                                                                                |                                                                                                                                                   
aTGACTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGa                                                                                                                                                >  2:1711036/1‑150 (MQ=255)
 tGACTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCAGCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAg                                                                                                                                               <  2:1466088/150‑1 (MQ=255)
   aCTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAggg                                                                                                                                             >  1:1338461/1‑150 (MQ=255)
    cTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGa                                                                                                                                            <  1:89698/150‑1 (MQ=255)
       aGGATGACCGCTATTACATGATTTTCAGCGTGGTCGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCCCCCAGCAACCACCACACCGAGGTGACCTTGTCTCGTTTGTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAAc                                                                                                                                         <  1:35985/150‑1 (MQ=255)
            gACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGtccc                                                                                                                                    <  1:1212350/150‑1 (MQ=255)
               cGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGccc                                                                                                                                     <  1:144816/146‑1 (MQ=255)
               cGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGccc                                                                                                                                     >  2:144816/1‑146 (MQ=255)
                    ttACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCta                                                                                                                            <  2:1762395/150‑1 (MQ=255)
                                 aGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGcaaacaaa                                                                                                               >  1:1714076/1‑150 (MQ=255)
                                  gCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                              <  2:1872221/150‑1 (MQ=255)
                                        ttGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATg                                                                                                        >  2:1940152/1‑150 (MQ=255)
                                           ctgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGccc                                                                                                     <  1:1940152/150‑1 (MQ=255)
                                            tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                    >  2:1116278/1‑150 (MQ=255)
                                            tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                    >  2:1904038/1‑150 (MQ=255)
                                            tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                    >  1:334271/1‑150 (MQ=255)
                                                ctgGTACTCGATCTGGCGGCGACTTATGTGGTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgttg                                                                                                >  1:1810001/1‑150 (MQ=255)
                                                 tgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGa                                                                                                                                                <  2:1045980/101‑1 (MQ=255)
                                                 tgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGa                                                                                                                                                >  1:1045980/1‑101 (MQ=255)
                                                    tACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGt                                                                                            <  1:1553529/150‑1 (MQ=255)
                                                           tCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAg                                                                                     >  1:95004/1‑150 (MQ=255)
                                                                 ggcgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACAt                                                                               >  1:1899835/1‑150 (MQ=255)
                                                                 gacgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACAt                                                                               >  1:125284/3‑150 (MQ=255)
                                                                  gcgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATg                                                                              <  2:334271/150‑1 (MQ=255)
                                                                      cTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                          <  2:1338461/150‑1 (MQ=255)
                                                                       ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                         >  2:1811033/1‑150 (MQ=255)
                                                                       ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                         >  1:297177/1‑150 (MQ=255)
                                                                       ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                         <  1:1711036/150‑1 (MQ=255)
                                                                          tgtgTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgtt                                                                                                 >  2:1594708/1‑123 (MQ=255)
                                                                          tgtgTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgtt                                                                                                 <  1:1594708/123‑1 (MQ=255)
                                                                            tgtTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACaa                                                                    >  1:791847/1‑150 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCtatata                                                                                                                        >  2:916230/1‑92 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCtatata                                                                                                                        <  1:916230/92‑1 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                              >  1:1903040/1‑102 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                              <  2:1903040/102‑1 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                          <  1:1228159/138‑1 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                          >  2:1228159/1‑138 (MQ=255)
                                                                                  tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCa                                                              >  2:1032348/1‑150 (MQ=255)
                                                                                     aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTCGt                                                                                            <  1:1937478/117‑1 (MQ=255)
                                                                                     aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGt                                                                                            >  2:1937478/1‑117 (MQ=255)
                                                                                     aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaa                                                           >  1:1488950/1‑150 (MQ=255)
                                                                                      cGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaa                                                          <  1:1811033/150‑1 (MQ=255)
                                                                                        ccATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAg                                                        <  2:125284/150‑1 (MQ=255)
                                                                                        ccATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAg                                                        <  2:1488950/150‑1 (MQ=255)
                                                                                         cATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAgg                                                       >  2:1322237/1‑150 (MQ=255)
                                                                                          aTCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAggg                                                      <  1:1904038/150‑1 (MQ=255)
                                                                                          aTCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAggg                                                      >  2:1896442/1‑150 (MQ=255)
                                                                                             cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                   >  1:1252357/1‑150 (MQ=255)
                                                                                             cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                   >  1:1252714/1‑150 (MQ=255)
                                                                                             cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaaa                                                         >  2:1830572/1‑144 (MQ=255)
                                                                                             cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaaa                                                         <  1:1830572/144‑1 (MQ=255)
                                                                                             cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                   <  2:1714076/150‑1 (MQ=255)
                                                                                              ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCa                                                                                      >  1:1458017/1‑114 (MQ=255)
                                                                                              ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCa                                                                                      <  2:1458017/114‑1 (MQ=255)
                                                                                              ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAg                                                  >  2:412966/1‑150 (MQ=255)
                                                                                                    cccGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGgctgct                                            >  2:216400/1‑150 (MQ=255)
                                                                                                      cGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTg                                          <  1:1322237/150‑1 (MQ=255)
                                                                                                         tAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCa                                       >  2:1040955/1‑150 (MQ=255)
                                                                                                                  ggCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTccc                              >  1:1745878/1‑150 (MQ=255)
                                                                                                                       tACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGttt                         <  1:412966/150‑1 (MQ=255)
                                                                                                                        aCTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                        <  2:1252714/150‑1 (MQ=255)
                                                                                                                        aCTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                        <  2:1252357/150‑1 (MQ=255)
                                                                                                                        aCTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                        <  1:216400/150‑1 (MQ=255)
                                                                                                                        aCTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAAAAAAAAAAAAAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                        <  2:1810001/150‑1 (MQ=255)
                                                                                                                              aaTCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATgg                  <  1:1040955/150‑1 (MQ=255)
                                                                                                                                             aGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTgg                                         >  1:663805/1‑112 (MQ=255)
                                                                                                                                             aGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTgg                                         <  2:663805/112‑1 (MQ=255)
                                                                                                                                              gggCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCgtgt  >  2:1417036/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATGACTCAGGATGACCGCTATTACATGATTTTCAGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGT  >  NC_002947/3401856‑3402147

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