Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,339,930 C→T 5.1% G46S (GGT→AGT)  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,9300CT 5.1% 257.7 / 2.2 78G46S (GGT→AGT) mlaFABC transporter maintaining OM lipid asymmetry, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (2/2);  total (41/37)
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.51e-01

CCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTT  >  NC_000913/3339783‑3340076
                                                                                                                                                   |                                                                                                                                                  
ccccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGCCATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGCTCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTAcc                                                                                                                                                    >  5:229607/1‑148 (MQ=255)
  ccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                                                                                                                                                  <  6:327949/148‑1 (MQ=255)
   cGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGAt                                                                                                                                                 >  5:29608/1‑148 (MQ=255)
          aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCg                                                                                                                                            <  4:150087/146‑1 (MQ=255)
          aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCg                                                                                                                                            <  5:34557/146‑1 (MQ=255)
                 aTGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCa                                                                                                                                   <  8:137054/148‑1 (MQ=255)
                   gCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGgcccct                                                                                                                                  >  8:362411/1‑145 (MQ=255)
                      cATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGAAGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATc                                                                                                                              >  6:39286/1‑148 (MQ=255)
                         ccGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGcc                                                                                                                           >  5:325925/1‑148 (MQ=255)
                         ccGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGcc                                                                                                                           >  1:39342/1‑148 (MQ=255)
                            tttGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGt                                                                                                                         >  3:197730/1‑147 (MQ=255)
                            tttGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGt                                                                                                                         >  7:148667/1‑147 (MQ=255)
                              tGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGa                                                                                                                       <  6:330135/147‑1 (MQ=255)
                              tGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGa                                                                                                                       <  6:29608/147‑1 (MQ=255)
                                  cACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTc                                                                                                                  >  4:301322/1‑148 (MQ=255)
                                      gTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTc                                                                                                              >  2:386058/1‑148 (MQ=255)
                                      gTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTc                                                                                                              >  4:298634/1‑148 (MQ=255)
                                         taCAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTcgc                                                                                                            <  6:325925/147‑1 (MQ=255)
                                           cAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCgg                                                                                                          <  2:149396/147‑1 (MQ=255)
                                                 cgAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTc                                                                                                   <  5:101/148‑1 (MQ=255)
                                                    aCGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAgg                                                                                                >  6:189652/1‑148 (MQ=255)
                                                      gagaCATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGa                                                                                              <  4:197730/148‑1 (MQ=255)
                                                       agaCATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGaa                                                                                             <  5:189652/148‑1 (MQ=255)
                                                        gaCATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGaaa                                                                                            <  1:114473/148‑1 (MQ=255)
                                                            tCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAatat                                                                                         >  4:151475/1‑147 (MQ=255)
                                                             cGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAAtatta                                                                                       >  3:304742/1‑148 (MQ=255)
                                                             cGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGCCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGGCAGGGAAAtatta                                                                                       >  2:185193/1‑148 (MQ=255)
                                                                 ggAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGa                                                                                   <  8:211031/148‑1 (MQ=255)
                                                                      atTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATg                                                                              >  5:235208/1‑148 (MQ=255)
                                                                       ttCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGc                                                                             <  2:39342/148‑1 (MQ=255)
                                                                        tctcACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATg                                                                              >  8:179474/1‑146 (MQ=255)
                                                                        tctcACAATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATg                                                                              >  3:43349/1‑146 (MQ=255)
                                                                         ctcACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCa                                                                            <  2:141802/147‑1 (MQ=255)
                                                                            aCCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGa                                                                        >  6:306603/1‑148 (MQ=255)
                                                                            aCCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGGCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGa                                                                        <  3:294234/148‑1 (MQ=255)
                                                                             ccATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGAt                                                                       >  6:181409/1‑148 (MQ=255)
                                                                               aTCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAg                                                                           >  4:69672/1‑142 (MQ=255)
                                                                               aTCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAg                                                                           <  3:69672/142‑1 (MQ=255)
                                                                               aTCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTg                                                                     <  3:353211/148‑1 (MQ=255)
                                                                               aTCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTg                                                                     <  1:386058/148‑1 (MQ=255)
                                                                                 cGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAat                                                                                           <  4:217775/124‑1 (MQ=255)
                                                                                 cGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAat                                                                                           >  3:217775/1‑124 (MQ=255)
                                                                                   aaaaGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCAc                                                                 <  6:235208/148‑1 (MQ=255)
                                                                                    aaaGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCAcg                                                                <  3:298634/148‑1 (MQ=255)
                                                                                         aTCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTa                                                            >  8:267556/1‑147 (MQ=255)
                                                                                         aTCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGGCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTa                                                            >  7:328160/1‑147 (MQ=255)
                                                                                             cACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACt                                                       >  8:56303/1‑148 (MQ=255)
                                                                                              aCCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATc                                                                                     <  5:293250/117‑1 (MQ=255)
                                                                                              aCCATGATCTGGGGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTaaa                                                          <  2:270660/144‑1 (MQ=255)
                                                                                                     tCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATcgc                                               <  6:229607/148‑1 (MQ=255)
                                                                                                      cTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGGCAGGGCAATATTATCGAAGCTGCAGCGATTGCCACGCGTAAAACTGACATcgcg                                              >  8:264107/1‑148 (MQ=255)
                                                                                                          tGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACTGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCata                                          >  8:263231/1‑148 (MQ=255)
                                                                                                          tGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACTGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCGCGCGGCACGGTCAGGGAAATATTATCGAATATGCAGAGATTGACAAGCGTAAATCTGAAATCGCGCata                                          >  8:263232/1‑148 (MQ=255)
                                                                                                           gCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCatat                                         >  6:373994/1‑148 (MQ=255)
                                                                                                             gATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCg                                       <  1:174222/148‑1 (MQ=255)
                                                                                                             gATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGAGCGCCGTGAGCTTCCCTCGCGGCACGGTCAGTGAACTATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCg                                       <  7:264107/148‑1 (MQ=255)
                                                                                                                   cccGCCAATCAGACGGAGTAGCGTCGTTTTACTGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTaa                                  <  7:263231/147‑1 (MQ=255)
                                                                                                                   cccGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTaa                                  <  3:301322/147‑1 (MQ=255)
                                                                                                                   cccGCAAATCAGACGGAGTAGCGTCGTTTTACTGATGCCCGATTGCACCATGATCGACGTTATCTTCCATCGCGGCACGGTCATGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTaa                                  <  7:263232/147‑1 (MQ=255)
                                                                                                                    ccGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAAt                                >  1:202584/1‑148 (MQ=255)
                                                                                                                      gCCAATCAGACGGAGTAGCGTCGTTTTACCGACGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTc                              >  5:95627/1‑148 (MQ=255)
                                                                                                                       ccAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCg                             <  8:148667/148‑1 (MQ=255)
                                                                                                                           tCAGACGGAGTAGCGTCGTTTTACCCATGCCCGATGGCCACATTATCGAAGTGATATTACCTCGCGTCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCca                          >  2:255260/1‑147 (MQ=255)
                                                                                                                             aGACGGAGTAGCGTCGTTTTACCGATGCCCGATTGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAg                       <  7:267556/148‑1 (MQ=255)
                                                                                                                             aGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAg                       <  7:56303/148‑1 (MQ=255)
                                                                                                                               aCGGAGTAGCGTCGTTTTAACGATGCCCGATGACCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTAAGGTAAATATTATCGAAGATGCAGCGAGTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGAc                     >  6:156714/1‑148 (MQ=255)
                                                                                                                                cGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATTATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAACACTGCGATCGCGCATATCGACTAAATTCGCCACAGACt                    >  8:40752/1‑148 (MQ=255)
                                                                                                                                 ggAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTg                   <  5:39286/148‑1 (MQ=255)
                                                                                                                                    gTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTc                >  5:69309/1‑148 (MQ=255)
                                                                                                                                        cgtcgtTTTACCGATGCCCGATGGCCCCATGTTCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCAt             >  6:51342/1‑147 (MQ=255)
                                                                                                                                        cgtcgtTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCAt             >  4:160946/1‑147 (MQ=255)
                                                                                                                                         gtcgtTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATaa           <  6:69309/148‑1 (MQ=255)
                                                                                                                                         gtcgtTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATa            >  6:292569/1‑147 (MQ=255)
                                                                                                                                           cgtTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAAtt         >  7:141352/1‑148 (MQ=255)
                                                                                                                                            gtTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTc        <  4:352962/148‑1 (MQ=255)
                                                                                                                                             ttttACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCa       >  1:116741/1‑148 (MQ=255)
                                                                                                                                             ttttACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCa       >  2:12781/1‑148 (MQ=255)
                                                                                                                                              tttACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATcgcg                                              <  5:196683/108‑1 (MQ=255)
                                                                                                                                              tttACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATcgcg                                              >  6:196683/1‑108 (MQ=255)
                                                                                                                                                  ccGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCtt  <  3:151475/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTT  >  NC_000913/3339783‑3340076

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