Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,217 C→A 6.2% G192W (GGG→TGG)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,2170CA 6.2% 474.7 / 15.5 129G192W (GGG→TGG) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (68/53);  new base A (4/4);  total (72/57)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 135 total aligned reads displayed.

TGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTG  >  NC_000913/1580069‑1580366
                                                                                                                                                    |                                                                                                                                                     
tGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCg                                                                                                                                                      >  1:197537/1‑150 (MQ=255)
          aaaTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCa                                                                                                                                            <  1:1472571/150‑1 (MQ=255)
                cACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGa                                                                                                                                      <  2:1731074/150‑1 (MQ=255)
                cACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGa                                                                                                                                      <  2:1731073/150‑1 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:1024655/1‑150 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:646937/1‑150 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:1380591/1‑150 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:646934/1‑150 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:504826/1‑150 (MQ=255)
                   ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTAGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                   >  1:271170/1‑150 (MQ=255)
                    cacaTCGTTTTTAACCCATTTCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaaa                                                                                                                                  <  2:2164017/150‑1 (MQ=255)
                    cacaTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaaa                                                                                                                                  <  2:1277478/150‑1 (MQ=255)
                      caTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAAtt                                                                                                                                <  2:192926/150‑1 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  2:845172/1‑150 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  1:1739720/1‑150 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  2:318919/1‑150 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  2:569802/1‑150 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  2:1659820/1‑150 (MQ=255)
                        tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc                                                                                                                              >  2:1000075/1‑150 (MQ=255)
                         gttttttAACACAGTGCATAAAAATGGTTTAAATAAACTTGAGATTATCCTTGTTATGGACAGAAAAATCAATGATCCTTAATGTGTTCTCACTATGACCACTGAAATCTATATTCGGAGTGCCTGTTTATAGAAATTCGTTAAATTGCa                                                                                                                             <  2:1427662/148‑1 (MQ=255)
                         cGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCa                                                                                                                             <  2:1170847/150‑1 (MQ=255)
                         cGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCa                                                                                                                             <  2:2478837/150‑1 (MQ=255)
                                  cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTAc                                                                                                                    >  2:19251/1‑150 (MQ=255)
                                  cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTAc                                                                                                                    >  2:1358640/1‑150 (MQ=255)
                                  cccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTAc                                                                                                                    >  2:585144/1‑150 (MQ=255)
                                   ccATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACt                                                                                                                   >  1:2241707/1‑150 (MQ=255)
                                      ttGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGcc                                                                                                                <  1:1454003/150‑1 (MQ=255)
                                      ttGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGcc                                                                                                                <  1:802919/150‑1 (MQ=255)
                                        gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGa                                                                                                              <  1:2488971/150‑1 (MQ=255)
                                        gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGa                                                                                                              <  2:1971022/150‑1 (MQ=255)
                                           tAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAg                                                                                                           >  1:66339/1‑150 (MQ=255)
                                              aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGAtt                                                                                                        >  2:296685/1‑150 (MQ=255)
                                              aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGAtt                                                                                                        >  2:296436/1‑150 (MQ=255)
                                                    ttGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaa                                                                                                  >  1:1054481/1‑150 (MQ=255)
                                                    ttGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaa                                                                                                  >  2:648392/1‑150 (MQ=255)
                                                    ttGACATAAACTTGCCATAAGCCGTAGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCCCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaa                                                                                                  >  1:1054504/1‑150 (MQ=255)
                                                     tGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaaa                                                                                                 >  2:1350919/1‑150 (MQ=255)
                                                     tGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaaa                                                                                                 >  2:2625982/1‑150 (MQ=255)
                                                     tGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTaaaaa                                                                                                 >  2:1427037/1‑150 (MQ=255)
                                                         aTACACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGa                                                                                             <  2:2098894/150‑1 (MQ=255)
                                                         aTAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGa                                                                                             <  2:82476/150‑1 (MQ=255)
                                                         aTAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGa                                                                                             <  2:1463451/150‑1 (MQ=255)
                                                         aTAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGa                                                                                             <  2:181972/150‑1 (MQ=255)
                                                         aTAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGAATATGTTTACTGCCGATAGATTTTAAAAAATGa                                                                                             <  2:307774/150‑1 (MQ=255)
                                                          tAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGAt                                                                                            <  2:1622172/150‑1 (MQ=255)
                                                           aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATa                                                                                           <  2:1936383/150‑1 (MQ=255)
                                                           aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATa                                                                                           <  2:2184605/150‑1 (MQ=255)
                                                           aaaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATa                                                                                           <  2:1170717/150‑1 (MQ=255)
                                                            aaCTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATaa                                                                                          <  2:2379095/150‑1 (MQ=255)
                                                               ttGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAAcc                                                                                       >  1:1613182/1‑150 (MQ=255)
                                                               ttGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAAcc                                                                                       >  1:1994898/1‑150 (MQ=255)
                                                               ttGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAAcc                                                                                       >  1:1109639/1‑150 (MQ=255)
                                                                  ccATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCa                                                                                    >  1:172875/1‑150 (MQ=255)
                                                                       aGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAgg                                                                               >  2:389291/1‑150 (MQ=255)
                                                                       aGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAgg                                                                               >  1:2260638/1‑150 (MQ=255)
                                                                          cGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTa                                                                            >  1:1148394/1‑150 (MQ=255)
                                                                          cGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTa                                                                            >  1:2514054/1‑150 (MQ=255)
                                                                                 ggCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                     <  1:1016485/150‑1 (MQ=255)
                                                                                 ggCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                     <  1:1016384/150‑1 (MQ=255)
                                                                                 ggCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                     <  1:1538479/150‑1 (MQ=255)
                                                                                 ggCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                     <  1:202810/150‑1 (MQ=255)
                                                                                 ggCACAGAAAAATAAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACa                                                                     <  1:1495933/150‑1 (MQ=255)
                                                                                  gCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACAt                                                                    >  2:2027448/1‑150 (MQ=255)
                                                                                  gCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACAt                                                                    >  2:615876/1‑150 (MQ=255)
                                                                                          aaataaatGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGtta                                                            <  2:566278/150‑1 (MQ=255)
                                                                                          aaaTCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGtta                                                            <  2:566269/150‑1 (MQ=255)
                                                                                                        gAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTaa                                              >  1:848452/1‑150 (MQ=255)
                                                                                                         aaTGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAAc                                             <  1:1418213/150‑1 (MQ=255)
                                                                                                               ttCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTa                                       <  1:523579/150‑1 (MQ=255)
                                                                                                                 ctcACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaaa                                     <  1:2220684/150‑1 (MQ=255)
                                                                                                                 ctcACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaaa                                     <  2:197537/150‑1 (MQ=255)
                                                                                                                  tcACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGt                                                 >  2:1145463/1‑137 (MQ=255)
                                                                                                                  tcACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGt                                                 <  1:1145463/137‑1 (MQ=255)
                                                                                                                            cACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGat                          >  1:2614498/1‑150 (MQ=255)
                                                                                                                            cACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGat                          >  1:2614499/1‑150 (MQ=255)
                                                                                                                            cACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGat                          >  2:1531371/1‑150 (MQ=255)
                                                                                                                            cACTGAAATCAAAAATCGGCGACCCCGTTACTAGCAATTCGATAAATTGCATAAGTTAACAGCCTAGAGATATAAAAAAATGATAAACCTCAAGAGGGTAATGGAAAGTAACGAAATTAATTACTGTAAAAGAATAAAACTCTAATTGaa                          >  2:1597541/1‑149 (MQ=255)
                                                                                                                             aCTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGata                         <  1:1966723/150‑1 (MQ=255)
                                                                                                                             aCTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGCATTAAACTCTACTTGata                         <  1:216710/150‑1 (MQ=255)
                                                                                                                                   aTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                   >  1:2671080/1‑150 (MQ=255)
                                                                                                                                   aTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                   >  1:2323686/1‑150 (MQ=255)
                                                                                                                                   aTCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAttt                   >  1:166861/1‑150 (MQ=255)
                                                                                                                                    tCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTACTGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                  <  2:1994898/150‑1 (MQ=255)
                                                                                                                                    tCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                  <  2:1613182/150‑1 (MQ=255)
                                                                                                                                    tCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                  <  2:1109639/150‑1 (MQ=255)
                                                                                                                                    tCAATATTCGGCGTCCAGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                  <  2:33102/150‑1 (MQ=255)
                                                                                                                                          ttCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGc            >  1:2267927/1‑150 (MQ=255)
                                                                                                                                          ttCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGc            >  2:1403781/1‑150 (MQ=255)
                                                                                                                                          ttCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGc            >  2:1807044/1‑150 (MQ=255)
                                                                                                                                          ttCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGc            >  2:1938711/1‑150 (MQ=255)
                                                                                                                                            cGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc          >  1:2119091/1‑150 (MQ=255)
                                                                                                                                            cGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc          >  1:2417371/1‑150 (MQ=255)
                                                                                                                                            cGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc          >  1:459657/1‑150 (MQ=255)
                                                                                                                                            cGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc          >  1:2417666/1‑150 (MQ=255)
                                                                                                                                            cGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACAGTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc          >  1:1218040/1‑150 (MQ=255)
                                                                                                                                                 tCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATgg     <  1:2211462/150‑1 (MQ=255)
                                                                                                                                                   ccGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgct   >  2:1863094/1‑150 (MQ=255)
                                                                                                                                                   ccGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgct   >  1:1394940/1‑150 (MQ=255)
                                                                                                                                                    cGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgctg  >  2:1065331/1‑150 (MQ=255)
                                                                                                                                                    cGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGgctg  >  2:1065313/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTG  >  NC_000913/1580069‑1580366

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