Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,368,691 G→A 5.2% A205T (GCA→ACA)  pspA → regulatory protein for phage‑shock‑protein operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,368,6910GA 5.2% 249.4 / 2.4 77A205T (GCA→ACA) pspAregulatory protein for phage‑shock‑protein operon
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base A (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTCCGTTAAC  >  NC_000913/1368544‑1368832
                                                                                                                                                   |                                                                                                                                             
cTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGatgatg                                                                                                                                               >  1:17586/1‑148 (MQ=255)
    aCAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAAt                                                                                                                                           >  7:57109/1‑148 (MQ=255)
          gCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGa                                                                                                                                     <  6:72998/148‑1 (MQ=255)
             aaCTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGaacaa                                                                                                                                 >  5:24804/1‑149 (MQ=255)
              aCTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGaacaac                                                                                                                                <  8:250421/149‑1 (MQ=255)
                  gATGACGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTg                                                                                                                              <  5:113853/147‑1 (MQ=255)
                  gATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTgg                                                                                                                             >  3:136049/1‑148 (MQ=255)
                   aTGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGACCAACTGGca                                                                                                                           >  2:317823/1‑149 (MQ=255)
                    tGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGcac                                                                                                                          <  8:197926/149‑1 (MQ=255)
                       aGCAATGGCGCGTTGCGAATCTTTCGAACGTCGTATTGACCAGATGGCAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGAGCAGTTGGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACaa                                                                                                                        <  6:254129/148‑1 (MQ=255)
                            tGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTaaaa                                                                                                                  <  8:57109/149‑1 (MQ=255)
                                  gTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaa                                                                                                            <  3:202381/149‑1 (MQ=255)
                                  gTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaa                                                                                                             >  8:197852/1‑148 (MQ=255)
                                   tttCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACTATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaaa                                                                                                           <  5:269508/149‑1 (MQ=255)
                                   tttCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaaa                                                                                                           <  6:124915/149‑1 (MQ=255)
                                      cGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATACAATCAGCGAACAACTGGCACAATTAAAAg                                                                                                                 >  1:317903/1‑140 (MQ=255)
                                      cGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATACAATCAGCGAACAACTGGCACAATTAAAAg                                                                                                                 >  5:302917/1‑140 (MQ=255)
                                      cGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATACAATCAGCGAACAACTGGCACAATTAAAAg                                                                                                                 <  6:302917/140‑1 (MQ=255)
                                      cGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATACAATCAGCGAACAACTGGCACAATTAAAAg                                                                                                                 <  2:317903/140‑1 (MQ=255)
                                       gAATCTTTCGAACGTCGTATTGACCAGATGGAATCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATg                                                                                                         >  3:174094/1‑147 (MQ=255)
                                       gAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATg                                                                                                         >  5:106344/1‑147 (MQ=255)
                                       gAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGaa                                                                                                       <  6:24804/149‑1 (MQ=255)
                                         aTCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTTGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGc                                                                                                     <  5:295671/149‑1 (MQ=255)
                                         aTCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGc                                                                                                     <  8:305891/149‑1 (MQ=255)
                                              tCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGa                                                                                                 <  4:136049/148‑1 (MQ=255)
                                                  aCGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAcaatc                                                                                            <  3:160242/149‑1 (MQ=255)
                                                   cgtcgtATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACGGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAcaatc                                                                                            >  3:209999/1‑148 (MQ=255)
                                                    gtcgtATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAcaatcaa                                                                                          <  6:106344/149‑1 (MQ=255)
                                                     tcgtATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAcaatcaat                                                                                         <  1:317823/149‑1 (MQ=255)
                                                       gtATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCaataa                                                                                       >  3:231694/1‑149 (MQ=255)
                                                       gtATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCaataa                                                                                       <  7:197852/149‑1 (MQ=255)
                                                               cAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATc                                                                                >  8:169581/1‑148 (MQ=255)
                                                                aGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCa                                                                              <  8:216217/149‑1 (MQ=255)
                                                                 gATTGAAGCGGAAGAAGAAAGACACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCAGAACAGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAg                                                                             <  4:209999/149‑1 (MQ=255)
                                                                      aagcggaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAggcggc                                                                        <  2:17586/149‑1 (MQ=255)
                                                                        gcggaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGc                                                                                                                                       <  6:188352/84‑1 (MQ=255)
                                                                        gcggaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGc                                                                                                                                       >  5:188352/1‑84 (MQ=255)
                                                                         cggaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTc                                                                     >  2:182384/1‑149 (MQ=255)
                                                                           gaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGc                                                                                                                                       >  2:76267/1‑81 (MQ=255)
                                                                           gaagcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGc                                                                                                                                       <  1:76267/81‑1 (MQ=255)
                                                                             agcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGaa                                                                 >  2:20417/1‑149 (MQ=255)
                                                                              gcAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAAc                                                                >  3:87014/1‑149 (MQ=255)
                                                                                aGAAAGCCACAGATTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAAcgc                                                              <  1:123107/149‑1 (MQ=255)
                                                                                 gAATGCCACAGGTTCGTTAAACAAAAATCGCTGGACGTTCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAAtt                                                                                                                      <  8:287539/92‑1 (MQ=255)
                                                                                 gAAAGCCACATCTTCGCTAAACAAAAATCGCTGGACGATCACTATGCCGAACTGAAAGCCGATGATGAAATCAGCGAACCCCTGCCACAAtt                                                                                                                      >  7:287539/1‑92 (MQ=255)
                                                                                       cacaGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGaa                                                                 >  1:324102/1‑139 (MQ=255)
                                                                                       cacaGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGaa                                                                 <  2:324102/139‑1 (MQ=255)
                                                                                        acaGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGcgccgccg                                                       >  2:301947/1‑148 (MQ=255)
                                                                                                ggTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGt                                               <  4:174094/148‑1 (MQ=255)
                                                                                                 gTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCTGCGTCCGAACGCGCCGCCGCTCATCGTCt                                             <  4:87014/149‑1 (MQ=255)
                                                                                                  tAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCt                                             >  5:337879/1‑148 (MQ=255)
                                                                                                  tAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCAGAACGCGCCGCCGCTCATCGACt                                             >  4:209877/1‑148 (MQ=255)
                                                                                                    aaCAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTaa                                           <  1:301947/148‑1 (MQ=255)
                                                                                                     aCAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCa                                                   <  3:209030/139‑1 (MQ=255)
                                                                                                     aCAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCa                                                   >  4:209030/1‑139 (MQ=255)
                                                                                                       aaaaaTCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGaa                                                                                                       >  4:68506/1‑85 (MQ=255)
                                                                                                       aaaaaTCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGaa                                                                                                       <  3:68506/85‑1 (MQ=255)
                                                                                                          aaTCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAg                                          >  8:134788/1‑143 (MQ=255)
                                                                                                          aaTCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAg                                          <  7:134788/143‑1 (MQ=255)
                                                                                                                tGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAg                                       <  5:247901/140‑1 (MQ=255)
                                                                                                                tGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAg                                       >  6:247901/1‑140 (MQ=255)
                                                                                                                  gACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTAt                                <  8:214393/145‑1 (MQ=255)
                                                                                                                  gACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTAt                                >  7:214393/1‑145 (MQ=255)
                                                                                                                          gTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAg                                                                             >  6:9137/1‑92 (MQ=255)
                                                                                                                          gTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAg                                                                             <  5:9137/92‑1 (MQ=255)
                                                                                                                           tttGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTAttt                   >  8:22900/1‑149 (MQ=255)
                                                                                                                              gCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTg                >  1:324468/1‑149 (MQ=255)
                                                                                                                                 gAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTgg               >  1:189812/1‑147 (MQ=255)
                                                                                                                                 gAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTgg               <  2:189812/147‑1 (MQ=255)
                                                                                                                                  aaCTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGc              <  1:145335/147‑1 (MQ=255)
                                                                                                                                  aaCTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGc              >  2:145335/1‑147 (MQ=255)
                                                                                                                                  aaCTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTa            >  2:56508/1‑149 (MQ=255)
                                                                                                                                   aCTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTAt           <  3:235065/149‑1 (MQ=255)
                                                                                                                                    cTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCAACGACTAAGCAGTAATTATGAGCGCGCTATTTCTGGCTAtt          >  5:67330/1‑149 (MQ=255)
                                                                                                                                     tGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTc         >  5:183621/1‑149 (MQ=255)
                                                                                                                                       aaaGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTCCg       >  5:37658/1‑149 (MQ=255)
                                                                                                                                         aGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTCCGt      >  6:211644/1‑148 (MQ=255)
                                                                                                                                             gatgatGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTCCGTTAAc  >  3:61345/1‑148 (MQ=255)
                                                                                                                                                gatgCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAgc                            >  3:173042/1‑119 (MQ=255)
                                                                                                                                                gatgCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAgc                            <  4:173042/119‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATAACATCCAGGCGGCGTCCGAACGCGCCGCCGCTCATCGTCTAAGGAGTACTTATGAGCGCGCTATTTCTGGCTATTCCGTTAAC  >  NC_000913/1368544‑1368832

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: