Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,333,058 A→G 6.9% L21P (CTG→CCG) ‡ 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,0580AG 6.9% 335.3 / 11.1 101L21P (CTG→CCG) ‡rpmA50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base A (44/50);  new base G (3/4);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 106 total aligned reads displayed.

TTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAG  >  NC_000913/3332910‑3333206
                                                                                                                                                    |                                                                                                                                                    
ttCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                      <  4:220942/149‑1 (MQ=255)
ttCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                      >  3:437569/1‑149 (MQ=255)
ttCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                      >  2:879945/1‑149 (MQ=255)
ttCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCa                                                                                                                                                      >  2:525505/1‑149 (MQ=255)
     aTTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACTGATTCGCCACCTAAACGCTTAACGCCCAGGCGt                                                                                                                                                 <  2:1216884/149‑1 (MQ=255)
       ttCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGttt                                                                                                                                               <  1:1073545/149‑1 (MQ=255)
       ttCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGttt                                                                                                                                               <  1:810212/149‑1 (MQ=255)
          aCTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTAAACGCCTGGGCGTTAAGc                                                                                                                                            <  2:1060535/149‑1 (MQ=255)
           cTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCt                                                                                                                                           <  1:654803/149‑1 (MQ=255)
            tttACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAGTTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCtt                                                                                                                                          >  4:850384/1‑149 (MQ=255)
             ttACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTc                                                                                                                                         <  1:879945/149‑1 (MQ=255)
              tACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCt                                                                                                                                        <  1:525505/149‑1 (MQ=255)
                       cTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                               <  3:850384/149‑1 (MQ=255)
                        tttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAc                                                                                                                              >  1:1086732/1‑149 (MQ=255)
                         ttAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAcc                                                                                                                             >  2:430211/1‑149 (MQ=255)
                            gCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                               >  2:750277/1‑144 (MQ=255)
                            gCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                               <  1:750277/144‑1 (MQ=255)
                            gCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAc                                                                                                                              >  4:386633/1‑145 (MQ=255)
                            gCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGAc                                                                                                                              <  3:386633/145‑1 (MQ=255)
                              aaaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTAAACGCCTGGGCGTTAAGc                                                                                                                                            >  4:828240/1‑129 (MQ=255)
                              aaaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTAAACGCCTGGGCGTTAAGc                                                                                                                                            <  3:828240/129‑1 (MQ=255)
                               aaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATcgc                                                                                                                                 <  1:500140/139‑1 (MQ=255)
                               aaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATcgc                                                                                                                                 >  2:500140/1‑139 (MQ=255)
                               aaCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACg                                                                                                                       >  4:943261/1‑149 (MQ=255)
                                aCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGtt                                                                                                                          <  1:184900/145‑1 (MQ=255)
                                aCAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGtt                                                                                                                          >  2:184900/1‑145 (MQ=255)
                                    agTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTgg                                                                                                                  >  1:556311/1‑149 (MQ=255)
                                     gtgtGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGa                                                                                                                 >  2:544940/1‑149 (MQ=255)
                                      tgtgGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGc                                                                                                                                            >  3:255781/1‑121 (MQ=255)
                                      tgtgGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGc                                                                                                                                            <  4:255781/121‑1 (MQ=255)
                                          gTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAgccgc                                                                                                            >  1:838638/1‑149 (MQ=255)
                                          gTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAgccgc                                                                                                            >  2:190574/1‑149 (MQ=255)
                                             aCGACCGCAACCTACGTTAGAGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAg                                                                                                         >  2:1374527/1‑149 (MQ=255)
                                               gACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTc                                                                                                                                         <  2:1314378/115‑1 (MQ=255)
                                               gACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTc                                                                                                                                         >  1:1314378/1‑115 (MQ=255)
                                               gACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGcc                                                                                                       >  2:1292376/1‑149 (MQ=255)
                                                   gCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtttt                                                                                                   <  1:430211/149‑1 (MQ=255)
                                                    cAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCttttt                                                                                                  <  1:544940/149‑1 (MQ=255)
                                                      aCCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAt                                                                                                <  1:190574/149‑1 (MQ=255)
                                                      aCCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAt                                                                                                <  1:1292376/149‑1 (MQ=255)
                                                        cTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtgt                                                                                              <  2:556311/149‑1 (MQ=255)
                                                          aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGc                                                                                            <  2:1086732/149‑1 (MQ=255)
                                                          aCGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGc                                                                                            <  4:437569/149‑1 (MQ=255)
                                                                gcgcCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTg                                                                                      >  3:278929/1‑149 (MQ=255)
                                                                  gcCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATc                                                                                                                                   <  4:1320570/102‑1 (MQ=255)
                                                                  gcCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATc                                                                                                                                   >  3:1320570/1‑102 (MQ=255)
                                                                  gcCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtg                                                                                               <  3:1311343/138‑1 (MQ=255)
                                                                  gcCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtg                                                                                               >  4:1311343/1‑138 (MQ=255)
                                                                        gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTAc                                                                                                                        <  1:1166461/107‑1 (MQ=255)
                                                                        gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTAc                                                                                                                        >  2:1166461/1‑107 (MQ=255)
                                                                        gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtgt                                                                                              >  2:238198/1‑133 (MQ=255)
                                                                        gTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAtgt                                                                                              <  1:238198/133‑1 (MQ=255)
                                                                         tGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctctc                                                                             >  2:874845/1‑149 (MQ=255)
                                                                         tGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtctctc                                                                             <  1:874845/149‑1 (MQ=255)
                                                                            aaTTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcctc                                                                          <  4:278929/149‑1 (MQ=255)
                                                                             aTTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAt                                                                                                <  3:1460378/126‑1 (MQ=255)
                                                                             aTTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTAt                                                                                                >  4:1460378/1‑126 (MQ=255)
                                                                                    aCCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtttt                                                                                                   >  4:649063/1‑116 (MQ=255)
                                                                                    aCCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCtttt                                                                                                   <  3:649063/116‑1 (MQ=255)
                                                                                     ccACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAg                                                                 <  3:943261/149‑1 (MQ=255)
                                                                                         gTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                               <  1:1455552/83‑1 (MQ=255)
                                                                                         gTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGa                                                                                                                               >  2:1455552/1‑83 (MQ=255)
                                                                                          ttGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCTGGGCGTTAAGc                                                                                                                                            <  2:865510/69‑1 (MQ=255)
                                                                                          ttGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCTGGGCGTTAAGc                                                                                                                                            >  1:865510/1‑69 (MQ=39)
                                                                                          ttGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCt                                                            <  2:838638/149‑1 (MQ=255)
                                                                                            gACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGa                                                          >  3:1058066/1‑149 (MQ=255)
                                                                                             acgaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGAt                                                         <  1:1042912/149‑1 (MQ=255)
                                                                                               gaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGa                                                                                                                                      <  2:646411/70‑1 (MQ=255)
                                                                                               gaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGa                                                                                                                                      >  1:646411/1‑70 (MQ=255)
                                                                                               gaacgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGc                                                       >  1:605903/1‑149 (MQ=255)
                                                                                                 acgaTGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCa                                                     <  1:1374527/149‑1 (MQ=255)
                                                                                                   gatgatGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTAc                                                                                                                        <  2:891555/80‑1 (MQ=255)
                                                                                                   gatgatGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTAc                                                                                                                        >  1:891555/1‑80 (MQ=255)
                                                                                                     tgatgCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCTGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTaa                                                 <  2:605903/149‑1 (MQ=255)
                                                                                                                gCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGttt                                                                                                                                               <  1:459410/44‑1 (MQ=255)
                                                                                                                gCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGttt                                                                                                                                               >  2:459410/1‑44 (MQ=255)
                                                                                                                gCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCa                                      >  2:703171/1‑149 (MQ=255)
                                                                                                                 ccAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCttt                                                                                                    >  3:535074/1‑86 (MQ=255)
                                                                                                                 ccAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCttt                                                                                                    <  4:535074/86‑1 (MQ=255)
                                                                                                                     aaCGTATTCGCCACCGAAACGCTTAACGCCCAGGCGt                                                                                                                                                 <  4:1253537/37‑1 (MQ=255)
                                                                                                                     aaCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGt                                                                                                                                                 >  3:1253537/1‑37 (MQ=255)
                                                                                                                     aaCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCAtt                                                                                        >  3:441534/1‑94 (MQ=255)
                                                                                                                     aaCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCAtt                                                                                        <  4:441534/94‑1 (MQ=255)
                                                                                                                     aaCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGaa                                 >  1:148560/1‑149 (MQ=255)
                                                                                                                      aCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcc                                                                            >  3:1469909/1‑105 (MQ=255)
                                                                                                                      aCGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTctcc                                                                            <  4:1469909/105‑1 (MQ=255)
                                                                                                                              gCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTa                                                                  <  3:1002259/107‑1 (MQ=255)
                                                                                                                              gCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTa                                                                  >  4:1002259/1‑107 (MQ=255)
                                                                                                                                  ccGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAtttt                                             >  1:1392684/1‑124 (MQ=255)
                                                                                                                                  ccGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAtttt                                             <  2:1392684/124‑1 (MQ=255)
                                                                                                                                  ccGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATgg                    <  1:703171/149‑1 (MQ=255)
                                                                                                                                      aaCGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGc                                                       <  2:1311578/110‑1 (MQ=255)
                                                                                                                                      aaCGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGc                                                       >  1:1311578/1‑110 (MQ=255)
                                                                                                                                         gCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCctgct             >  2:82356/1‑149 (MQ=255)
                                                                                                                                         gCTAAACGCCTGGGCGTTAAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAt                                                <  1:906593/114‑1 (MQ=255)
                                                                                                                                         gCTAAACGCCTGGGCGTTAAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAAt                                                >  2:906593/1‑114 (MQ=255)
                                                                                                                                             aaCGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtct                                                                                >  1:888970/1‑78 (MQ=255)
                                                                                                                                             aaCGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAAtct                                                                                <  2:888970/78‑1 (MQ=255)
                                                                                                                                               cGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTAc       <  4:1058066/149‑1 (MQ=255)
                                                                                                                                                    aGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAg  >  2:100805/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACGTTAGCGCCAGCGTGGAATTTGGTACCACGTTGACGAACGATGATGCTACCCGCCAGAACGGATTCGCCACCGAAACGCTTAACGCCCAGGCGTTTAGCTTCTGAATCGCGACCGTTACGTGTGGAGCCGCCAGCCTTTTTATGTGCCATTTGAAATCTCTCCTCAGGTCTTAGGCGCTGATGCCAGTAATTTTCACATCAGTGAACCACTGACGATGGCCCTGCTGCTTACGATAG  >  NC_000913/3332910‑3333206

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