Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,791,886 C→A 7.2% T45T (ACG→ACT ydiV ← anti‑FlhDC factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,791,8860CA 7.2% 275.3 / 8.7 83T45T (ACG→ACTydiVanti‑FlhDC factor
Reads supporting (aligned to +/- strand):  ref base C (37/40);  new base A (4/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

AATAATAATGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCT  >  NC_000913/1791737‑1792035
                                                                                                                                                     |                                                                                                                                                     
aataataatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTACCAATGTTGATACTCTGTAAGTTGCGCTATAACCCGAAtc                                                                                                                                                       >  2:1616038/1‑150 (MQ=255)
  taataatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACtcgt                                                                                                                                                     >  1:1063807/1‑150 (MQ=255)
     taatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCgg                                                                                                                                                  <  1:895466/150‑1 (MQ=255)
      aatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCggg                                                                                                                                                 >  1:897883/1‑150 (MQ=255)
      aatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCggg                                                                                                                                                 <  1:1952184/150‑1 (MQ=255)
      aatGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCggg                                                                                                                                                 <  1:1568993/150‑1 (MQ=255)
        tGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGAt                                                                                                                                               <  2:1963895/150‑1 (MQ=255)
        tGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACGAGTCgggtt                                                                                                                                               <  1:270475/150‑3 (MQ=255)
           aGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACg                                                                                                                                            <  1:132760/150‑1 (MQ=255)
           aGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACGAGTCGGGTTata                                                                                                                                            <  1:778947/150‑3 (MQ=255)
            gCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGa                                                                                                                                           >  1:293537/1‑150 (MQ=255)
             ccACTTGAGGCGTTAAATTAAGCCAGGAAAAGAGTTTATGTTCAATAAAAAAATCTTGGCATGAATTTACTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGaa                                                                                                                                          >  2:1622581/1‑150 (MQ=255)
              cACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAAc                                                                                                                                         >  2:1689448/1‑150 (MQ=255)
              cACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAAc                                                                                                                                         <  1:38243/150‑1 (MQ=255)
              cACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAAc                                                                                                                                         <  1:103825/150‑1 (MQ=255)
              cACTTGTGGAGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACACGACTCGTCGGGATACGAAc                                                                                                                                         >  2:517912/1‑150 (MQ=255)
                  tgtgGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTa                                                                                                                                     <  1:1681571/150‑1 (MQ=255)
                    tgGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAAcc                                                                                                                                        <  2:389264/145‑1 (MQ=255)
                    tgGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAAcc                                                                                                                                        >  1:389264/1‑145 (MQ=255)
                     ggCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCa                                                                                                                                  <  1:747435/150‑1 (MQ=255)
                      gCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCAt                                                                                                                                 <  1:1499399/150‑1 (MQ=255)
                            aaTTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCAc                                                                                                                           <  1:779959/150‑1 (MQ=255)
                            aaTTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCAc                                                                                                                           >  2:215888/1‑150 (MQ=255)
                               tAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTgg                                                                                                                        <  1:1387401/150‑1 (MQ=255)
                                 aGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGgag                                                                                                                      <  1:1933295/150‑1 (MQ=255)
                                  gCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGgaga                                                                                                                     <  1:319504/150‑1 (MQ=255)
                                      ggCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATg                                                                                                                 >  1:1020710/1‑150 (MQ=255)
                                        cAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATggg                                                                                                               >  2:1949095/1‑150 (MQ=255)
                                        cAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATggg                                                                                                               >  2:1611826/1‑150 (MQ=255)
                                          aaGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTAcc                                                                                                                                   <  2:465104/128‑1 (MQ=255)
                                          aaGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTAcc                                                                                                                                   >  1:465104/1‑128 (MQ=255)
                                          aaGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTa                                                                                                             <  2:1763568/150‑1 (MQ=255)
                                              gTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCGATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTa                                                                                                         >  2:134422/1‑150 (MQ=255)
                                              gTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTa                                                                                                         <  2:347119/150‑1 (MQ=255)
                                                    gTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGATCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCaa                                                                                                   >  2:1540099/1‑150 (MQ=255)
                                                         aTAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCg                                                                                              <  1:584913/150‑1 (MQ=255)
                                                            aaaaaaTGTTGGCATGATTTTAGTAATTCTAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACa                                                                                           >  1:1659611/1‑150 (MQ=255)
                                                            aaaaaaTGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACa                                                                                           >  2:1941656/1‑150 (MQ=255)
                                                            aaaaaaTGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACa                                                                                           >  2:368467/1‑150 (MQ=255)
                                                            aaaaaaTGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACa                                                                                           >  2:1164136/1‑150 (MQ=255)
                                                                aaTGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCAATGGAGAAATGGGTAATTAACTCAACGCCGACAAGGa                                                                                       <  1:517912/150‑1 (MQ=255)
                                                                       gCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTg                                                                                >  1:1338256/1‑150 (MQ=255)
                                                                        cATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGTACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTgg                                                                               <  2:897883/150‑1 (MQ=255)
                                                                        cATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTACTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGTTACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGAATTGCTgg                                                                               >  2:1420458/1‑150 (MQ=255)
                                                                          tGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGtt                                                                             >  2:1195539/1‑150 (MQ=255)
                                                                           gATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTg                                                                            <  1:844656/150‑1 (MQ=255)
                                                                            aTTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGCCTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGt                                                                           >  2:855113/1‑150 (MQ=255)
                                                                             ttttAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGTAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGtc                                                                          <  2:293537/150‑1 (MQ=255)
                                                                             ttttAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGtc                                                                          <  1:368467/150‑1 (MQ=255)
                                                                                tAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTctct                                                                       <  1:215888/150‑1 (MQ=255)
                                                                                     aTTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCg                                                                  <  2:880299/150‑1 (MQ=255)
                                                                                       tCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGc                                                                >  1:171917/1‑150 (MQ=255)
                                                                                          aaCTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAgg                                                             >  1:109940/1‑150 (MQ=255)
                                                                                              gCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGCACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAAt                                                         >  2:1164284/1‑150 (MQ=255)
                                                                                              gCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAAt                                                         <  1:1689448/150‑1 (MQ=255)
                                                                                               ctctGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATa                                                        >  1:1782618/1‑150 (MQ=255)
                                                                                                tctGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAg                                                       >  1:1904324/1‑150 (MQ=255)
                                                                                                 ctGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCTGGATACGATCCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACTCCTACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGc                                                      <  2:594755/150‑1 (MQ=255)
                                                                                                 ctGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGc                                                      <  1:609902/150‑1 (MQ=255)
                                                                                                 ctGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGc                                                      <  1:1949095/150‑1 (MQ=255)
                                                                                                   gAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCaa                                                    >  2:1221861/1‑150 (MQ=255)
                                                                                                    aaaaTAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAAt                                                   >  1:1826467/1‑150 (MQ=255)
                                                                                                       ataatTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAg                                                <  2:633960/150‑1 (MQ=255)
                                                                                                          atTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                                             <  2:1904324/150‑1 (MQ=255)
                                                                                                            tGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGa                                           >  1:793790/1‑150 (MQ=255)
                                                                                                                 aTGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAg                                      <  2:1782618/150‑1 (MQ=255)
                                                                                                                  tGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                                     <  2:1338256/150‑1 (MQ=255)
                                                                                                                     tGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATtc                                  >  1:1883376/1‑150 (MQ=255)
                                                                                                                      gCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATtct                                 <  1:1195539/150‑1 (MQ=255)
                                                                                                                       ctcctcGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATtctc                                >  2:1962233/1‑150 (MQ=255)
                                                                                                                        tcctcGGTAAGTTTCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGATAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTcc                               <  1:1616038/150‑1 (MQ=255)
                                                                                                                        tcctcGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTcc                               <  1:1164136/150‑1 (MQ=255)
                                                                                                                            cGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaa                           >  1:1153230/1‑150 (MQ=255)
                                                                                                                                  gTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCt                     <  1:665615/150‑1 (MQ=255)
                                                                                                                                       gcTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATg                <  1:1941656/150‑1 (MQ=255)
                                                                                                                                          gtatcccGACGAGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa             >  2:243707/5‑150 (MQ=255)
                                                                                                                                          gtatcccGACGAGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa             >  1:229081/5‑150 (MQ=255)
                                                                                                                                          gtatcccGACGAGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa             >  1:1427363/5‑150 (MQ=255)
                                                                                                                                          ctatcccGACGAGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa             >  1:1123560/5‑150 (MQ=255)
                                                                                                                                           taACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTaa            >  2:188539/1‑150 (MQ=255)
                                                                                                                                           taACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAACATTCTCCAAAAAAATCTTCATGGTaa            >  1:732420/1‑150 (MQ=255)
                                                                                                                                            aaCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTGCAAAAAAATCTTCATGGTAAc           >  1:850998/1‑150 (MQ=255)
                                                                                                                                             aCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATc                                                                   >  2:572411/1‑93 (MQ=255)
                                                                                                                                             aCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATc                                                                   <  1:572411/93‑1 (MQ=255)
                                                                                                                                                 gACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTcc      >  2:1485339/1‑150 (MQ=255)
                                                                                                                                                     cgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct  >  2:258555/1‑150 (MQ=255)
                                                                                                                                                     cgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct  >  1:560909/1‑150 (MQ=255)
                                                                                                                                                     cgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct  >  2:85846/1‑150 (MQ=255)
                                                                                                                                                     cgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct  >  1:1864226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATAATAATGTAGCCACTTGTGGCGTTAAATTAAGCCAGGCAAAGAGTTTATGTTGAATAAAAAAATGTTGGCATGATTTTAGTAATTCCAACTGCTCTGAAAATAATTGCCAATGTTGCTCCTCGGTAAGTTGCGCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCT  >  NC_000913/1791737‑1792035

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