Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,038 T→C 5.4% E51E (GAA→GAG rlpA ← septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,0380TC 5.4% 241.9 / 3.6 74E51E (GAA→GAGrlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base C (2/2);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-02

TGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCAC  >  NC_000913/664891‑665179
                                                                                                                                                   |                                                                                                                                             
tGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCACCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGtt                                                                                                                                              >  1:724189/1‑149 (MQ=255)
 gCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCTTCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTc                                                                                                                                             <  2:152717/149‑1 (MQ=255)
   ttCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGa                                                                                                                                           >  2:749372/1‑149 (MQ=255)
   ttCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGa                                                                                                                                           >  4:181638/1‑149 (MQ=255)
        cAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAAcgc                                                                                                                                      <  3:528956/149‑1 (MQ=255)
        cAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGCATACTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAAcgc                                                                                                                                      <  4:518318/149‑1 (MQ=255)
          gagGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCgg                                                                                                                                    >  4:158654/1‑149 (MQ=255)
          gagGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCgg                                                                                                                                    >  4:308642/1‑149 (MQ=255)
                gTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGc                                                                                                                              >  1:170444/1‑149 (MQ=255)
                     gTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTGTGTAGCTTTTACCGTCGCGCTGGTAATACTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGc                                                                                                                         >  1:112045/1‑149 (MQ=255)
                        aCTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTaa                                                                                                                      <  1:307577/149‑1 (MQ=255)
                              gggTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGCGACGGATCTTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCaa                                                                                                                <  3:689940/149‑1 (MQ=255)
                              gggTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCaa                                                                                                                <  1:280767/149‑1 (MQ=255)
                              gggTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCaa                                                                                                                <  3:17787/149‑1 (MQ=255)
                              gggTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCCCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCaa                                                                                                                <  4:462898/149‑1 (MQ=255)
                                   cGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACtata                                                                                                           <  2:390509/149‑1 (MQ=255)
                                   cGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACtata                                                                                                           <  2:562243/149‑1 (MQ=255)
                                    ggCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAg                                                                                                          >  3:413161/1‑149 (MQ=255)
                                    ggCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAg                                                                                                          >  4:395900/1‑149 (MQ=255)
                                        tcatAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGcc                                                                                                      >  1:121963/1‑149 (MQ=255)
                                           tAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGtt                                                                                                   >  1:724655/1‑149 (MQ=255)
                                           tAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGtt                                                                                                   <  4:244853/149‑1 (MQ=255)
                                            aGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGc                                                                                                       >  1:169963/1‑144 (MQ=255)
                                            aGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGc                                                                                                       <  2:169963/144‑1 (MQ=255)
                                                     ccAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCa                                                                                         <  2:88051/149‑1 (MQ=255)
                                                             gCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGCTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCgg                                                                                 >  4:556779/1‑149 (MQ=255)
                                                                tGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTat                                                                              <  2:121963/149‑2 (MQ=255)
                                                                  gCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                            >  1:114200/1‑149 (MQ=255)
                                                                  gCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                                            >  1:282649/1‑149 (MQ=255)
                                                                   cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                                                           >  3:386845/1‑149 (MQ=255)
                                                                    tAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACtt                                                                          >  2:205099/1‑149 (MQ=255)
                                                                      aaTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                                        <  3:113182/149‑1 (MQ=255)
                                                                       aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAcc                                                                       <  2:24467/149‑1 (MQ=255)
                                                                       aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAcc                                                                       <  3:181638/149‑1 (MQ=255)
                                                                          gagaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTc                                                                    >  1:599972/1‑149 (MQ=255)
                                                                           agaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                                   <  4:760877/149‑1 (MQ=255)
                                                                             aCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGt                                                                 <  2:170444/149‑1 (MQ=255)
                                                                             aCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGt                                                                 <  3:395900/149‑1 (MQ=255)
                                                                              cGGATCCTGCACGATTTTGTAGCTTTTATAGACGCGCTCGTAATCCTGATTTGCCGTCGCGTGCAGTTGTTCGATACGCGGGTCCACCCCGTTAATTTCAACTATATGTTCGTTACATACCGCAGGCTGCGGAACACTTACCGTCTgat                                                                <  2:112045/149‑3 (MQ=255)
                                                                               ggATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTg                                                               >  3:628660/1‑149 (MQ=255)
                                                                                  tCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                                                            <  4:628660/149‑1 (MQ=255)
                                                                                       cACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                       >  2:704042/1‑149 (MQ=255)
                                                                                       cACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                                                       >  3:653185/1‑149 (MQ=255)
                                                                                          gATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                    >  3:597015/1‑149 (MQ=255)
                                                                                             tttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc                                                 <  4:82729/149‑1 (MQ=255)
                                                                                               tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCtt                                               <  3:158352/149‑1 (MQ=255)
                                                                                                gTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg                                              <  1:455652/149‑1 (MQ=255)
                                                                                                   gCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAc                                           <  1:749372/149‑1 (MQ=255)
                                                                                                   gCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAc                                           <  4:209410/149‑1 (MQ=255)
                                                                                                     ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGtt                                                                                                   <  4:760620/91‑1 (MQ=255)
                                                                                                     ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGtt                                                                                                   >  3:760620/1‑91 (MQ=255)
                                                                                                     ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                                         >  3:605723/1‑149 (MQ=255)
                                                                                                     ttttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                                         >  1:591950/1‑149 (MQ=255)
                                                                                                          ccGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc                                    <  1:712511/149‑1 (MQ=255)
                                                                                                            gTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCTGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGa                                  <  1:525598/149‑1 (MQ=255)
                                                                                                            gTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGa                                  <  3:158654/149‑1 (MQ=255)
                                                                                                               gcgcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCa                               <  4:597015/149‑1 (MQ=255)
                                                                                                                     gaaaTCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGtt                                                                >  1:376094/3‑110 (MQ=255)
                                                                                                                     gTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGtt                                                                <  2:376094/110‑1 (MQ=255)
                                                                                                                     gTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTg                         >  1:201328/1‑149 (MQ=255)
                                                                                                                         tCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAgg                                                                                                         >  2:681596/1‑65 (MQ=255)
                                                                                                                         tCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAgg                                                                                                         <  1:681596/65‑1 (MQ=255)
                                                                                                                           cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGa                   >  1:500890/1‑149 (MQ=255)
                                                                                                                           cTGATTTGCCGTCGCGTTCAGTGGCTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                         <  3:703802/95‑1 (MQ=255)
                                                                                                                           cTGATTTGCCGTCGCGTTCAGTGGCTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                         >  4:703802/1‑95 (MQ=255)
                                                                                                                            tGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatgt                                                                            >  1:344136/1‑88 (MQ=255)
                                                                                                                            tGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatgt                                                                            <  2:344136/91‑4 (MQ=255)
                                                                                                                                 tGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGa             <  3:264580/149‑1 (MQ=255)
                                                                                                                                  gCCGTCGCGTTCAGTGGCTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGAt            <  3:556779/149‑1 (MQ=255)
                                                                                                                                       cgcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGa       <  1:205099/149‑1 (MQ=255)
                                                                                                                                            tCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCAc  <  2:599972/149‑1 (MQ=255)
                                                                                                                                               gTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc          >  4:497979/1‑138 (MQ=255)
                                                                                                                                               gTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc          <  3:497979/138‑1 (MQ=255)
                                                                                                                                                  gTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCCGACCATCATCGCTTGTACATGCCGCGa                                  <  4:13607/111‑1 (MQ=255)
                                                                                                                                                  gTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCCGACCATCATCGCTTGTACATGCCGCGa                                  >  3:13607/1‑111 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCAC  >  NC_000913/664891‑665179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: