Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,333,057 C→T 5.7% L21L (CTG→CTA) ‡ rpmA ← 50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,0570CT 5.7% 417.2 / 9.6 123L21L (CTG→CTA) ‡rpmA50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base C (60/56);  new base T (4/3);  total (64/59)
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.41e-01

Only 100 of 126 total aligned reads displayed.

ACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATA  >  NC_000913/3332920‑3333205
                                                                                                                                         |                                                                                                                                                    
aCTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGc                                                                                                                                           >  4:1861923/1‑149 (MQ=255)
aCTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTAAACGCCTGGGCGTTAAGc                                                                                                                                           <  3:1095205/149‑1 (MQ=255)
        gTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGa                                                                                                                                     <  4:1086727/147‑1 (MQ=255)
          cTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATcg                                                                                                                                 >  4:775169/1‑149 (MQ=255)
             cTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              <  3:1021165/149‑1 (MQ=255)
              tttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              <  3:1467357/148‑1 (MQ=255)
              tttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              >  4:1467357/1‑148 (MQ=255)
              tttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAc                                                                                                                             <  3:1582570/149‑1 (MQ=255)
               ttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAcc                                                                                                                            <  3:1861923/149‑1 (MQ=255)
                 aGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTa                                                                                                                                             <  2:500809/130‑1 (MQ=255)
                 aGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTa                                                                                                                                             >  1:500809/1‑130 (MQ=255)
                  gCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGtt                                                                                                                         >  1:210512/1‑149 (MQ=255)
                    aaaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTAc                                                                                                                       <  3:1006925/149‑1 (MQ=255)
                     aaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACg                                                                                                                      >  2:1487386/1‑149 (MQ=255)
                      aCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     >  1:1569942/1‑149 (MQ=255)
                         gagTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCAAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGtgtg                                                                                                                  <  2:1301289/149‑1 (MQ=255)
                          agTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTc                                                                                                                                        >  3:978579/1‑126 (MQ=255)
                          agTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              >  4:1571505/1‑136 (MQ=255)
                           gtgtGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGa                                                                                                                <  2:210512/149‑1 (MQ=255)
                               ggTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                     <  1:1468073/108‑1 (MQ=255)
                               ggTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                     >  2:1468073/1‑108 (MQ=255)
                                   aCGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAg                                                                                                        >  1:1137357/1‑149 (MQ=255)
                                     gACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCTGGGCGttaag                                                                                                                                            >  3:1393462/1‑108 (MQ=255)
                                     gACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCTGGGCGTTAAg                                                                                                                                            <  4:1393462/111‑1 (MQ=255)
                                       ccGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtt                                                                                                    >  1:837181/1‑149 (MQ=255)
                                          cAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCa                                                                                                         <  1:1621838/141‑1 (MQ=255)
                                          cAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCa                                                                                                         >  2:1621838/1‑141 (MQ=255)
                                               tACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTAAACGCCTGGGCGTTAAGc                                                                                                                                           >  1:1665151/1‑102 (MQ=255)
                                                aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGc                                                                                           <  4:666728/149‑1 (MQ=255)
                                                aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGc                                                                                           <  3:1357846/149‑1 (MQ=255)
                                                aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACAAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              >  4:11558/1‑114 (MQ=255)
                                                aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACAAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                              <  3:11558/114‑1 (MQ=255)
                                                   ttAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCAt                                                                                        >  2:1833295/1‑149 (MQ=255)
                                                         ccAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtgt                                                                                             <  4:1610215/138‑1 (MQ=255)
                                                         ccAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtgt                                                                                             >  3:1610215/1‑138 (MQ=255)
                                                         ccAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGaaa                                                                                  <  1:1487386/149‑1 (MQ=255)
                                                              gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctct                                                                             >  3:220265/1‑149 (MQ=255)
                                                              gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctct                                                                             <  4:1390717/149‑1 (MQ=255)
                                                              gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctct                                                                             >  4:1729603/1‑149 (MQ=255)
                                                               tGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctctc                                                                            <  2:1569942/149‑1 (MQ=255)
                                                                ggAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcc                                                                           <  2:837181/149‑1 (MQ=255)
                                                                ggAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcc                                                                           <  4:549422/149‑1 (MQ=255)
                                                                ggAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcc                                                                           <  3:1729603/149‑1 (MQ=255)
                                                                 gAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcct                                                                          >  4:284504/1‑149 (MQ=255)
                                                                 gAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcct                                                                          >  3:1617064/1‑149 (MQ=255)
                                                                        gTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATcg                                                                                                                                 <  1:1581189/87‑1 (MQ=255)
                                                                        gTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATcg                                                                                                                                 >  2:1581189/1‑87 (MQ=255)
                                                                         tACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCtt                                                                  >  4:1706290/1‑149 (MQ=255)
                                                                          aCCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTa                                                                 <  3:775169/149‑1 (MQ=255)
                                                                            cACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAgg                                                               >  3:1550421/1‑149 (MQ=255)
                                                                              cGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtt                                                                                                    >  4:364423/1‑110 (MQ=255)
                                                                              cGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtt                                                                                                    <  3:364423/110‑1 (MQ=255)
                                                                                   acgaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTa                                                                 <  1:1220534/140‑1 (MQ=255)
                                                                                   acgaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGAt                                                        <  2:434437/149‑1 (MQ=255)
                                                                                       acgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCa                                                    <  1:1833295/149‑1 (MQ=255)
                                                                                         gatgatGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAg                                                                >  1:1753004/1‑135 (MQ=255)
                                                                                         gatgatGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAg                                                                <  2:1753004/135‑1 (MQ=255)
                                                                                         gatgatGCTACCCGCAAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGt                                                  >  4:1402647/1‑149 (MQ=255)
                                                                                          atgatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     >  3:876855/1‑81 (MQ=255)
                                                                                          atgatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     >  3:529837/1‑81 (MQ=255)
                                                                                          atgatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     <  4:529837/81‑1 (MQ=255)
                                                                                          atgatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     <  4:876855/81‑1 (MQ=255)
                                                                                            gatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCt                                                                                                                                       >  3:1715291/1‑61 (MQ=255)
                                                                                            gatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCt                                                                                                                                       <  4:1715291/61‑1 (MQ=255)
                                                                                             atgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAtt                                              <  4:1617064/149‑1 (MQ=255)
                                                                                             atgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGATTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAtt                                              <  2:1137357/149‑1 (MQ=255)
                                                                                              tgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCttt                                                                                                   <  2:476648/95‑1 (MQ=255)
                                                                                              tgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCttt                                                                                                   >  1:476648/1‑95 (MQ=255)
                                                                                              tgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAttt                                             <  3:284504/149‑1 (MQ=255)
                                                                                                 tACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACgt                                                                                                                     <  3:480819/74‑1 (MQ=255)
                                                                                                 tACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcct                                                                          >  4:265030/1‑117 (MQ=255)
                                                                                                 tACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTca                                          >  3:1176746/1‑149 (MQ=255)
                                                                                                 tACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTca                                          <  3:1663152/149‑1 (MQ=255)
                                                                                                   cccGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTcaca                                        >  1:575200/1‑149 (MQ=255)
                                                                                                     cGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTaa                                                <  3:1067718/139‑1 (MQ=255)
                                                                                                      gCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAt                                               >  2:1036408/1‑139 (MQ=255)
                                                                                                      gCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAt                                               <  1:1036408/139‑1 (MQ=255)
                                                                                                         aGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGAGTTTAGCTTCTGAATCGCGACCGTTATGTGTGCAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTg                                  <  4:1176746/149‑1 (MQ=255)
                                                                                                          gAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTaa                                                <  4:968927/134‑1 (MQ=255)
                                                                                                          gAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTaa                                                >  3:968927/1‑134 (MQ=255)
                                                                                                          gAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGa                                 >  2:660588/1‑149 (MQ=255)
                                                                                                          gAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGa                                 >  1:155253/1‑149 (MQ=255)
                                                                                                           aaCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGaa                                >  2:523412/1‑149 (MQ=255)
                                                                                                            aCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAAc                               >  1:66086/1‑149 (MQ=255)
                                                                                                               gATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCAc                            <  4:1550421/149‑1 (MQ=255)
                                                                                                                   cGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGAc                        >  2:925823/1‑149 (MQ=255)
                                                                                                                       aCCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAg                                                                       >  2:492079/1‑98 (MQ=255)
                                                                                                                       aCCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAg                                                                       <  1:492079/98‑1 (MQ=255)
                                                                                                                       aCCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATg                    >  1:1466550/1‑149 (MQ=255)
                                                                                                                       aCCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATg                    >  4:417104/1‑149 (MQ=255)
                                                                                                                        ccGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATgg                   <  3:1402647/149‑1 (MQ=255)
                                                                                                                          gAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGt                                                                     <  1:1205312/97‑1 (MQ=255)
                                                                                                                          gAAACGCTTAACGCCCAGGCGTTTAGATTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGt                                                                     >  2:1205312/1‑97 (MQ=255)
                                                                                                                               gCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCctgct            <  1:443913/149‑1 (MQ=255)
                                                                                                                                 taaaCGCCTGGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCctgct            >  4:88721/3‑147 (MQ=255)
                                                                                                                                   aaCGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCtgctgc          <  1:1404936/147‑1 (MQ=255)
                                                                                                                                   aaCGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCtgctgc          >  2:1404936/1‑147 (MQ=255)
                                                                                                                                    aCGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTa       <  2:1635180/149‑1 (MQ=255)
                                                                                                                                       cccAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGa    <  2:66086/149‑1 (MQ=255)
                                                                                                                                       cccAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGa    <  1:925823/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
ACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATA  >  NC_000913/3332920‑3333205

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