Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,339,828:1 +T 8.5% coding (238/810 nt) 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,8281.T 8.5% 282.6 / 22.4 71coding (238/810 nt)mlaFABC transporter maintaining OM lipid asymmetry, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (40/25);  new base T (2/4);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAC  >  NC_000913/3339683‑3339975
                                                                                                                                                  |                                                                                                                                                   
ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑gcg                                                                                                                                                  >  1:232479/1‑149 (MQ=255)
ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑gcg                                                                                                                                                  >  3:416028/1‑149 (MQ=255)
    tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTGATCCGTTTGCGCACTGTATACAG‑GCGCGaa                                                                                                                                              >  1:309038/1‑149 (MQ=255)
    tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGaa                                                                                                                                              >  3:286783/1‑149 (MQ=255)
     catcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAAc                                                                                                                                             <  4:125995/149‑1 (MQ=255)
      atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACg                                                                                                                                            >  4:212454/1‑149 (MQ=255)
      atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACg                                                                                                                                            >  4:332774/1‑149 (MQ=255)
       tcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACgg                                                                                                                                           <  1:34897/149‑2 (MQ=255)
         aCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACgaga                                                                                                                                         <  2:289526/149‑1 (MQ=255)
          ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACggat                                                                                                                                         >  2:72679/1‑145 (MQ=255)
          ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACggat                                                                                                                                         <  1:72679/148‑4 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAat                                                                                                                             >  1:135469/1‑149 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAat                                                                                                                             >  2:186641/1‑149 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAat                                                                                                                             >  4:50186/1‑149 (MQ=255)
                           ggCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtctc                                                                                                                       <  4:416028/149‑1 (MQ=255)
                            gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCa                                                                                                                      <  4:372736/149‑1 (MQ=255)
                                          gTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATGACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAt                                                                                                        >  4:60195/1‑149 (MQ=255)
                                                      aGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATc                                                                                            <  3:332774/149‑1 (MQ=255)
                                                         ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTgg                                                                                         >  1:348281/1‑149 (MQ=255)
                                                           aTAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTg                                                                                       <  1:186641/149‑1 (MQ=255)
                                                              ggCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATGACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGa                                                                                    <  3:60195/149‑1 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   >  2:72181/1‑149 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   <  2:56557/149‑1 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   >  3:33318/1‑149 (MQ=255)
                                                                cGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAtt                                                                                  <  3:50186/149‑1 (MQ=255)
                                                                 gACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         <  2:30196/109‑1 (MQ=255)
                                                                 gACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         >  1:30196/1‑109 (MQ=255)
                                                                    gTTGTCAAATACGTTCATATCATTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACg                                                                                                                                            <  3:143844/87‑1 (MQ=25)
                                                                    gTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACg                                                                                                                                            >  4:143844/1‑87 (MQ=37)
                                                                       gTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGc                                                                           <  3:212454/149‑1 (MQ=255)
                                                                        tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACga                                                                                                                                           <  4:80836/84‑1 (MQ=255)
                                                                        tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACga                                                                                                                                           >  3:80836/1‑84 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                        >  2:402135/1‑133 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                        <  1:402135/133‑1 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCaa                                                                        >  1:456740/1‑149 (MQ=255)
                                                                            aTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATc                                                                      >  1:158925/1‑149 (MQ=255)
                                                                             tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                           >  3:379309/1‑111 (MQ=255)
                                                                             tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                           <  4:379309/111‑1 (MQ=255)
                                                                                 ttCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACg                                                                 >  4:300465/1‑149 (MQ=255)
                                                                                     tatCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGt                                                             >  1:384742/1‑149 (MQ=255)
                                                                                      atCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTa                                                            >  4:108497/1‑149 (MQ=255)
                                                                                       tCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAg                                                           >  4:267472/1‑149 (MQ=255)
                                                                                         aGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcg                                                         >  2:429885/1‑149 (MQ=255)
                                                                                          gTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                                        >  4:273549/1‑149 (MQ=255)
                                                                                           tGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtc                                                       >  3:17042/1‑149 (MQ=255)
                                                                                            gAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcg                                                      <  2:232479/149‑1 (MQ=255)
                                                                                             aacaacGCCCCGGACCGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcgt                                                     >  4:125793/1‑149 (MQ=255)
                                                                                                     cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                            <  2:380537/86‑1 (MQ=255)
                                                                                                     cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                            >  1:380537/1‑86 (MQ=255)
                                                                                                      ccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGAt                                            <  3:108497/149‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          >  4:229791/1‑145 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          <  1:453857/145‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          >  2:453857/1‑145 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          <  3:229791/145‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGa                                      >  4:430013/1‑149 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         >  1:51822/1‑64 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         <  2:51822/64‑1 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                               >  2:281077/1‑122 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                               <  1:281077/122‑1 (MQ=255)
                                                                                                                 tAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGcc                                 >  4:84713/1‑149 (MQ=255)
                                                                                                                  aaCATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGccc                                >  2:265220/1‑149 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                >  3:61186/1‑92 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                <  4:61186/92‑1 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGc                       >  1:307493/1‑149 (MQ=255)
                                                                                                                               gTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTg                   <  1:429885/149‑1 (MQ=255)
                                                                                                                                 ttGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGAt                 <  3:84713/149‑1 (MQ=255)
                                                                                                                                     gcACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTc             >  3:302990/1‑149 (MQ=255)
                                                                                                                                         tGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                            <  4:169630/82‑1 (MQ=255)
                                                                                                                                         tGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                            >  3:169630/1‑82 (MQ=255)
                                                                                                                                          gTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTcg        >  2:98535/1‑149 (MQ=255)
                                                                                                                                               cAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCa   >  1:371349/1‑149 (MQ=255)
                                                                                                                                                aG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAc  >  4:185536/1‑149 (MQ=255)
                                                                                                                                                aG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAc  <  2:348281/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAC  >  NC_000913/3339683‑3339975

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