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

TTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAC  >  NC_000913/3339683‑3339975
                                                                                                                                                  |                                                                                                                                                   
ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGgcg                                                                                                                                                  >  1:232479/1‑149 (MQ=255)
ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGgcg                                                                                                                                                  >  3:416028/1‑149 (MQ=255)
    tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTGATCCGTTTGCGCACTGTATACAGGCGCGaa                                                                                                                                              >  1:309038/1‑149 (MQ=255)
    tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGaa                                                                                                                                              >  3:286783/1‑149 (MQ=255)
     catcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAAc                                                                                                                                             <  4:125995/149‑1 (MQ=255)
      atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACg                                                                                                                                            >  4:212454/1‑149 (MQ=255)
      atcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACg                                                                                                                                            >  4:332774/1‑149 (MQ=255)
       tcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACgg                                                                                                                                           <  1:34897/149‑2 (MQ=255)
         aCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACgaga                                                                                                                                         <  2:289526/149‑1 (MQ=255)
          ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACggat                                                                                                                                         >  2:72679/1‑145 (MQ=255)
          ccGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACggat                                                                                                                                         <  1:72679/148‑4 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                             >  1:135469/1‑149 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                             >  2:186641/1‑149 (MQ=255)
                     aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAat                                                                                                                             >  4:50186/1‑149 (MQ=255)
                           ggCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtctc                                                                                                                       <  4:416028/149‑1 (MQ=255)
                            gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCa                                                                                                                      <  4:372736/149‑1 (MQ=255)
                                          gTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATGACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAt                                                                                                        >  4:60195/1‑149 (MQ=255)
                                                      aGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATc                                                                                            <  3:332774/149‑1 (MQ=255)
                                                         ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTgg                                                                                         >  1:348281/1‑149 (MQ=255)
                                                           aTAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTg                                                                                       <  1:186641/149‑1 (MQ=255)
                                                              ggCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATGACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGa                                                                                    <  3:60195/149‑1 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   >  2:72181/1‑149 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   <  2:56557/149‑1 (MQ=255)
                                                               gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                                   >  3:33318/1‑149 (MQ=255)
                                                                cGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAtt                                                                                  <  3:50186/149‑1 (MQ=255)
                                                                 gACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         <  2:30196/109‑1 (MQ=255)
                                                                 gACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         >  1:30196/1‑109 (MQ=255)
                                                                    gTTGTCAAATACGTTCATATCATTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACg                                                                                                                                            <  3:143844/87‑1 (MQ=25)
                                                                    gTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGCCTGTATACAGTGCGCAAACg                                                                                                                                            >  4:143844/1‑87 (MQ=37)
                                                                       gTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGc                                                                           <  3:212454/149‑1 (MQ=255)
                                                                        tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACga                                                                                                                                           <  4:80836/84‑1 (MQ=255)
                                                                        tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACga                                                                                                                                           >  3:80836/1‑84 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                        >  2:402135/1‑133 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                                        <  1:402135/133‑1 (MQ=255)
                                                                          aaaTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCaa                                                                        >  1:456740/1‑149 (MQ=255)
                                                                            aTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATc                                                                      >  1:158925/1‑149 (MQ=255)
                                                                             tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                           >  3:379309/1‑111 (MQ=255)
                                                                             tACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAg                                                                                                           <  4:379309/111‑1 (MQ=255)
                                                                                 ttCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACg                                                                 >  4:300465/1‑149 (MQ=255)
                                                                                     tatCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGt                                                             >  1:384742/1‑149 (MQ=255)
                                                                                      atCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTa                                                            >  4:108497/1‑149 (MQ=255)
                                                                                       tCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAg                                                           >  4:267472/1‑149 (MQ=255)
                                                                                         aGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcg                                                         >  2:429885/1‑149 (MQ=255)
                                                                                          gTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                                        >  4:273549/1‑149 (MQ=255)
                                                                                           tGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtc                                                       >  3:17042/1‑149 (MQ=255)
                                                                                            gAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcg                                                      <  2:232479/149‑1 (MQ=255)
                                                                                             aacaacGCCCCGGACCGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcgt                                                     >  4:125793/1‑149 (MQ=255)
                                                                                                     cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                            <  2:380537/86‑1 (MQ=255)
                                                                                                     cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                                            >  1:380537/1‑86 (MQ=255)
                                                                                                      ccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGAt                                            <  3:108497/149‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          >  4:229791/1‑145 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          <  1:453857/145‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          >  2:453857/1‑145 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGc                                          <  3:229791/145‑1 (MQ=255)
                                                                                                            tGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGa                                      >  4:430013/1‑149 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         >  1:51822/1‑64 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATtc                                                                                                                         <  2:51822/64‑1 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                               >  2:281077/1‑122 (MQ=255)
                                                                                                              aaaTAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGa                                                               <  1:281077/122‑1 (MQ=255)
                                                                                                                 tAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGcc                                 >  4:84713/1‑149 (MQ=255)
                                                                                                                  aaCATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGccc                                >  2:265220/1‑149 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                >  3:61186/1‑92 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTg                                                                                <  4:61186/92‑1 (MQ=255)
                                                                                                                           aTCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGc                       >  1:307493/1‑149 (MQ=255)
                                                                                                                               gTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTTACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTg                   <  1:429885/149‑1 (MQ=255)
                                                                                                                                 ttGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGAt                 <  3:84713/149‑1 (MQ=255)
                                                                                                                                     gcACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTc             >  3:302990/1‑149 (MQ=255)
                                                                                                                                         tGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                            <  4:169630/82‑1 (MQ=255)
                                                                                                                                         tGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCg                                                                            >  3:169630/1‑82 (MQ=255)
                                                                                                                                          gTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTcg        >  2:98535/1‑149 (MQ=255)
                                                                                                                                               cAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCa   >  1:371349/1‑149 (MQ=255)
                                                                                                                                                aGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAc  >  4:185536/1‑149 (MQ=255)
                                                                                                                                                aGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAc  <  2:348281/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAGGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTGATCTTCCCTCGCGGCAC  >  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: