Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 28,399 C→G 5.5% A9G (GCC→GGC)  dapB → 4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,3990CG 5.5% 310.7 / 6.7 91A9G (GCC→GGC) dapB4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (45/41);  new base G (3/2);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAG  >  NC_000913/28251‑28548
                                                                                                                                                    |                                                                                                                                                     
cTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGcc                                                                                                                                                      >  2:1322978/1‑150 (MQ=255)
   tCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATc                                                                                                                                                   >  1:882918/1‑150 (MQ=255)
      cTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATcgcg                                                                                                                                                <  2:372250/150‑1 (MQ=255)
      cTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCAAGATGCAAACATCCGCGTTGCCATcgcg                                                                                                                                                >  1:549420/1‑150 (MQ=255)
        gATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCggg                                                                                                                                              >  1:798136/1‑150 (MQ=255)
        gATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGATTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCggg                                                                                                                                              <  2:233572/150‑1 (MQ=255)
          tatCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAg                                                                                                                                            <  2:1048221/150‑1 (MQ=255)
          tatCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGATGGCAACGCggatg                                                                                                                                            <  1:728149/150‑4 (MQ=255)
             ccctccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCg                                                                                                                                         >  2:300282/1‑150 (MQ=255)
              cctccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCgg                                                                                                                                        >  1:801869/1‑150 (MQ=255)
              cctccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCgg                                                                                                                                        >  1:842602/1‑150 (MQ=255)
              cctccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCGGGGTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGATGGCAACGCggatgtttg                                                                                                                                        <  2:661053/150‑8 (MQ=255)
                 ccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCggggg                                                                                                                                     >  2:255354/1‑150 (MQ=255)
                  cctGCTTATCATTAATTTCTAATTATCAGCCTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTCACTCTGAAAACGGGCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAAGATGCGCGTTGCCATCGCGGGAGCCgggggg                                                                                                                                    <  1:914075/150‑1 (MQ=255)
                    tGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCg                                                                                                                                  <  2:437934/150‑1 (MQ=255)
                        tATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATg                                                                                                                              >  2:1454913/1‑150 (MQ=255)
                             ttAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGgccg                                                                                                                         >  1:473578/1‑150 (MQ=255)
                               aaTTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGgccgcc                                                                                                                       <  2:1451865/150‑1 (MQ=255)
                               aaTTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGgccgcc                                                                                                                       <  1:630462/150‑1 (MQ=255)
                                   tCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGtt                                                                                                                   <  2:38184/150‑1 (MQ=255)
                                     tAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGa                                                                                                                 <  1:255354/150‑1 (MQ=255)
                                        ttATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAGCGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTc                                                                                                              <  2:829055/150‑1 (MQ=255)
                                          aTCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAg                                                                                                            >  2:808600/1‑150 (MQ=255)
                                            cAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAggc                                                                                                          >  2:350991/1‑150 (MQ=255)
                                            cAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAggc                                                                                                          >  2:1086606/1‑150 (MQ=255)
                                               cGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGcc                                                                                                                                                      >  2:863426/1‑103 (MQ=255)
                                               cGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGcc                                                                                                                                                      <  1:863426/103‑1 (MQ=255)
                                                         tggcGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTa                                                                                             <  1:350991/150‑1 (MQ=255)
                                                            cggcgTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGCTGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTagag                                                                                          >  2:635385/1‑150 (MQ=255)
                                                                 tAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGt                                                                                     >  2:232042/1‑150 (MQ=255)
                                                                 tAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGt                                                                                     <  2:902690/150‑1 (MQ=255)
                                                                    cGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCa                                                                                  >  1:53014/1‑150 (MQ=255)
                                                                     gATGCACTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAg                                                                                 >  1:526227/1‑150 (MQ=255)
                                                                       tGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGtt                                                                               <  2:798136/150‑1 (MQ=255)
                                                                       tGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGtt                                                                               >  1:1093697/1‑150 (MQ=255)
                                                                        gcgcTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTg                                                                              <  2:9830/150‑1 (MQ=255)
                                                                           cTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgc                                                                           <  2:1535735/150‑1 (MQ=255)
                                                                           cTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgc                                                                           <  1:1454913/150‑1 (MQ=255)
                                                                               ttACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctg                                                                       <  2:53368/150‑1 (MQ=255)
                                                                                tACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctgc                                                                      >  2:1420444/1‑150 (MQ=255)
                                                                                tACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctgc                                                                      >  2:1066645/1‑150 (MQ=255)
                                                                                 aCTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctgcg                                                                     >  2:1159022/1‑150 (MQ=255)
                                                                                 aCTCTGAAAACCGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGCGAGCCGCGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGCgcgatgcg                                                                     >  2:30741/1‑150 (MQ=255)
                                                                                   tctGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctgcgct                                                                   >  2:777133/1‑150 (MQ=255)
                                                                                    ctGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGcgctgcgctg                                                                  >  1:447931/1‑150 (MQ=255)
                                                                                    cagAAAAAGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCAATCGCGGGAGCCGGGGGGCGTATGGGACGCCAGTTGATTCCGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGctattcgctg                                                                  >  2:1195892/3‑150 (MQ=255)
                                                                                     tGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTgg                                                                 <  2:1177721/150‑1 (MQ=255)
                                                                                      gAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGa                                                                <  1:1497948/150‑1 (MQ=255)
                                                                                        aaaCGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGc                                                              >  1:264283/1‑150 (MQ=255)
                                                                                            ggTCTATGCAAATTAACAAATGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGa                                                          <  2:1205451/150‑1 (MQ=255)
                                                                                            ggTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGAGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGa                                                          >  1:721646/1‑150 (MQ=255)
                                                                                               cTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAgg                                                       <  2:836684/150‑1 (MQ=255)
                                                                                               cTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAgg                                                       <  1:808600/150‑1 (MQ=255)
                                                                                                tATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGa                                                      >  2:242168/1‑150 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    >  2:380722/1‑150 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    <  2:1224968/150‑1 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    <  2:801869/150‑1 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    <  2:1352036/150‑1 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    >  2:836083/1‑150 (MQ=255)
                                                                                                  tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                    >  2:1413656/1‑150 (MQ=255)
                                                                                                   gCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtct                                                   >  2:122740/1‑150 (MQ=255)
                                                                                                         tAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTa                                             <  2:1445019/150‑1 (MQ=255)
                                                                                                          aaCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTAc                                            <  2:526227/150‑1 (MQ=255)
                                                                                                           aCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACt                                           <  1:1065971/150‑1 (MQ=255)
                                                                                                            cAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTg                                          >  2:243720/1‑150 (MQ=255)
                                                                                                             aaaaGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTgg                                         >  1:250415/1‑150 (MQ=255)
                                                                                                               aaGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGc                                       <  1:790428/150‑1 (MQ=255)
                                                                                                                agagAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCa                                      >  1:891943/1‑150 (MQ=255)
                                                                                                                 gagaATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCac                                     >  1:800984/1‑149 (MQ=255)
                                                                                                                  agaATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGc                                    >  1:140621/1‑150 (MQ=255)
                                                                                                                   gaATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCg                                   >  2:1458386/1‑150 (MQ=255)
                                                                                                                     aTAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGAc                                 <  1:232042/150‑1 (MQ=255)
                                                                                                                        gCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGcc                              >  1:1338033/1‑150 (MQ=255)
                                                                                                                            tGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTg                          <  2:1093697/150‑1 (MQ=255)
                                                                                                                            tGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTg                          <  2:622045/150‑1 (MQ=255)
                                                                                                                              cATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAg                        >  2:1516637/1‑150 (MQ=255)
                                                                                                                              cATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAg                        >  1:1285516/1‑150 (MQ=255)
                                                                                                                                 gatgCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTg                     >  2:810102/1‑150 (MQ=255)
                                                                                                                                 gatgCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTg                     >  2:562922/1‑150 (MQ=255)
                                                                                                                                   tgCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGc                   >  1:1099992/1‑150 (MQ=255)
                                                                                                                                    gCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGcc                  <  2:1155838/150‑1 (MQ=255)
                                                                                                                                     cAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGTAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCg                 <  2:140621/150‑1 (MQ=255)
                                                                                                                                      aaaCATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCgg                <  1:1159022/150‑1 (MQ=255)
                                                                                                                                      aaaCATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCgg                >  1:648925/1‑150 (MQ=255)
                                                                                                                                       ggctcccGCGATGGCAACGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGa               >  1:320306/6‑150 (MQ=255)
                                                                                                                                        gctcccGCGATGGCAACGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAg              >  1:1102789/5‑150 (MQ=255)
                                                                                                                                         cATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTTTGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGc             <  1:677248/150‑1 (MQ=255)
                                                                                                                                         cATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGc             <  1:836083/150‑1 (MQ=255)
                                                                                                                                         cATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGc             <  1:615003/150‑1 (MQ=255)
                                                                                                                                           cccGCGATGGCAACGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCg           >  1:779747/2‑150 (MQ=255)
                                                                                                                                             cgcgTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCggg         <  1:242779/150‑1 (MQ=255)
                                                                                                                                              gcgTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGa        >  1:894361/1‑150 (MQ=255)
                                                                                                                                               cgTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGaa       <  1:1458386/150‑1 (MQ=255)
                                                                                                                                                gTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGaaa      <  1:70326/150‑1 (MQ=255)
                                                                                                                                                    ccATCGCGGTAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAg  <  1:1195892/150‑1 (MQ=255)
                                                                                                                                                    ccATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAg  <  2:53014/150‑1 (MQ=255)
                                                                                                                                                    ccATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAg  <  1:1322978/150‑1 (MQ=255)
                                                                                                                                                    ccATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAg  <  2:250415/150‑1 (MQ=255)
                                                                                                                                                    ccATCGCGGGAGCCGGGGGGAGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACac  >  1:517920/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAG  >  NC_000913/28251‑28548

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