Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 28,396 T→A 5.3% V8D (GTT→GAT)  dapB → 4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,3960TA 5.3% 281.9 / 4.8 94V8D (GTT→GAT) dapB4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base T (47/42);  new base A (3/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

ATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAA  >  NC_000913/28247‑28544
                                                                                                                                                     |                                                                                                                                                    
aTGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGt                                                                                                                                                      >  1:836684/1‑150 (MQ=255)
    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:842602/1‑150 (MQ=255)
                  cctccctGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCgg                                                                                                                                    >  1:801869/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                                                                                                                   <  1:630462/150‑1 (MQ=255)
                                   aaTTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGgccgcc                                                                                                                   <  2:1451865/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                                                                                                                                                  <  1:863426/103‑1 (MQ=255)
                                                   cGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGcc                                                                                                                                                  >  2:863426/1‑103 (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:1066645/1‑150 (MQ=255)
                                                                                    tACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctgc                                                                  >  2:1420444/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                                                   <  1:808600/150‑1 (MQ=255)
                                                                                                   cTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAgg                                                   <  2:836684/150‑1 (MQ=255)
                                                                                                    tATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGa                                                  >  2:242168/1‑150 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                <  2:1352036/150‑1 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                >  2:836083/1‑150 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                >  2:1413656/1‑150 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                <  2:1224968/150‑1 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                <  2:801869/150‑1 (MQ=255)
                                                                                                      tGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                >  2:380722/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                    >  1:1285516/1‑150 (MQ=255)
                                                                                                                                  cATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAg                    >  2:1516637/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)
                                                                                                                                                     |                                                                                                                                                    
ATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAA  >  NC_000913/28247‑28544

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: