Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,588 A→T 5.4% intergenic (‑53/+16) guaA ← / ← guaB GMP synthetase (glutamine aminotransferase)/IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5880AT 5.4% 424.1 / 10.2 121intergenic (‑53/+16)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (49/65);  new base T (4/3);  total (53/68)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 121 total aligned reads displayed.

AACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCC  >  NC_000913/2632440‑2632735
                                                                                                                                                    |                                                                                                                                                   
 acacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGt                                                                                                                                                   >  2:134014/1‑150 (MQ=255)
 acacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGt                                                                                                                                                   >  2:134016/1‑150 (MQ=255)
   acCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCg                                                                                                                                                 >  1:1217717/1‑150 (MQ=255)
   acCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCg                                                                                                                                                 >  1:351532/1‑150 (MQ=255)
     ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCTAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg                                                                                                                                               <  1:777247/150‑1 (MQ=255)
     ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCggg                                                                                                                                               <  2:709840/150‑1 (MQ=255)
     ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCggg                                                                                                                                               <  2:1368981/150‑1 (MQ=255)
        gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGa                                                                                                                                            >  1:1370143/1‑150 (MQ=255)
            acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAagag                                                                                                                                        >  2:337528/1‑150 (MQ=255)
            acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAagag                                                                                                                                        >  2:1158967/1‑150 (MQ=255)
               cacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                                     <  2:251708/150‑1 (MQ=255)
               cacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                                     <  2:251726/150‑1 (MQ=255)
               cacgcGGCGCGCAACCACTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATCAATGATTACCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAAAATAAATCGCCCGACATGAAGTCGGGCGAAGAGAAt                                                                                                                                     <  2:852007/150‑1 (MQ=255)
                acgcGGCGCGCAACCTTTTGTGTTTACTGAGAACCGAAGTCCATATTTAGGATGCGATGCTTATGAATGTTTTCCGTCTTTGACGCTTATTCCGAGGCAAGTTAAACATATAATATAAATCGCCCGACATGATGTCGGGCGAAGAGAATc                                                                                                                                    <  1:1725385/150‑1 (MQ=255)
                 cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                                                                                   <  2:70353/150‑1 (MQ=255)
                 cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                                                                                   <  2:1362956/150‑1 (MQ=255)
                   cggcgCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAgg                                                                                                                                 <  1:196413/150‑1 (MQ=255)
                    ggcgCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                                                                                <  2:839290/150‑1 (MQ=255)
                             ccAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                                                                       >  2:1234824/1‑150 (MQ=255)
                                   gAGTGTACTGAGAACCGAAGTCCAGACTGAGGATGCGATGCTTATGACTGTTTTCCGTCATGGACGCTTATTCCGAGGCCAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGtt                                                                                                                 <  2:513458/150‑1 (MQ=255)
                                   gAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGtt                                                                                                                 <  2:513484/150‑1 (MQ=255)
                                       gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  2:1613506/150‑1 (MQ=255)
                                       gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  1:1623571/150‑1 (MQ=255)
                                       gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  1:153678/150‑1 (MQ=255)
                                       gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  1:153538/150‑1 (MQ=255)
                                       gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  1:1417504/150‑1 (MQ=255)
                                        tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                                                                            >  1:875370/1‑150 (MQ=255)
                                        tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                                                                            >  1:601004/1‑150 (MQ=255)
                                             agaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActct                                                                                                       >  2:545340/1‑150 (MQ=255)
                                              gaACCGAAGTCCCGACTGAGGATGCGATGCTTATGACTGTTTTCCGTCATTGACGCTTATTCCGCGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                                      <  1:218035/150‑1 (MQ=255)
                                                aCCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTa                                                                                                    <  2:1227816/150‑1 (MQ=255)
                                                     aGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                                                                               <  2:1339027/150‑1 (MQ=255)
                                                     aGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                                                                               <  2:1371197/150‑1 (MQ=255)
                                                      gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                              <  1:1567094/150‑1 (MQ=255)
                                                      gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                              <  1:1518399/150‑1 (MQ=255)
                                                          aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCa                                                                                          >  2:828768/1‑150 (MQ=255)
                                                          aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCa                                                                                          <  1:447817/150‑1 (MQ=255)
                                                           gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAc                                                                                         >  2:248737/1‑150 (MQ=255)
                                                           gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAc                                                                                         >  2:40575/1‑150 (MQ=255)
                                                              tGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtc                                                                                      <  1:236711/150‑1 (MQ=255)
                                                                aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                                                                    >  1:1489398/1‑150 (MQ=255)
                                                                aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                                                                    <  1:1585999/150‑1 (MQ=255)
                                                                  gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGa                                                                                  <  1:634313/150‑1 (MQ=255)
                                                                  gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGa                                                                                  >  1:1603722/1‑150 (MQ=255)
                                                                    tgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAAc                                                                                <  2:351532/150‑1 (MQ=255)
                                                                      cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTACTGGTCACGTCGTGAACGt                                                                              <  1:1039514/150‑1 (MQ=255)
                                                                      cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                                                                              <  1:409928/150‑1 (MQ=255)
                                                                      cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                                                                              <  1:1608093/150‑1 (MQ=255)
                                                                      cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                                                                              <  1:1608091/150‑1 (MQ=255)
                                                                       gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                                                                             >  2:1429266/1‑150 (MQ=255)
                                                                           cTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCtt                                                                         >  1:308662/1‑150 (MQ=255)
                                                                             tATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATACATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                                                                       <  1:368655/150‑1 (MQ=255)
                                                                              aTGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCTTGAACGTGGCTTTcc                                                                      <  2:1465042/150‑1 (MQ=255)
                                                                                 aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGa                                                                   >  1:1710119/1‑150 (MQ=255)
                                                                                        tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                            <  1:78277/150‑1 (MQ=255)
                                                                                          cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGATTGCCCGct                                                          >  2:665179/1‑149 (MQ=255)
                                                                                          cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCa                                                          >  1:1709385/1‑150 (MQ=255)
                                                                                          cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCa                                                          >  2:527620/1‑150 (MQ=255)
                                                                                            tCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAcc                                                        >  2:1625183/1‑150 (MQ=255)
                                                                                               ttGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                                     >  1:1016283/1‑150 (MQ=255)
                                                                                               ttGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                                     >  1:1508652/1‑150 (MQ=255)
                                                                                                tGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTg                                                    >  2:1252311/1‑150 (MQ=255)
                                                                                                    gCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATAc                                                <  2:826943/150‑1 (MQ=255)
                                                                                                    gCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATAc                                                <  2:827106/150‑1 (MQ=255)
                                                                                                      ttATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGt                                              >  1:925397/1‑150 (MQ=255)
                                                                                                      ttATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGt                                              >  2:880620/1‑150 (MQ=255)
                                                                                                       tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTa                                             >  1:258222/1‑150 (MQ=255)
                                                                                                                   aGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                                 >  1:1036142/1‑150 (MQ=255)
                                                                                                                   aGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                                 >  1:494094/1‑150 (MQ=255)
                                                                                                                     tGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                               >  2:356398/1‑150 (MQ=255)
                                                                                                                       aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGTACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             <  2:1650558/150‑1 (MQ=255)
                                                                                                                       aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             <  2:609166/150‑1 (MQ=255)
                                                                                                                       aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             <  2:571609/150‑1 (MQ=255)
                                                                                                                       aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGACTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             <  1:69179/150‑1 (MQ=255)
                                                                                                                         aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTAc                           <  2:523628/150‑1 (MQ=255)
                                                                                                                           aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACg                          <  2:225240/149‑1 (MQ=255)
                                                                                                                           aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACg                          >  1:225240/1‑149 (MQ=255)
                                                                                                                            gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCa                        <  2:875370/150‑1 (MQ=255)
                                                                                                                            gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCa                        <  2:337304/150‑1 (MQ=255)
                                                                                                                            gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCa                        <  2:337291/150‑1 (MQ=255)
                                                                                                                              taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGt                      <  2:623335/150‑1 (MQ=255)
                                                                                                                                  ataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgt                  >  1:1162773/1‑150 (MQ=255)
                                                                                                                                  ataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgt                  <  1:200859/150‑1 (MQ=255)
                                                                                                                                          cccGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          >  2:875022/1‑150 (MQ=255)
                                                                                                                                          cccGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          >  2:1505241/1‑150 (MQ=255)
                                                                                                                                          cccGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          >  2:898914/1‑150 (MQ=255)
                                                                                                                                          cccGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTGTAGTACGCAGTTCGTCGATAGTa          >  2:268342/1‑150 (MQ=255)
                                                                                                                                          cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          <  2:1724520/150‑1 (MQ=255)
                                                                                                                                          cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          <  2:543711/150‑1 (MQ=255)
                                                                                                                                          cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          <  2:543642/150‑1 (MQ=255)
                                                                                                                                          cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa          >  2:1163152/1‑150 (MQ=255)
                                                                                                                                          cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCGTTAGTACGCAGTTCGTCGAGAGTa          >  2:1163205/1‑150 (MQ=255)
                                                                                                                                           ccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAc         >  1:73358/1‑150 (MQ=255)
                                                                                                                                                aTGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAg    >  1:1490456/1‑150 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCTGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGcc  <  2:1335357/150‑1 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                            <  1:1424362/60‑1 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                            >  2:1424362/1‑60 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGcc  <  2:1034876/150‑1 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGcc  <  2:1045072/150‑1 (MQ=255)
                                                                                                                                                  gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGcc  <  2:170767/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCC  >  NC_000913/2632440‑2632735

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: