Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,604,638 C→T 5.2% R198C (CGC→TGC)  lsrD → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,6380CT 5.2% 248.2 / 2.3 76R198C (CGC→TGC) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (38/38)
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.11e-01

CCGATGGCGTTTACAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTG  >  NC_000913/1604491‑1604781
                                                                                                                                                   |                                                                                                                                               
ccGTTGGCGTTTACAGATTTCGCTAACCTGGATGTGATGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTGCTTTTCTGGCTCTGGCTGCATAATACGCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCcgc                                                                                                                                                 <  5:362983/148‑1 (MQ=255)
    tGGCGTTTACAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTg                                                                                                                                            <  6:279288/149‑1 (MQ=255)
         tttAGAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTGCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCGCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCt                                                                                                                                       <  4:57213/149‑1 (MQ=255)
         tttACAGAATACGCAAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATCTCTCCTCGTTTACTGGCTCTGGCTGCATAAAACCCAAGCCGGCCGTAATGTCTGTTTGATTGGGCAAAGCCAGCGCGTGgcgca                                                                                                                                       >  4:116241/1‑148 (MQ=255)
            aCAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTta                                                                                                                                    >  4:262520/1‑149 (MQ=255)
             cAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTta                                                                                                                                    <  3:352415/148‑1 (MQ=255)
              aGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTtata                                                                                                                                  <  2:19561/149‑1 (MQ=255)
              aGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTtata                                                                                                                                  <  4:347520/149‑1 (MQ=255)
               gATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAg                                                                                                                                 >  3:311226/1‑149 (MQ=255)
                 tttCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGTGCAAAGCCCGCGCGTGGCGCTTTATAgcg                                                                                                                               <  4:304285/149‑1 (MQ=255)
                   tCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAgcgc                                                                                                                              <  1:138423/148‑1 (MQ=255)
                     gCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGa                                                                                                                            <  7:315180/148‑1 (MQ=255)
                         aCCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTc                                                                                                                         >  7:313275/1‑147 (MQ=255)
                           cTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAACGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGt                                                                                                                     <  1:79037/149‑1 (MQ=255)
                           cTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGt                                                                                                                     <  7:86383/149‑1 (MQ=255)
                                    cTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGTCGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTAc                                                                                                            <  2:270495/149‑1 (MQ=255)
                                        gACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGTGCATGGCGCTTTATAGCGCGAtt                                                                                                                          <  2:197773/131‑1 (MQ=255)
                                        gACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGTGCATGGCGCTTTATAGCGCGAtt                                                                                                                          >  1:197773/1‑131 (MQ=255)
                                        gACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTa                                                                                                        >  7:389386/1‑149 (MQ=255)
                                           tCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTAtgt                                                                                                     >  2:242107/1‑149 (MQ=255)
                                             cccGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTAtgt                                                                                                     <  6:47888/147‑1 (MQ=255)
                                             cccGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTgc                                                                                                   >  6:155112/1‑149 (MQ=255)
                                             cccGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTgc                                                                                                   >  8:382373/1‑149 (MQ=255)
                                                    cgcgGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCGCTGGCTTCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGCCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTAt                                                                                            <  2:247877/149‑1 (MQ=255)
                                                     gCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATg                                                                                           <  1:82990/149‑1 (MQ=255)
                                                        gATTACCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCa                                                                                        >  4:93653/1‑149 (MQ=255)
                                                          ttATCTTCCTGATATGTCTCCTCGTTTTCTGGGTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATg                                                                                      <  3:56361/149‑1 (MQ=255)
                                                            aTCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGAc                                                                                    >  3:283084/1‑149 (MQ=255)
                                                               ttCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACggg                                                                                 >  7:138279/1‑149 (MQ=255)
                                                                tCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGACATGACgggg                                                                                >  4:141648/1‑149 (MQ=255)
                                                                  cTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCt                                                                              >  1:309612/1‑149 (MQ=255)
                                                                    gATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGTGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTgg                                                                            >  4:73966/1‑149 (MQ=255)
                                                                     atatGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGc                                                                           <  3:184887/149‑1 (MQ=255)
                                                                        tGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGc                                                                               >  4:5593/1‑142 (MQ=255)
                                                                        tGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGc                                                                               <  3:5593/142‑1 (MQ=255)
                                                                             cctcGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCgg                                                                   >  7:221779/1‑149 (MQ=255)
                                                                               tcGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCt                                                                                                          <  2:35922/108‑1 (MQ=255)
                                                                               tcGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCt                                                                                                          >  1:35922/1‑108 (MQ=255)
                                                                                 gTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTcgc                                                               >  7:210533/1‑149 (MQ=255)
                                                                                    ttCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTcgccg                                                             <  3:96759/148‑1 (MQ=255)
                                                                                      ctggctctggctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCt                                                           <  1:261810/148‑1 (MQ=255)
                                                                                         gctctggctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCtgtg                                                        <  8:313275/148‑1 (MQ=255)
                                                                                         gctctggctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCtgtg                                                        <  8:210533/148‑1 (MQ=255)
                                                                                         gctctggctGCATAAAACCAATGCCGGACGTATTGTGTTTTTGATTGGTCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCtgtg                                                        <  4:311226/148‑1 (MQ=255)
                                                                                               gctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTgg                                                 >  3:40847/1‑149 (MQ=255)
                                                                                                ctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCACTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTgg                                                 >  1:394707/1‑148 (MQ=255)
                                                                                                 tGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTa                                               >  1:117842/1‑149 (MQ=255)
                                                                                                   cATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCg                                                                    >  2:242737/1‑126 (MQ=255)
                                                                                                   cATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCg                                                                    <  1:242737/126‑1 (MQ=255)
                                                                                                     tAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCg                                            >  6:330322/1‑148 (MQ=255)
                                                                                                           ccATGCCGGACGTAATGTG‑TTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTgg                                    >  4:135654/1‑149 (MQ=255)
                                                                                                            cATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTgg                                    >  2:288978/1‑149 (MQ=255)
                                                                                                            cATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTgg                                    <  7:382373/149‑1 (MQ=255)
                                                                                                            cATGCCGGACGTAATGGGTTTTTGATTGGGCATAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGTTATCGTATTTTgg                                    <  3:116241/149‑1 (MQ=255)
                                                                                                             aTGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGt                                   <  1:242107/149‑1 (MQ=255)
                                                                                                                ccGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCa                                <  2:117842/149‑1 (MQ=255)
                                                                                                                ccGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCa                                <  3:262520/149‑1 (MQ=255)
                                                                                                                        aaTGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTc                        >  5:36199/1‑149 (MQ=255)
                                                                                                                         aTGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTcc                       >  1:254042/1‑149 (MQ=255)
                                                                                                                         aTGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTcc                       <  4:200185/149‑1 (MQ=255)
                                                                                                                             gtTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGAtc                   >  6:381362/1‑149 (MQ=255)
                                                                                                                             gtTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGAtc                   >  2:356402/1‑149 (MQ=255)
                                                                                                                              tttttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGct                                                                                                >  5:67882/1‑71 (MQ=255)
                                                                                                                              tttttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGct                                                                                                <  6:67882/71‑1 (MQ=255)
                                                                                                                              tttttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGAtct                  >  3:36727/1‑149 (MQ=255)
                                                                                                                                tttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCg                >  1:307466/1‑149 (MQ=255)
                                                                                                                                 ttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCg                <  3:93653/148‑1 (MQ=255)
                                                                                                                                   gATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGt              <  4:283084/148‑1 (MQ=255)
                                                                                                                                    aTTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGc            >  1:45215/1‑149 (MQ=255)
                                                                                                                                    aTTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGc            <  2:254042/149‑1 (MQ=255)
                                                                                                                                     ttGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcg           >  3:243663/1‑149 (MQ=255)
                                                                                                                                       gggCAAAGCCCGTGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGGCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgt          <  3:73966/148‑1 (MQ=255)
                                                                                                                                        ggCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgtcg        >  5:103203/1‑149 (MQ=255)
                                                                                                                                        ggCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCACTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgtc         <  2:394707/148‑1 (MQ=255)
                                                                                                                                         gCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgtcgt       >  3:137341/1‑149 (MQ=255)
                                                                                                                                         gCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgtcgt       >  5:18797/1‑149 (MQ=255)
                                                                                                                                              gCCCGCGCGTGGCGCGTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTg  <  6:310414/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGATGGCGTTTACAGATTTCGCTAACCTGGATGTGCTGGGACTCCCCGTTCCGCTGATTATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTG  >  NC_000913/1604491‑1604781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑