Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,629,495 A→T 5.3% L15Q (CTG→CAG)  rbbA ← ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,4950AT 5.3% 308.5 / 5.1 94L15Q (CTG→CAG) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base A (38/51);  new base T (3/2);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGG  >  NC_000913/3629346‑3629643
                                                                                                                                                     |                                                                                                                                                    
cAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCa                                                                                                                                                      >  1:1349322/1‑150 (MQ=255)
 aaTGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAg                                                                                                                                                     <  2:1474740/150‑1 (MQ=255)
  aTGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGt                                                                                                                                                    <  1:1239854/150‑1 (MQ=255)
    gACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTGACGGCCGCCAGTTg                                                                                                                                                  <  2:236513/150‑1 (MQ=255)
    gACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTg                                                                                                                                                  <  2:482054/150‑1 (MQ=255)
     aCGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgc                                                                                                                                                 >  2:770137/1‑150 (MQ=255)
      cgcgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcg                                                                                                                                                <  2:1208785/150‑1 (MQ=255)
          ggCACCGGAAATCAACTACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                            <  2:1064396/150‑1 (MQ=255)
          ggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                            <  2:198236/150‑1 (MQ=255)
          ggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                            >  1:458159/1‑150 (MQ=255)
           gCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACa                                                                                                                                           <  2:122973/150‑1 (MQ=255)
              ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcgc                                                                                                                                               <  2:1320822/143‑1 (MQ=255)
              ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcgc                                                                                                                                               >  1:1320822/1‑143 (MQ=255)
              ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                        <  2:174586/150‑1 (MQ=255)
              ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                        >  1:600928/1‑150 (MQ=255)
               cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAACCCGACCATACGGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                                                                       >  2:684838/1‑150 (MQ=255)
                    aTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacg                                                                                                                                  >  2:889919/1‑150 (MQ=255)
                     tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCTGTCGCGCAACTGGcgggcgtg                                                                                                                                 <  1:986909/150‑5 (MQ=255)
                     tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                 <  1:691399/150‑1 (MQ=255)
                      cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                >  2:531352/1‑150 (MQ=255)
                       aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                                               >  1:904585/1‑150 (MQ=255)
                         cGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                                                             >  2:209843/1‑150 (MQ=255)
                              acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                        >  2:911407/1‑150 (MQ=255)
                              acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                        >  1:1128867/1‑150 (MQ=255)
                              acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                        >  1:1241650/1‑150 (MQ=255)
                               caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                                                       <  1:89373/150‑1 (MQ=255)
                               caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                                                       <  2:931651/150‑1 (MQ=255)
                                 aGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                     <  1:531352/150‑1 (MQ=255)
                                   cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                                                                   >  1:421956/1‑150 (MQ=255)
                                     cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATg                                                                                                                 >  1:1026098/1‑150 (MQ=255)
                                       aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                                               <  2:1387557/150‑1 (MQ=255)
                                       aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                                               <  2:84705/150‑1 (MQ=255)
                                        cTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                              >  2:794745/1‑150 (MQ=255)
                                        cTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                              <  2:904585/150‑1 (MQ=255)
                                          tCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                                                            <  1:1016181/150‑1 (MQ=255)
                                              gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGc                                                                                                        <  1:884788/150‑1 (MQ=255)
                                              gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGc                                                                                                        <  1:943867/150‑1 (MQ=255)
                                                cGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCgg                                                                                                      <  1:976239/150‑1 (MQ=255)
                                                    gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                                                                  >  2:449151/1‑150 (MQ=255)
                                                    gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                                                                  >  2:590323/1‑150 (MQ=255)
                                                     tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                                                                 <  2:944425/150‑1 (MQ=255)
                                                           ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                           >  1:529355/1‑150 (MQ=255)
                                                             aaTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                     <  2:1198199/122‑1 (MQ=255)
                                                             aaTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                     >  1:1198199/1‑122 (MQ=255)
                                                                cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                      <  1:911407/150‑1 (MQ=255)
                                                                cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                      <  1:1241034/150‑1 (MQ=255)
                                                                 aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                                                                     <  1:495502/150‑1 (MQ=255)
                                                                  gCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                    >  2:1358164/1‑150 (MQ=255)
                                                                    ccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                                                  >  1:575200/1‑150 (MQ=255)
                                                                         cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                                                                             <  1:209843/150‑1 (MQ=255)
                                                                         cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                                                                             >  1:987400/1‑150 (MQ=255)
                                                                              aGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACgg                                                                        <  2:271176/150‑1 (MQ=255)
                                                                                    ccGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                  >  1:1137207/1‑150 (MQ=255)
                                                                                      ggAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                >  2:1278170/1‑150 (MQ=255)
                                                                                       gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                               <  1:934159/150‑1 (MQ=255)
                                                                                       gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                               <  2:1349322/150‑1 (MQ=255)
                                                                                         atatCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcat                                                             <  2:1008866/150‑1 (MQ=255)
                                                                                           atCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc                                                           <  2:1137207/150‑1 (MQ=255)
                                                                                            tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattca                                                          >  1:134202/1‑150 (MQ=255)
                                                                                             cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                                                                          >  2:1000317/1‑133 (MQ=255)
                                                                                             cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                                                                          <  1:1000317/133‑1 (MQ=255)
                                                                                                gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                      <  2:1241650/150‑1 (MQ=255)
                                                                                                gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                      <  1:1027221/150‑1 (MQ=255)
                                                                                                gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                      >  2:498387/1‑150 (MQ=255)
                                                                                                  gTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCa                                                    <  2:1533453/150‑1 (MQ=255)
                                                                                                    gATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc                                                  >  2:39317/1‑150 (MQ=255)
                                                                                                     atatTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc                                                 >  1:72272/1‑150 (MQ=255)
                                                                                                          gTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCt                                            >  2:328907/1‑150 (MQ=255)
                                                                                                           ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                           <  1:1358164/150‑1 (MQ=255)
                                                                                                           ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                           <  1:449151/150‑1 (MQ=255)
                                                                                                           ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                           <  1:859559/150‑1 (MQ=255)
                                                                                                            tCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTAc                                          >  1:1094836/1‑150 (MQ=255)
                                                                                                             cAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCCGGCCACGGAAGTTCTTCATTCACCCGCACGCCCGCGAcc                                         >  2:985931/1‑149 (MQ=255)
                                                                                                                   aaCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGc                                   <  2:575200/150‑1 (MQ=255)
                                                                                                                    aCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCa                                  >  2:1395587/1‑150 (MQ=255)
                                                                                                                       gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAAc                               <  2:239942/150‑1 (MQ=255)
                                                                                                                       gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAAc                               <  1:1397591/150‑1 (MQ=255)
                                                                                                                           ttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGt                           >  1:702849/1‑150 (MQ=255)
                                                                                                                              cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                        <  1:770137/150‑1 (MQ=255)
                                                                                                                                tAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGa                      <  2:600928/150‑1 (MQ=255)
                                                                                                                                   tGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACat                   >  1:766899/1‑150 (MQ=255)
                                                                                                                                       ggCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTc               <  1:684838/150‑1 (MQ=255)
                                                                                                                                        gCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTcc              <  2:529355/150‑1 (MQ=255)
                                                                                                                                        gCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTcc              <  1:889919/150‑1 (MQ=255)
                                                                                                                                         cTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa             <  1:842931/150‑1 (MQ=255)
                                                                                                                                            aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt          <  2:1433484/150‑1 (MQ=255)
                                                                                                                                             gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCCCGAGGTCGTCAGGCCACGGAAGTTCTTCATTCCCCCGCACCCACGCTACGCCCGGCAACCCGGTTTTGACATATTCCAGAtg         >  1:828511/3‑150 (MQ=255)
                                                                                                                                             gtcccgccTGTCGCGCAACTGGCGGG‑CGTGAGCCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgc        <  1:52089/148‑1 (MQ=38)
                                                                                                                                              tcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGAGTCATTGCGGCAACCTCACCGCGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCAGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgc        >  2:1063208/2‑150 (MQ=255)
                                                                                                                                               cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct       >  2:1187018/1‑150 (MQ=255)
                                                                                                                                               cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct       >  1:488857/1‑150 (MQ=255)
                                                                                                                                               cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct       >  1:1267777/1‑150 (MQ=255)
                                                                                                                                               cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct       >  2:1313809/1‑150 (MQ=255)
                                                                                                                                                   ccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGctgctg   <  2:1128867/150‑1 (MQ=255)
                                                                                                                                                    cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTgg  <  2:72272/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGG  >  NC_000913/3629346‑3629643

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