Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,629,506 A→T 9.4% P11P (CCT→CCA rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,5060AT 9.4% 337.9 / 20.1 105P11P (CCT→CCArbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base A (53/42);  new base T (6/4);  total (59/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 105 total aligned reads displayed.

ACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAATTC  >  NC_000913/3629359‑3629651
                                                                                                                                                    |                                                                                                                                                 
aCCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAg                                                                                                                                                  >  1:643782/1‑149 (MQ=255)
  cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                                >  4:508833/1‑149 (MQ=255)
  cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                                >  4:784829/1‑149 (MQ=255)
    gAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgg                                                                                                                                              <  1:557755/149‑1 (MQ=255)
    aaaaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCTCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgg                                                                                                                                              >  2:540425/2‑149 (MQ=255)
      aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCTGTCGCGCAACTGGcgggc                                                                                                                                            <  2:810088/149‑2 (MQ=255)
      aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCTGTCGCGCAACTGGcgggc                                                                                                                                            <  2:292749/149‑2 (MQ=255)
        tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacg                                                                                                                                          >  3:480551/1‑149 (MQ=255)
          aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                        >  1:297081/1‑149 (MQ=255)
           aCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                                                       <  2:511348/149‑1 (MQ=255)
           aCGACAACAAGCTCGACTTCCCGACACCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTCGCGCAACTGGcgggcgtga                                                                                                                                        <  4:3776/149‑6 (MQ=37)
             gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                                                                      <  1:845262/148‑1 (MQ=255)
             gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                                                                      >  2:845262/1‑148 (MQ=255)
             gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                                                                     >  1:400596/1‑149 (MQ=255)
                aacaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAg                                                                                                                                  >  3:489617/1‑149 (MQ=255)
                 acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGt                                                                                                                                 >  3:393027/1‑149 (MQ=255)
                 acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGt                                                                                                                                 >  1:886300/1‑149 (MQ=255)
                    aGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGAGCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                                                                              >  4:103891/1‑149 (MQ=255)
                     gCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                             >  4:106071/1‑149 (MQ=255)
                      cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTACATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAg                                                                                                                            >  1:97489/1‑149 (MQ=255)
                       tCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                                                                           >  3:451295/1‑149 (MQ=255)
                        cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCTGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                                                          <  2:809499/149‑1 (MQ=255)
                          aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGc                                                                                                                        <  4:668824/149‑1 (MQ=255)
                          aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGc                                                                                                                        <  4:480551/149‑1 (MQ=255)
                            ttCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                      <  2:252279/149‑1 (MQ=255)
                            ttCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                      <  4:385413/149‑1 (MQ=255)
                              cccGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                                                                    >  2:380265/1‑149 (MQ=255)
                               ccGACGCCGTCCGGGCCAATCAGCCCGACAATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCAt                                                                                                                   >  4:346146/1‑149 (MQ=255)
                                 gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTg                                                                                                                 >  4:127056/1‑149 (MQ=255)
                                    gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                         >  3:660292/1‑122 (MQ=255)
                                    gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                         <  4:660292/122‑1 (MQ=255)
                                       gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                                                                           >  2:508849/1‑149 (MQ=255)
                                       gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                                                                           <  1:384242/149‑1 (MQ=255)
                                              ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                                                          >  3:393082/1‑127 (MQ=255)
                                              ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                                                          <  4:393082/127‑1 (MQ=255)
                                              ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacc                                                                                                    >  4:541952/1‑149 (MQ=255)
                                               cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                                   >  1:697866/1‑149 (MQ=255)
                                               cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                                   <  4:802605/149‑1 (MQ=255)
                                                 aTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTGGCGCGACAGGCGGCACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACg                                                                                                 >  1:721131/1‑149 (MQ=255)
                                                    aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                              <  3:508833/149‑1 (MQ=255)
                                                     gCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                                                                             >  4:407533/1‑149 (MQ=255)
                                                      cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                            <  2:524881/149‑1 (MQ=255)
                                                      cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                            >  3:575629/1‑149 (MQ=255)
                                                      cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                            <  2:676127/149‑1 (MQ=255)
                                                       ccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcg                                                                                           >  3:925642/1‑149 (MQ=255)
                                                        cGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                                                          <  3:905433/149‑1 (MQ=255)
                                                         gACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtc                                                                                         <  1:521496/149‑1 (MQ=255)
                                                         gACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGAGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtc                                                                                         >  3:762384/1‑149 (MQ=255)
                                                                    gggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                              <  3:117398/149‑1 (MQ=255)
                                                                      gCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGt                                                                            >  2:767175/1‑149 (MQ=255)
                                                                          gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                        <  1:707289/149‑1 (MQ=255)
                                                                          gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                        <  4:575629/149‑1 (MQ=255)
                                                                           aaTATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                                       <  3:333095/149‑1 (MQ=255)
                                                                               tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc                                                                   >  3:657972/1‑149 (MQ=255)
                                                                                 gagaGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattcaa                                                                 >  3:717612/1‑148 (MQ=255)
                                                                                 gagaGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAc                                                                 <  3:794591/149‑1 (MQ=255)
                                                                                   gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAccc                                                               >  2:66643/1‑149 (MQ=255)
                                                                                          atTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCa                                                        >  3:97234/1‑149 (MQ=255)
                                                                                                    gcAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGccc                                              <  3:574292/149‑1 (MQ=255)
                                                                                                      aaCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCgg                                            >  2:355546/1‑149 (MQ=255)
                                                                                                       aCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGc                                           >  2:42291/1‑149 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                          >  1:327372/1‑108 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                          <  2:327372/108‑1 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 <  4:393027/149‑1 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 <  1:380265/149‑1 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 >  1:354072/1‑149 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 <  2:886300/149‑1 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 >  1:171961/1‑149 (MQ=255)
                                                                                                                 cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                                 <  2:237613/149‑1 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCTGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCCTTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCTGCCAAACCGGtttt                                >  2:759018/1‑149 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                              >  4:826582/1‑103 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                              <  3:826582/103‑1 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                <  2:293019/149‑1 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                >  4:479788/1‑149 (MQ=255)
                                                                                                                  aTAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                >  2:557575/1‑149 (MQ=255)
                                                                                                                    aaTGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGa                              <  2:513030/149‑1 (MQ=255)
                                                                                                                      tGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACa                            >  1:29776/1‑149 (MQ=255)
                                                                                                                            ttcccgtcccgccTGTCGCGCAACTGGCTGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTc                       >  2:831589/8‑149 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                     >  3:633603/1‑149 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                     >  2:190894/1‑149 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                     <  1:723854/149‑1 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                     <  2:643782/149‑1 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                     >  4:123404/1‑149 (MQ=255)
                                                                                                                             tCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAGGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGACACGGAAGGTCTTCATTAACCCGCACCCACGCTACGCCCGGCAAACCAGTTTTGACATATTCCa                     >  4:752738/1‑149 (MQ=255)
                                                                                                                               aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                <  4:285043/136‑1 (MQ=255)
                                                                                                                               aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                                >  3:285043/1‑136 (MQ=255)
                                                                                                                                 gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt                  >  2:629231/3‑149 (MQ=255)
                                                                                                                                 gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt                  >  2:388003/3‑149 (MQ=255)
                                                                                                                                  tcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtg                 >  3:11869/2‑149 (MQ=255)
                                                                                                                                   cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtt                                  <  3:558789/131‑1 (MQ=255)
                                                                                                                                   cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtt                                  >  4:558789/1‑131 (MQ=255)
                                                                                                                                    ccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct               >  2:799214/1‑149 (MQ=255)
                                                                                                                                          gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa         <  1:190894/149‑1 (MQ=255)
                                                                                                                                          gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa         <  4:451295/149‑1 (MQ=255)
                                                                                                                                          gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa         <  2:697866/149‑1 (MQ=255)
                                                                                                                                            tGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGt       <  4:395324/149‑1 (MQ=255)
                                                                                                                                                cgACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAAtt   >  1:744794/1‑149 (MQ=255)
                                                                                                                                                 gACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTTACATATTCCAGATGCTGCTGGAGTAATTc  >  4:121030/1‑149 (MQ=255)
                                                                                                                                                  aCAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                              <  2:473607/72‑1 (MQ=255)
                                                                                                                                                  aCAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                              >  1:473607/1‑72 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGTAATTC  >  NC_000913/3629359‑3629651

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