Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,589 C→A 7.6% A198S (GCA→TCA) ‡ 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,5890CA 7.6% 194.4 / 6.2 66A198S (GCA→TCA) ‡malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

GCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAA  >  NC_000913/1698444‑1698737
                                                                                                                                                 |                                                                                                                                                    
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   >  1:1026033/1‑150 (MQ=255)
                                                                                                                                               tGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGaa   >  2:1237726/1‑150 (MQ=255)
                                                                                                                                                gCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCCCGGGAAGCGaaa  >  1:725917/1‑150 (MQ=255)
                                                                                                                                                gCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGaaa  >  1:1467115/1‑150 (MQ=255)
                                                                                                                                                gCACGGGGTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGCACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGaaa  >  2:1211902/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAA  >  NC_000913/1698444‑1698737

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: