Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,495 A→G S68G (AGC→GGC)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,4950AG99.1% 422.2 / ‑4.3 114S68G (AGC→GGC) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (49/64);  total (50/64)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

Only 100 of 114 total aligned reads displayed.

CTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTA  >  NC_000913/3937348‑3937642
                                                                                                                                                   |                                                                                                                                                   
cTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGc                                                                                                                                                   >  2:1418925/1‑150 (MQ=255)
 ttAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCAt                                                                                                                                                 >  1:400908/1‑151 (MQ=255)
  tAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGc                                                                                                                                                   <  2:586920/148‑1 (MQ=255)
   aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCAtt                                                                                                                                                <  1:101393/150‑1 (MQ=255)
   aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTg                                                                                                                                               >  2:1340145/1‑151 (MQ=255)
   aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTg                                                                                                                                               >  2:481121/1‑151 (MQ=255)
   aaTCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTg                                                                                                                                               <  1:1364814/151‑1 (MQ=255)
     tCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTgg                                                                                                                                              >  2:1451481/1‑150 (MQ=255)
      cTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTg                                                                                                                                            >  2:1410562/1‑151 (MQ=255)
        tCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTg                                                                                                                                            <  2:440422/149‑1 (MQ=255)
         cAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTg                                                                                                                                            >  2:940301/1‑148 (MQ=255)
          aaTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTgag                                                                                                                                          >  1:1007720/1‑149 (MQ=255)
          aaTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGagag                                                                                                                                        <  2:1022711/151‑1 (MQ=255)
               tttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGt                                                                                                                                     <  2:1288432/149‑1 (MQ=255)
                ttCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCg                                                                                                                                  >  2:1395425/1‑151 (MQ=255)
                 tCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGc                                                                                                                                 >  2:148621/1‑151 (MQ=255)
                        ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCcagcag                                                                                                                           <  1:458402/150‑1 (MQ=255)
                        ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCcagca                                                                                                                            <  2:1321049/149‑1 (MQ=255)
                        ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCcagca                                                                                                                            <  1:341836/149‑1 (MQ=255)
                        ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCcagca                                                                                                                            <  1:780309/149‑1 (MQ=255)
                        ccAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCagcagc                                                                                                                          >  1:585500/1‑151 (MQ=255)
                           ggCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCg                                                                                                                       <  1:708383/151‑1 (MQ=255)
                            gCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGc                                                                                                                      <  1:179843/151‑1 (MQ=255)
                              gAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGcc                                                                                                                     >  2:507209/1‑150 (MQ=255)
                              gAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGcc                                                                                                                     >  1:669412/1‑150 (MQ=255)
                               aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCAc                                                                                                                   >  1:743901/1‑151 (MQ=255)
                               aaaCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCAc                                                                                                                   >  1:848288/1‑151 (MQ=255)
                                 aCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTg                                                                                                                 <  1:1191221/151‑1 (MQ=255)
                                 aCCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTg                                                                                                                 <  2:1267080/151‑1 (MQ=255)
                                  ccGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCAAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                                                                                <  2:349017/151‑1 (MQ=255)
                                   cGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                                               >  2:711062/1‑151 (MQ=255)
                                   cGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                                               <  1:1138946/151‑1 (MQ=255)
                                   cGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                                               <  1:393490/151‑1 (MQ=255)
                                    gTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATa                                                                                                              <  2:1324266/151‑1 (MQ=255)
                                          ggTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTg                                                                                                        >  2:96936/1‑151 (MQ=255)
                                           gTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGa                                                                                                       <  1:1330997/151‑1 (MQ=255)
                                            tAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGat                                                                                                      >  2:1082559/1‑151 (MQ=255)
                                             aaCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                                                                                     <  1:791091/151‑1 (MQ=255)
                                                 aCTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTAc                                                                                                 <  1:462044/151‑1 (MQ=255)
                                                  cTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACt                                                                                                >  1:617124/1‑151 (MQ=255)
                                                   tATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTc                                                                                               >  1:982436/1‑151 (MQ=255)
                                                        ggTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                                                                          <  2:10789/151‑1 (MQ=255)
                                                        ggTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                                                                          <  2:607704/151‑1 (MQ=255)
                                                         gTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                                                                          >  2:895170/1‑150 (MQ=255)
                                                         gTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCa                                                                                         <  2:35946/151‑1 (MQ=255)
                                                          ttGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                                                                          >  1:916077/1‑149 (MQ=255)
                                                           tGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAg                                                                                        <  1:507267/150‑1 (MQ=255)
                                                             cATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGt                                                                                     >  2:308918/1‑151 (MQ=255)
                                                               tttGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTg                                                                                    >  2:116825/1‑150 (MQ=255)
                                                                ttGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGCGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGa                                                                                   <  2:92549/150‑1 (MQ=255)
                                                                 tGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATc                                                                                 >  2:917165/1‑151 (MQ=255)
                                                                  ggcggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATc                                                                                 <  1:22615/150‑1 (MQ=255)
                                                                     ggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCaaa                                                                              <  2:98814/150‑1 (MQ=255)
                                                                       cAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                                           >  1:1402677/1‑151 (MQ=255)
                                                                       cAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                                           <  2:59283/151‑1 (MQ=255)
                                                                        aaaGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCg                                                                          >  2:1107591/1‑151 (MQ=255)
                                                                        aaaGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCg                                                                          >  2:1255708/1‑151 (MQ=255)
                                                                          aGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCg                                                                          >  1:876636/1‑149 (MQ=255)
                                                                            gcgcGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATc                                                                      <  1:437903/151‑1 (MQ=255)
                                                                              gcgAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCaa                                                                    <  1:865851/151‑1 (MQ=255)
                                                                                gAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACa                                                                  <  1:1218176/151‑1 (MQ=255)
                                                                                 aaTCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAAc                                                                   >  2:1413617/1‑149 (MQ=255)
                                                                                 aaTCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAg                                                                 <  2:287165/151‑1 (MQ=255)
                                                                                    cAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtg                                                              >  2:728695/1‑151 (MQ=255)
                                                                                    cAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgt                                                               <  1:413141/150‑1 (MQ=255)
                                                                                          gtgGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGgcg                                                         <  2:682557/150‑1 (MQ=255)
                                                                                           tgGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCt                                                       >  2:318862/1‑151 (MQ=255)
                                                                                            ggCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGTTg                                                      <  2:1422501/151‑1 (MQ=255)
                                                                                               tgctgGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAtt                                                   >  1:612901/1‑151 (MQ=255)
                                                                                                gctgGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAt                                                    <  2:792163/149‑1 (MQ=255)
                                                                                                   gggCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTg                                               <  2:315740/151‑1 (MQ=255)
                                                                                                     gCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGtt                                             <  2:415928/151‑1 (MQ=255)
                                                                                                           cGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGc                                       <  1:1287689/151‑1 (MQ=255)
                                                                                                             gTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGa                                     <  2:1398362/151‑1 (MQ=255)
                                                                                                                  aaTATCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGAGTGGCGCTGATTTTTGTTAATGGCGAAGGTg                                >  1:303000/1‑151 (MQ=255)
                                                                                                                   atatCGCGTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTga                               <  1:739848/151‑1 (MQ=255)
                                                                                                                     atCGCGTTTATTGCCTGTACGGGGGATGACGGCATTGGGGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTaa                                           >  1:358154/1‑137 (MQ=255)
                                                                                                                       cgcgTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAAt                           <  1:791978/151‑1 (MQ=255)
                                                                                                                        gcgTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAAt                           <  1:906064/150‑1 (MQ=255)
                                                                                                                         cgTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGt                         >  2:1338755/1‑151 (MQ=255)
                                                                                                                         cgTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGt                         <  1:1116208/151‑1 (MQ=255)
                                                                                                                          gTTTATTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTc                        >  1:93153/1‑151 (MQ=255)
                                                                                                                              aTTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCAta                     <  2:1004156/150‑2 (MQ=255)
                                                                                                                              aTTGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCg                    >  1:353557/1‑151 (MQ=255)
                                                                                                                                tGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCgg                   <  1:1290335/150‑1 (MQ=255)
                                                                                                                                tGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCgg                   <  1:422885/150‑1 (MQ=255)
                                                                                                                                tGCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCgg                   <  1:116842/150‑1 (MQ=255)
                                                                                                                                 gCCTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTa                 <  1:1395532/151‑1 (MQ=255)
                                                                                                                                  ccTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTAt                >  1:1095021/1‑151 (MQ=255)
                                                                                                                                   cTGTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTAtt               >  2:551153/1‑151 (MQ=255)
                                                                                                                                     gTACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCa             <  2:1435549/151‑1 (MQ=255)
                                                                                                                                      tACGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCAt            >  1:937309/1‑151 (MQ=255)
                                                                                                                                       aCGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCAt            <  1:32595/150‑1 (MQ=255)
                                                                                                                                        cGGGTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGc          >  1:474091/1‑151 (MQ=255)
                                                                                                                                         gggTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGcc         >  1:1409703/1‑151 (MQ=255)
                                                                                                                                          ggTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCg        <  2:669340/151‑1 (MQ=255)
                                                                                                                                          ggTGATGACGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCg        <  1:68987/151‑1 (MQ=255)
                                                                                                                                                gaCGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCt   >  2:734114/1‑150 (MQ=255)
                                                                                                                                                gaCGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCt   <  1:481176/150‑1 (MQ=255)
                                                                                                                                                gaCGGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTa  <  1:968579/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTA  >  NC_000913/3937348‑3937642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑