Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,584,863 T→A 7.2% K542* (AAA→TAA)  ydeP ← putative oxidoreductase YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,8630TA 7.2% 253.5 / 6.8 82K542* (AAA→TAA) ydePputative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base A (3/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAG  >  NC_000913/1584716‑1585010
                                                                                                                                                   |                                                                                                                                                   
aCTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTa                                                                                                                                                   <  2:1241374/150‑1 (MQ=255)
aCTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTa                                                                                                                                                   >  1:1127826/1‑150 (MQ=255)
  tctGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAac                                                                                                                                                 <  2:116598/150‑1 (MQ=255)
     ggCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAAcacg                                                                                                                                              <  1:290080/150‑1 (MQ=255)
     ggCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGAGTTCAGCATTACACCGGCGGGTTTTAAcacg                                                                                                                                              >  2:1281263/1‑150 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATTCGATCATAATCTTCCACCAGATACTACCAGGCTACCACGCTCTGGGGTAGTGATGCCTCCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATAACACCGGCGGGTTTTAAcaagccacg                                                                                                                                         >  2:1436311/1‑150 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                         <  1:935903/150‑4 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                         <  1:407355/150‑4 (MQ=255)
            caGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCAcgcg                                                                                                                                       <  2:267611/150‑1 (MQ=255)
            caGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCAcgcg                                                                                                                                       >  2:408210/1‑150 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                     >  1:548957/1‑150 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                     <  2:90742/150‑1 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                     <  1:199873/150‑1 (MQ=255)
                  cAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCAt                                                                                                                                 <  2:213236/150‑1 (MQ=255)
                    aTGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGa                                                                                                                               >  1:771630/1‑150 (MQ=255)
                      gTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTTTAACACGCCACGCGAGGCATGAat                                                                                                                             <  2:1577221/150‑1 (MQ=255)
                         aTTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcat                                                                                                                          >  1:1419449/1‑150 (MQ=255)
                            gCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGa                                                                                                                       >  2:78805/1‑150 (MQ=255)
                                 tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                                                                                                  <  1:95306/150‑1 (MQ=255)
                                 tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTCCGCGATTCCCCCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                                                                                                  >  2:221269/1‑150 (MQ=255)
                                   cGATCATATTCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                                                                                                >  2:1613516/1‑150 (MQ=255)
                                   cGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                                                                                                >  2:1057693/1‑150 (MQ=255)
                                      tCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATcc                                                                                                             <  2:288147/150‑1 (MQ=255)
                                                ccaccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTa                                                                                                   <  1:1302591/150‑1 (MQ=255)
                                                     aGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGcc                                                                                              <  1:951857/150‑1 (MQ=255)
                                                        tACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgc                                                                                           >  1:1071427/1‑150 (MQ=255)
                                                         aCTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcg                                                                                          <  1:548863/150‑1 (MQ=255)
                                                          cTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcgc                                                                                         <  2:1419449/150‑1 (MQ=255)
                                                            cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc                                                                                       >  2:590882/1‑150 (MQ=255)
                                                            cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc                                                                                       <  2:125926/150‑1 (MQ=255)
                                                                 gCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGttt                                                                                  <  1:446459/150‑1 (MQ=255)
                                                                  cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGtttt                                                                                 <  1:711298/150‑1 (MQ=255)
                                                                  cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGATACCGCCTGCGCACCGtttt                                                                                 >  2:847490/1‑150 (MQ=255)
                                                                    aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGtttttt                                                                               <  1:499670/150‑1 (MQ=255)
                                                                       aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGc                                                                            >  2:672586/1‑150 (MQ=255)
                                                                       aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGc                                                                            >  2:60325/1‑150 (MQ=255)
                                                                          ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAATGAATCCTCAACGGTTACCGCCTGCGCACAGTTTTTTTGCATg                                                                         <  2:644258/150‑1 (MQ=255)
                                                                                ggTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAAtc                                                                   >  1:106240/1‑150 (MQ=255)
                                                                                   aGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCTCACCGTTTTTTTGCATGTCAATCTCg                                                                <  2:771630/150‑1 (MQ=255)
                                                                                      gctgcCTGCGCGATTCCCGCGACCACTTCACACTCTGATTTCAGCATTACACCGGCGTGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTa                                                             <  1:844177/150‑1 (MQ=255)
                                                                                      gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTa                                                             >  1:1525379/1‑150 (MQ=255)
                                                                                        tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACg                                                           >  2:227816/1‑150 (MQ=255)
                                                                                        tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACg                                                           >  1:471222/1‑150 (MQ=255)
                                                                                            tgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTccc                                                       >  1:774349/1‑150 (MQ=255)
                                                                                               gccaTTCCCTCGACCACTGCACACTATGATTTAAGCATTACACAGGCGGGTTTTAACACGCAACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAgg                                                    >  2:594238/4‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                 <  2:106240/150‑1 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                 >  1:579706/1‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                 >  2:1304759/1‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                 >  1:1087865/1‑150 (MQ=255)
                                                                                                    tCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                               >  1:1649977/1‑150 (MQ=255)
                                                                                                     cccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGc                                              <  1:507316/150‑1 (MQ=255)
                                                                                                      ccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCa                                             >  2:1572589/1‑150 (MQ=255)
                                                                                                      ccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCa                                             >  2:1566771/1‑150 (MQ=255)
                                                                                                         cgACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                          >  2:411973/1‑150 (MQ=255)
                                                                                                             cACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAAtata                                      <  1:1057693/150‑1 (MQ=255)
                                                                                                              aCTGCACACAGTGATTTAAGCATTACACCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGGCATTGCATCCTCAACGGTTTCCGCCTGCGCACCGTTTTTTTTCATGTGAATCACGCTACGTCCCAGGCCCGGCAGAATATAg                                     >  2:1027800/1‑150 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCTGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                    <  1:1456546/150‑1 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                    >  2:1493932/1‑150 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                    >  1:222341/1‑150 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAg                                                     >  2:55173/1‑131 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAg                                                     <  1:55173/131‑1 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTa                                  <  2:343516/150‑1 (MQ=255)
                                                                                                                   acacTGTGATTTCAGCATTACACCGGCGGGTTATAACACGCCACGCGAGGCATGAATCATCGTCATTGAATCCTCATCTGTTAACGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCGAGGAGCGGCAGAATATAGCTAtt                                <  2:1577258/150‑2 (MQ=255)
                                                                                                                   acacTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATg                                <  1:408210/150‑1 (MQ=255)
                                                                                                                        ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCaa                                                                     <  1:421939/108‑1 (MQ=255)
                                                                                                                        ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCaa                                                                     >  2:421939/1‑108 (MQ=255)
                                                                                                                           aTTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCgg                        <  1:221269/150‑1 (MQ=255)
                                                                                                                             ttCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTc                      >  2:1175049/1‑150 (MQ=255)
                                                                                                                                aGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAAc                   <  1:227816/150‑1 (MQ=255)
                                                                                                                                 gCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACa                  >  1:1316812/1‑150 (MQ=255)
                                                                                                                                  cATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAg                 <  1:636520/150‑1 (MQ=255)
                                                                                                                                   aTTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGa                >  1:581204/1‑150 (MQ=255)
                                                                                                                                    ttACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt               <  1:1199105/150‑1 (MQ=255)
                                                                                                                                     tACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg              <  1:672586/150‑1 (MQ=255)
                                                                                                                                       gcgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag            >  2:1513479/5‑150 (MQ=255)
                                                                                                                                       gcgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag            >  1:1299214/5‑150 (MQ=255)
                                                                                                                                       cacCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag            <  1:219132/150‑1 (MQ=255)
                                                                                                                                        cgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCATCAACGGTTACCGCCTGCGCACCGTTTTTTTGCAAGTCAATCTCGCTACGACCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAAgaga           >  2:911987/4‑150 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg              >  1:899964/1‑146 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg              <  2:899964/146‑1 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCAAGAATCATCGACATTGAATCCTCAACGGTTACCGACTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATCCCGTGCGGTAAACAGATGagag          >  2:1173414/1‑150 (MQ=255)
                                                                                                                                               ggTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTa    <  1:1436311/150‑1 (MQ=255)
                                                                                                                                                 ttttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAg  <  2:1649977/150‑1 (MQ=255)
                                                                                                                                                 ttttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAg  <  2:1279421/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAG  >  NC_000913/1584716‑1585010

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