Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,587 T→A 8.2% A198A (GCA→GCT) ‡ malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,5870TA 8.2% 169.2 / 5.1 61A198A (GCA→GCT) ‡malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/3);  total (29/32)
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.46e-01

GCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAA  >  NC_000913/1698444‑1698736
                                                                                                                                               |                                                                                                                                                     
gCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACgg                                                                                                                                                 <  1:1686677/150‑1 (MQ=255)
   cTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGtt                                                                                                                                              <  2:728748/150‑1 (MQ=255)
     gATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgg                                                                                                                                            <  2:705914/150‑3 (MQ=255)
      aTGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                                           <  2:1618249/150‑1 (MQ=255)
      aTGGTCGGGGTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                                           <  2:1479032/150‑1 (MQ=255)
       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATg                                                                                                                                          >  2:449121/1‑150 (MQ=255)
       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttgggg                                                                                                                                          <  1:1108303/150‑5 (MQ=255)
        ggTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttggggg                                                                                                                                         <  1:1363604/150‑6 (MQ=255)
            gggTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGaa                                                                                                                                     <  1:364695/150‑1 (MQ=255)
                  tGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACAGCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                               <  2:71704/150‑1 (MQ=255)
                  tGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                               <  2:1716057/150‑1 (MQ=255)
                  tGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                               >  2:1175678/1‑150 (MQ=255)
                    aCGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGAACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGccc                                                                                                                             >  1:926065/1‑150 (MQ=255)
                           aGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                                                                                                      <  2:35091/150‑1 (MQ=255)
                              gcCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTGGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                                   <  1:1009737/150‑1 (MQ=255)
                                   gATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGTACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGAt                                                                                                              >  1:447658/1‑150 (MQ=255)
                                    aTAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCCGCCAGGCGATc                                                                                                             >  1:1433517/1‑150 (MQ=255)
                                       gCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                          >  1:427678/1‑150 (MQ=255)
                                          tCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCATCCCAGCCAGGCGATCCGCTGa                                                                                                       >  2:1307295/1‑150 (MQ=255)
                                              cgGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAAATTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACTGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGcc                                                                                                   <  1:1299796/150‑1 (MQ=255)
                                               ggCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGccc                                                                                                  <  2:1186649/150‑1 (MQ=255)
                                                gcttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCa                                                                                                 >  2:761038/1‑150 (MQ=255)
                                                            ctggAACTGCACTCCAACCCCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGAGCCGCTGAGGCCCAGTGCGAATgaga                                                                                     >  2:699364/1‑150 (MQ=255)
                                                                  gTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTcc                                                                               <  2:1145850/150‑1 (MQ=255)
                                                                     cACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                            >  1:1354293/1‑150 (MQ=255)
                                                                      aCTCCAACACCCAATCGCTATGAAACGGCCGGCCTAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTca                                                                           <  2:1078916/150‑1 (MQ=255)
                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTca                                                                           <  2:447658/150‑1 (MQ=255)
                                                                        tCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                         <  2:127622/150‑1 (MQ=255)
                                                                           aaCACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCaacaac                                                                      <  1:497448/150‑1 (MQ=255)
                                                                            acacCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACt                                                                     <  1:416984/150‑1 (MQ=255)
                                                                             cacCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtg                                                                    >  1:1628107/1‑150 (MQ=255)
                                                                               tccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCATCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                  <  1:1530382/149‑1 (MQ=255)
                                                                               cccAATCGCTATGAAACGGCATGCCAAATTTTAGTATAGTTGCACAATAGCCCCCAACCCGGTCTGCACGGGTTAATGATGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCTAATGAGATGCTCCGTCAACAACtgtg                                                                  <  1:748836/150‑1 (MQ=255)
                                                                                    tCGCTCTGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGcc                                                             <  2:1413287/150‑1 (MQ=255)
                                                                                         aTGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCGATTGCGAAGGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                        >  1:124947/1‑150 (MQ=255)
                                                                                          tGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtg                                                       >  1:1670764/1‑150 (MQ=255)
                                                                                           gAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgt                                                      >  1:902600/1‑150 (MQ=255)
                                                                                              aCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgttgt                                                   <  1:1175678/150‑1 (MQ=255)
                                                                                                    ggCCAAATTTTCGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGc                                             <  1:699364/150‑1 (MQ=255)
                                                                                                     gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCg                                            >  1:1164055/1‑150 (MQ=255)
                                                                                                        aaaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAc                                         >  2:528117/1‑150 (MQ=255)
                                                                                                         aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAcc                                        <  1:1618658/150‑1 (MQ=255)
                                                                                                         aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACAGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAcc                                        <  2:1588845/150‑1 (MQ=255)
                                                                                                          aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCg                                       >  1:556300/1‑150 (MQ=255)
                                                                                                                ttagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAACGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCaa                                 <  1:711121/149‑1 (MQ=255)
                                                                                                                  agagTTGCACAATAGCCCCCAACCCGTTCTGCCCGGGTTAATGATGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAca                               <  1:449121/150‑1 (MQ=255)
                                                                                                                  agagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAca                               >  2:799638/1‑150 (MQ=255)
                                                                                                                    agTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAcatc                             <  2:435527/150‑1 (MQ=255)
                                                                                                                             aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTCTCAACATCATCAAGATa                    >  1:1066479/1‑150 (MQ=255)
                                                                                                                             aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATa                    >  2:1713402/1‑150 (MQ=255)
                                                                                                                             aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATa                    >  2:195633/1‑150 (MQ=255)
                                                                                                                                gCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt                 >  1:112526/1‑150 (MQ=255)
                                                                                                                                ctcattaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt                 >  1:534265/7‑150 (MQ=255)
                                                                                                                                 tcattaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTg                >  1:551892/6‑150 (MQ=255)
                                                                                                                                  ccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGCTGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGGCCGGGCGAACCGTATCAACATCATCAAGATaagggg               >  2:1188445/1‑146 (MQ=255)
                                                                                                                                  ccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGGGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTgg               >  2:1478030/1‑150 (MQ=255)
                                                                                                                                  ccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTgg               <  1:519349/150‑1 (MQ=255)
                                                                                                                                        cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTCTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACggg         >  2:473520/1‑150 (MQ=255)
                                                                                                                                         ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGa        >  2:1586607/1‑150 (MQ=255)
                                                                                                                                           gTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAg      <  2:124947/150‑1 (MQ=255)
                                                                                                                                           gTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAg      >  2:576178/1‑150 (MQ=255)
                                                                                                                                               tGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGaa  >  2:1237726/1‑150 (MQ=255)
                                                                                                                                               tGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGaa  >  1:1026033/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAA  >  NC_000913/1698444‑1698736

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