Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 462,737 C→A 5.4% Q275K (CAG→AAG)  ppiD → periplasmic folding chaperone, has an inactive PPIase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913462,7370CA 5.4% 293.8 / 2.9 87Q275K (CAG→AAG) ppiDperiplasmic folding chaperone, has an inactive PPIase domain
Reads supporting (aligned to +/- strand):  ref base C (35/47);  new base A (2/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGCCGGAACAATTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGATGCCA  >  NC_000913/462591‑462878
                                                                                                                                                  |                                                                                                                                             
tGACGCCGGAACAATTCCTCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAg                                                                                                                                             <  2:56/149‑1 (MQ=255)
  aCGCCGGAACAATTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGAc                                                                                                                                           <  4:290562/149‑1 (MQ=255)
         aaCAATTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGAGGATATACAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACTAAAACt                                                                                                                                    <  6:30305/149‑1 (MQ=255)
           cAATTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGa                                                                                                                                  <  2:116208/149‑1 (MQ=255)
            aaTTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGaa                                                                                                                                 <  7:37938/149‑1 (MQ=255)
            aaTTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGaa                                                                                                                                 <  7:226490/149‑1 (MQ=255)
            aaTTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGaa                                                                                                                                 <  6:371296/149‑1 (MQ=255)
            aaTTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGaa                                                                                                                                 <  5:9769/149‑1 (MQ=255)
              ttCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTg                                                                                                                                   >  7:260089/1‑145 (MQ=255)
              ttCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTg                                                                                                                                   >  7:253664/1‑145 (MQ=255)
                    gTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGc                                                                                                                         <  3:286174/149‑1 (MQ=255)
                          tACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAg                                                                                                                    <  1:57830/148‑1 (MQ=255)
                           aCATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCAGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCg                                                                                                                  >  4:361789/1‑149 (MQ=255)
                            cATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                                                                                                 <  4:391981/149‑1 (MQ=255)
                            cATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCg                                                                                                                  <  2:226684/148‑1 (MQ=255)
                             aTCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                                                                                                 >  4:80260/1‑148 (MQ=255)
                             aTCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                                                                                                 >  1:90604/1‑148 (MQ=255)
                                 aGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAGTTCACCCAGACGCAGCGTACCCGCTACAGCATCATCCATAACAAAACTGAAGATGAAGCGAAATCGGTACtt                                                                                                            >  4:225670/1‑149 (MQ=255)
                                   cTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtga                                                                                                          >  7:203263/1‑149 (MQ=255)
                                    tGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgat                                                                                                         >  3:360248/1‑149 (MQ=255)
                                    tGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgat                                                                                                         <  2:369086/149‑1 (MQ=255)
                                     gatatccGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtga                                                                                                          <  1:153525/142‑1 (MQ=255)
                                      gATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgatga                                                                                                       <  4:360248/149‑1 (MQ=255)
                                      gATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgatga                                                                                                       <  1:101021/149‑1 (MQ=255)
                                      gATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgat                                                                                                         <  6:329358/147‑1 (MQ=255)
                                      gATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgat                                                                                                         <  6:217419/147‑1 (MQ=255)
                                      gATGCCGCAACGATGCAGCAACCGGTTAGCGATGCAGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTtgatga                                                                                                       <  3:361789/149‑1 (MQ=255)
                                          ccGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCt                                                                                                    <  5:1152/148‑1 (MQ=255)
                                            gCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGaa                                                                                                 <  8:203263/149‑1 (MQ=255)
                                            gCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGaa                                                                                                 >  7:79469/1‑149 (MQ=255)
                                            gCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGaa                                                                                                 >  1:128561/1‑149 (MQ=255)
                                                 gATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                                                                                            <  8:260089/149‑1 (MQ=255)
                                                 gATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                                                                                            <  8:253664/149‑1 (MQ=255)
                                                 gATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                                                                                            <  7:114376/149‑1 (MQ=255)
                                                      agcaACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCAAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGt                                                                                       >  4:347780/1‑149 (MQ=255)
                                                      agcaACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCAAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGt                                                                                       >  4:353535/1‑149 (MQ=255)
                                                          aCCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtt                                                                                   >  3:75194/1‑149 (MQ=255)
                                                          aCCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtt                                                                                   <  6:325746/149‑1 (MQ=255)
                                                          aCCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtt                                                                                   >  8:43255/1‑149 (MQ=255)
                                                               ttAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGACGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACGGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgct                                                                              >  1:33256/1‑149 (MQ=255)
                                                                  gCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCa                                                                           >  1:212529/1‑149 (MQ=255)
                                                                        cGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGc                                                                      <  8:79469/148‑1 (MQ=255)
                                                                         ggATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTAATTGATGACCTGAATAAAGGCGGTGATTTTGATGCATTAgac                                                                     >  8:110641/1‑146 (MQ=255)
                                                                          gATATCAAGAGCTACAACGACCAGCATCAGGATCAATTCACCCAGCCGCCGCGTACACGCTACAGCATCATCCAGCACAAAACTGAAGATGAAGCGAAAGCGGTACTTGAAGAGCTGAATAAAGGCTGTGATTTTGCTGCATTAGcaa                                                                    >  5:255517/1‑146 (MQ=255)
                                                                               ccAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgc                                                                               <  2:411296/132‑1 (MQ=255)
                                                                               ccAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgc                                                                               >  1:411296/1‑132 (MQ=255)
                                                                                 agagCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGGTGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGaaaaa                                                            >  3:384653/1‑149 (MQ=255)
                                                                                 agagCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGaaaaa                                                            >  1:384820/1‑149 (MQ=255)
                                                                                 agagCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGaaaaa                                                            >  6:201920/1‑149 (MQ=255)
                                                                                      tactacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGc                                                       >  3:219449/1‑149 (MQ=255)
                                                                                        ctacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGcc                                                      >  4:65043/1‑148 (MQ=255)
                                                                                        ctacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGc                                                       <  5:201920/147‑1 (MQ=255)
                                                                                           cGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGat                                                   <  4:272519/148‑1 (MQ=255)
                                                                                                  cATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTAtctc                                           <  4:368016/149‑1 (MQ=255)
                                                                                                  cATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTAtctc                                           <  2:128561/149‑1 (MQ=255)
                                                                                                   aTCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATctct                                          >  7:152076/1‑149 (MQ=255)
                                                                                                    tCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTg                                         >  6:96605/1‑149 (MQ=255)
                                                                                                       gtaTCAATTCTCCCAGCCGCAGCGTAGCCGCTACAGCATCATGCAGACCAAAAGTGAAGATGAAGCGTAAGCGGTTCTTGATGAGCTGTATAAAGGCGGTGATTTTGCTGCATTAGCCAATGAAAAATCTGCCGATATTATCTCTGCTc                                      <  2:77197/147‑1 (MQ=255)
                                                                                                       ggATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                                                                                            >  7:414040/1‑95 (MQ=255)
                                                                                                       ggATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                                                                                            <  8:414040/95‑1 (MQ=255)
                                                                                                       ggATCAAATCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTc                                      <  6:114172/149‑1 (MQ=255)
                                                                                                        gATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCg                                     >  1:192685/1‑149 (MQ=255)
                                                                                                         aTCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTc                                      <  2:212529/147‑1 (MQ=255)
                                                                                                         aTCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCg                                     >  7:144902/1‑148 (MQ=255)
                                                                                                            aaTTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAATATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAAc                                 <  2:90604/149‑1 (MQ=255)
                                                                                                            aaTTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTaa                                  <  2:384820/148‑1 (MQ=255)
                                                                                                              ttCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAAcgg                               >  8:5158/1‑149 (MQ=255)
                                                                                                                 aCCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGAtattat                                              >  2:166332/1‑131 (MQ=255)
                                                                                                                 aCCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGAtattat                                              <  1:166332/131‑1 (MQ=255)
                                                                                                                 aCCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAAcggcgg                            >  4:153114/1‑149 (MQ=255)
                                                                                                                  cccAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACggcggc                           >  7:147262/1‑149 (MQ=255)
                                                                                                                  cccAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACggcggc                           <  8:152076/149‑1 (MQ=255)
                                                                                                                     aGCCGCAGCGTACCCGCTACAGCATCATCAAGACCAAAACTGAAGATGAATCTAAATCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAATAAAAATCTGCCGATATTATCTCTGCTCGTAACTGCGGCGat                        <  6:255517/149‑1 (MQ=255)
                                                                                                                       ccGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGatat                      <  3:65043/149‑1 (MQ=255)
                                                                                                                        cGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATg                     >  8:263345/1‑149 (MQ=255)
                                                                                                                         gCAGCGTACCCGCTACAGCATCCTCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATgg                    >  1:117570/1‑149 (MQ=255)
                                                                                                                             cGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGggttg                <  8:144902/149‑1 (MQ=255)
                                                                                                                                 cccGCTACAGCATCATCAAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTa            <  3:353535/149‑1 (MQ=255)
                                                                                                                                 cccGCTACAGCATCATCAAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTa            <  3:347780/149‑1 (MQ=255)
                                                                                                                                  ccGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTag           <  7:110641/149‑1 (MQ=255)
                                                                                                                                    gCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTagaa         <  7:219696/149‑1 (MQ=255)
                                                                                                                                    gCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTag           >  1:380044/1‑147 (MQ=255)
                                                                                                                                    gCTACAGCAGCATCCAGACCAAAACTGAAGATTAAGCGAAAGCTGTACTTTATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTagaa         <  3:225670/149‑1 (MQ=255)
                                                                                                                                       aCAGCATCATCCAGACCAAAACTGAAGCTGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCAGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGAt      >  8:282365/1‑149 (MQ=255)
                                                                                                                                        cAGCATCATCCAGACCAAAACTGAAGGTGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGAt      <  4:384653/148‑1 (MQ=255)
                                                                                                                                         aGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGATGc    <  5:96605/149‑1 (MQ=255)
                                                                                                                                           catcatCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGATGCCa  <  7:43255/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGACGCCGGAACAATTCCGCGTGAGTTACATCAAGCTGGATGCCGCAACGATGCAGCAACCGGTTAGCGATGCGGATATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTGCATTAGCCAAAGAAAAATCTGCCGATATTATCTCTGCTCGTAACGGCGGCGATATGGGTTGGTTAGAAGATGCCA  >  NC_000913/462591‑462878

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: