Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,339,660 C→A 5.0% G136* (GGA→TGA)  mlaF ← ABC transporter maintaining OM lipid asymmetry, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,339,6600CA 5.0% 268.9 / 2.3 79G136* (GGA→TGA) mlaFABC transporter maintaining OM lipid asymmetry, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base A (2/2);  total (40/39)
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.74e-01

TTCACCAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATC  >  NC_000913/3339512‑3339808
                                                                                                                                                    |                                                                                                                                                    
ttCACCAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTcc                                                                                                                                                      >  1:41189/1‑149 (MQ=255)
     cAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGCTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGTCATTAGTTTAGCCGCTCCACg                                                                                                                                                   >  5:264142/1‑147 (MQ=255)
     cAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACg                                                                                                                                                   >  1:147661/1‑147 (MQ=255)
     cAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACg                                                                                                                                                   >  2:383676/1‑147 (MQ=255)
     cAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCa                                                                                                                                                 <  1:383676/149‑1 (MQ=255)
       gTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGc                                                                                                                                               <  4:275473/149‑1 (MQ=255)
       gTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGc                                                                                                                                               <  7:175880/149‑1 (MQ=255)
           gCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCa                                                                                                                                           >  5:302587/1‑149 (MQ=255)
             ccATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACg                                                                                                                                         >  5:285704/1‑149 (MQ=255)
             ccAGGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACAAGATGAGATACGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCAGCCCACCGGAAAGTTCAGAAGGCAGGAGTTTAGCCGCTCCACGCAGCCCCACg                                                                                                                                         >  4:193635/1‑149 (MQ=255)
                  gTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTc                                                                                                                                    >  5:272041/1‑149 (MQ=255)
                        ggATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCtt                                                                                                                              <  7:317122/149‑1 (MQ=255)
                               gCCCAACAAAAGGTTCATCAAACATTATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTtcatc                                                                                                                          <  4:96546/146‑1 (MQ=255)
                                ccccACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCACCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTtcatcatca                                                                                                                      <  6:51576/149‑1 (MQ=255)
                                cccAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTtcatcatc                                                                                                                       <  6:285704/148‑1 (MQ=255)
                                        aaGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTa                                                                                                              >  3:233443/1‑149 (MQ=255)
                                           gTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATg                                                                                                            >  7:144263/1‑148 (MQ=255)
                                            ttCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCCCGCCGCCACCCGGCCTCCAGCTTCATCATCACCGTACTATGc                                                                                                           >  5:79101/1‑148 (MQ=255)
                                              catcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGATTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTCCTATGcaa                                                                                                         >  8:30804/1‑148 (MQ=255)
                                               atcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCAACGCAGCCCCACg                                                                                                                                         <  5:173654/115‑1 (MQ=255)
                                               atcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCAACGCAGCCCCACg                                                                                                                                         <  5:167920/115‑1 (MQ=255)
                                               atcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCAACGCAGCCCCACg                                                                                                                                         >  6:167920/1‑115 (MQ=255)
                                               atcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCAACGCAGCCCCACg                                                                                                                                         >  6:173654/1‑115 (MQ=255)
                                                tcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGcaacaa                                                                                                      <  4:156635/149‑1 (MQ=255)
                                                tcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGcaacaa                                                                                                      >  5:332065/1‑149 (MQ=255)
                                                tcaAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGcaacaa                                                                                                      >  8:343092/1‑149 (MQ=255)
                                                tcaAACATGATGAGATCCGGCTCCAGCGCAAACGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGcaacaa                                                                                                      <  7:182464/149‑1 (MQ=255)
                                                  aaaCATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATg                                                                                                    <  2:27245/149‑1 (MQ=255)
                                                   aaCATGATGAGATCCTGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGCCATTAGATTAGCCGCTCCACCCAGCCCCACCGACTCCAGCTTCATCATCACCGTACTATGCAACAAAgg                                                                                                   >  4:100523/1‑149 (MQ=255)
                                                   aaCATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCTATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATgg                                                                                                   >  4:152925/1‑149 (MQ=255)
                                                   aaCATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTCTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCACCAATgg                                                                                                   >  1:58966/1‑149 (MQ=255)
                                                       tgatgaGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCAATAGTTTAGCCGCTCCACGCAGCCCCACGCCCTCCAGCTTCATCATCACCGTACTATGCAACAATcgcg                                                                                                 >  7:197398/1‑147 (MQ=255)
                                                         atgaGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCgg                                                                                              <  2:393234/148‑1 (MQ=255)
                                                         atgaGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCgg                                                                                              >  7:317404/1‑148 (MQ=255)
                                                            agaTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGaa                                                                                           <  3:125042/148‑1 (MQ=255)
                                                             gaTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGt                                                                                         >  2:78510/1‑149 (MQ=255)
                                                              aTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGtt                                                                                        >  3:135309/1‑149 (MQ=255)
                                                              aTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGtt                                                                                        <  4:233033/149‑1 (MQ=255)
                                                              aTCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGACGAGCCATCCCACCGGAAAGTTCAGATGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATGATCACCGTACTATGCAACAATGGCGCGGGAAGtt                                                                                        <  3:100523/149‑1 (MQ=255)
                                                                 cGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTgg                                                                                      <  6:175881/148‑1 (MQ=255)
                                                                 cGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTgg                                                                                      <  6:272041/148‑1 (MQ=255)
                                                                  ggCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGt                                                                                    <  1:148046/149‑1 (MQ=255)
                                                                  ggCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCAGTAGTTTAGCCGCGCCACGCCGCCCCACGGCCTCCAGCTTCATCATCACCGTACTAGGCCACAATGGCGCGGGAAGTTggg                                                                                     >  5:56531/1‑148 (MQ=255)
                                                                   gCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCTCGCAGCCCGACGGCCTCCAGCTTCATCATCACCGTAGTATGCAACAATGGCGCGGGAAGTTGGGTa                                                                                   <  1:71251/149‑1 (MQ=255)
                                                                     tCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCCTCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCGTCCCGCTTCATCATCACCGTACTATGCAACATTGGCGCGGGAAGTGGGGTATg                                                                                 >  2:210396/1‑149 (MQ=255)
                                                                       cAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGtt                                                                               >  6:148603/1‑149 (MQ=255)
                                                                        aGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTc                                                                              <  3:273765/149‑1 (MQ=255)
                                                                         gcgcAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTT‑GGTATGTTcgc                                                                            <  1:78510/149‑1 (MQ=255)
                                                                              aTCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAg                                                                        >  3:316504/1‑149 (MQ=255)
                                                                               tCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGt                                                                       <  1:112181/149‑1 (MQ=255)
                                                                                 gCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTg                                                                      >  3:10386/1‑148 (MQ=255)
                                                                                  cACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGa                                                                    >  2:311786/1‑149 (MQ=255)
                                                                                       gcCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGc                                                               >  8:385282/1‑149 (MQ=255)
                                                                                             gcTGCACGCGGCGCCATCCCACCGGAAAGTGCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCGACGGGCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGAAGtt                                                         <  3:193635/149‑1 (MQ=255)
                                                                                             gcTGCACGCCGCGCTATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGtt                                                         <  3:152925/149‑1 (MQ=255)
                                                                                             gcTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCTTTAGTTTAGGCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCTCCGTACTATGCAACATTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGtt                                                         <  8:344189/149‑1 (MQ=255)
                                                                                             gcTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGtt                                                         <  8:292247/149‑1 (MQ=255)
                                                                                                 cACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGt                                                                       >  5:160567/1‑131 (MQ=255)
                                                                                                 cACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGt                                                                       <  6:160567/131‑1 (MQ=255)
                                                                                                 cACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTc                                                      <  2:147661/148‑1 (MQ=255)
                                                                                                 cACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTc                                                      <  4:10386/148‑1 (MQ=255)
                                                                                                   cgccgcGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCaaa                                                   <  4:233443/149‑1 (MQ=255)
                                                                                                        cgcCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGt                                              <  6:222062/149‑1 (MQ=255)
                                                                                                         gcCATCCCACGGGAAAGTTCAGCAGGCATTAGTTTAGCCGCTCCACGCAGCGCGACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGATGGGGATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGtt                                             <  2:335904/149‑1 (MQ=255)
                                                                                                          ccATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATTGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTc                                            <  6:264142/149‑1 (MQ=255)
                                                                                                                 aCCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAg                                     >  5:98265/1‑149 (MQ=255)
                                                                                                                     gAACGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGAAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGAc                                                            <  4:285146/122‑1 (MQ=255)
                                                                                                                     gAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATATGTTCATATCAGTGaa                                 >  7:281126/1‑149 (MQ=255)
                                                                                                                          ttCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACg                            >  4:42731/1‑149 (MQ=255)
                                                                                                                             aGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGccc                         >  3:111686/1‑149 (MQ=255)
                                                                                                                                  gCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGAc                    <  7:30804/149‑1 (MQ=255)
                                                                                                                                  gCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGCAc                    >  2:12850/1‑149 (MQ=255)
                                                                                                                                  gCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGATCATATCAGTGAACAACGCCCCGGAc                    >  8:169219/1‑149 (MQ=255)
                                                                                                                                    aTTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGaacaac                             <  1:326363/138‑1 (MQ=255)
                                                                                                                                    aTTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGaacaac                             >  2:326363/1‑138 (MQ=255)
                                                                                                                                    aTTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACt                   >  2:6932/1‑148 (MQ=255)
                                                                                                                                    aTTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACt                   >  8:188091/1‑148 (MQ=255)
                                                                                                                                       aggttAGCCGCTCCCCGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGGATGGTCGCGCAGTGGATAGGCGACGTTGTCAACTACGTTCATATCAGTGAACAACGCCCCGGACTGaaa               <  6:56531/146‑1 (MQ=255)
                                                                                                                                             gCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACa          <  8:317404/148‑1 (MQ=255)
                                                                                                                                              ccGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATg        >  1:262626/1‑149 (MQ=255)
                                                                                                                                                    cACGCAGCCCCACGGCCTCCCGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATc  >  7:367324/1‑149 (MQ=255)
                                                                                                                                                    cACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCAt   <  4:135309/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCACCAGTACGCCCATGGTGATGGGATCTTGCCCAACAAAAGGTTCATCAAACATGATGAGATCCGGCTCCAGCGCAATCGCACGCGCCAGCGCTGCACGCCGCGCCATCCCACCGGAAAGTTCAGAAGGCATTAGTTTAGCCGCTCCACGCAGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATC  >  NC_000913/3339512‑3339808

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: ‑