Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,391,330 G→A 5.0% D60D (GAC→GAT rsgA ← ribosome small subunit‑dependent GTPase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,391,3300GA 5.0% 323.7 / 4.1 99D60D (GAC→GATrsgAribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base G (48/45);  new base A (3/2);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

CGCGTCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGCGTTCACGCGGC  >  NC_000913/4391184‑4391475
                                                                                                                                                  |                                                                                                                                                 
cgcgTCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTc                                                                                                                                                 <  6:359921/149‑1 (MQ=255)
cgcgGCAACACCGAGGTACGCGCATGCACCGCTTCCACGAGCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACggc                                                                                                                                                 <  4:91844/149‑3 (MQ=255)
   gTCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGc                                                                                                                                               >  5:30905/1‑148 (MQ=255)
    tCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCa                                                                                                                                             >  4:187091/1‑149 (MQ=255)
    tCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCa                                                                                                                                             >  4:293393/1‑149 (MQ=255)
      aaCACCGAGGTACGCTCATGCACCGCATCCACGATCCCTTTGACAATTACGCCTACCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATc                                                                                                                                           >  8:21262/1‑149 (MQ=255)
        cacCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATcgg                                                                                                                                         <  5:145155/149‑1 (MQ=255)
         acCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATcggc                                                                                                                                        <  4:42614/149‑1 (MQ=255)
         acCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTAGCCATcggc                                                                                                                                        >  2:36725/1‑149 (MQ=255)
            gAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGa                                                                                                                                     <  7:235547/149‑1 (MQ=255)
                   gCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCcaca                                                                                                                              >  4:334610/1‑149 (MQ=255)
                     tCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACAt                                                                                                                             >  8:27439/1‑148 (MQ=255)
                       aTGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAg                                                                                                                          >  6:270325/1‑149 (MQ=255)
                         gCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCg                                                                                                                        >  8:206288/1‑149 (MQ=255)
                           aCCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATcggc                                                                                                                                        >  5:3526/1‑131 (MQ=255)
                           aCCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATcggc                                                                                                                                        <  6:3526/131‑1 (MQ=255)
                               cTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATa                                                                                                                  <  5:78801/149‑1 (MQ=255)
                                  ccACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATAc                                                                                                                 <  1:172474/147‑1 (MQ=255)
                                   cACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCa                                                                                                               >  8:42848/1‑148 (MQ=255)
                                    aCGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACATCGCCATCGGCGGATTCCACATCAGCGTGCATACCaaa                                                                                                             >  8:280434/1‑149 (MQ=255)
                                    aCGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACATCGCCATCGGCGGATTCCACATCAGCGTGCATACCa                                                                                                               >  6:133077/1‑147 (MQ=255)
                                      gATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGc                                                                                                           >  5:181312/1‑149 (MQ=255)
                                      gATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGc                                                                                                           >  3:253332/1‑149 (MQ=255)
                                       aTCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGCTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCa                                                                                                                           >  5:142166/1‑132 (MQ=255)
                                       aTCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGCTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCa                                                                                                                           <  6:142166/132‑1 (MQ=255)
                                       aTCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTc                                                                                                                                                 <  4:347383/110‑1 (MQ=255)
                                       aTCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTc                                                                                                                                                 >  3:347383/1‑110 (MQ=255)
                                          ccTTTGACATTTACGCCTTCCGCCGCCTGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCt                                                                                                       <  1:172939/149‑1 (MQ=255)
                                           cTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTTGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCt                                                                                                       >  1:262738/1‑148 (MQ=255)
                                           cTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCt                                                                                                       >  1:167712/1‑148 (MQ=255)
                                             ttGACATTTACGCCTTCCGCCGCCTGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCCAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCtgat                                                                                                    >  1:371911/1‑149 (MQ=255)
                                             ttGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGACCGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCtgat                                                                                                    >  5:314406/1‑149 (MQ=255)
                                              tGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCtgat                                                                                                    <  6:266780/148‑1 (MQ=255)
                                                  aTTTCCGCCTTCCGCCGCCGGTGTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCACAGCGGCTGATGACGa                                                                                               <  6:409951/149‑1 (MQ=255)
                                                  aTTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGa                                                                                               >  1:195421/1‑149 (MQ=255)
                                                   tttACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCTCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACAt                                                                                                                             <  5:164394/118‑1 (MQ=255)
                                                   tttACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCTCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACAt                                                                                                                             >  6:164394/1‑118 (MQ=255)
                                                   tttACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCtga                                                                                                     >  6:359724/1‑142 (MQ=255)
                                                   tttACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCtga                                                                                                     <  5:359724/142‑1 (MQ=255)
                                                       cGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGGCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGGCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                                                           <  8:404741/148‑1 (MQ=255)
                                                         ccTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACATCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                                                        <  5:133077/149‑1 (MQ=255)
                                                         ccTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACATCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                                                        <  7:280434/149‑1 (MQ=255)
                                                          cTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                                                       >  7:113943/1‑149 (MQ=255)
                                                          cTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAGCGGCTGATGCCGCTACCTtcat                                                                                      >  6:222832/1‑149 (MQ=255)
                                                           ttCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGAc                                                                                                 >  4:246065/1‑138 (MQ=255)
                                                           ttCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGAc                                                                                                 <  3:246065/138‑1 (MQ=255)
                                                           ttCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                                                      >  8:49531/1‑149 (MQ=255)
                                                            tCCGCCGCCGGTTTACCCGGACGGCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGACCGTCGCCATCGGCGGATGCCACATCCGCGTGCGTACCAAAGCGGCTGATGACGCTACCGtcct                                                                                      >  5:109622/1‑146 (MQ=255)
                                                              cgccgccgGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                                                   >  2:217533/1‑149 (MQ=255)
                                                                ccgccgGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                                                 <  8:288992/149‑1 (MQ=255)
                                                                   ccgGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCg                                                                              <  3:187091/149‑1 (MQ=255)
                                                                      gTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCa                                                                           <  4:96396/149‑1 (MQ=255)
                                                                            ccGGACGCCAGACTACGCGGTCGCCGGTTACCAGAGAACGGATGGTACGGCGAATATTGCAGCTGTGAACTTCGACATCGGAGGAGTCCACATCAGCGTTCATACAAACGAGGCTGATGACTATAAATTCATCAGGATCGCCaaa                                                                         >  5:281398/1‑145 (MQ=255)
                                                                              ggACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                                                   >  8:115038/1‑149 (MQ=255)
                                                                               gACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                                                  >  8:307042/1‑149 (MQ=255)
                                                                               gACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                                                  <  2:174803/149‑1 (MQ=255)
                                                                               gACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTAAGGCGAATATTGCAGCGGTGAACGTCGGCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                                                  <  8:73427/149‑1 (MQ=255)
                                                                                 cGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                                                                <  5:301775/149‑1 (MQ=255)
                                                                                    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                                                             <  7:42848/149‑1 (MQ=255)
                                                                                     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATAcc                                                                                                                <  6:385384/97‑1 (MQ=255)
                                                                                     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGtcgtcg                                                            <  3:334610/149‑1 (MQ=255)
                                                                                        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGc                                                                                                                     >  5:285714/1‑89 (MQ=255)
                                                                                        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGc                                                                                                                     <  6:285714/89‑1 (MQ=255)
                                                                                        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAg                                                         <  7:49531/149‑1 (MQ=255)
                                                                                         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                                                        >  5:271298/1‑149 (MQ=255)
                                                                                           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTTGTAGTCg                                                      >  1:283732/1‑149 (MQ=255)
                                                                                           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                                                       <  4:253332/148‑1 (MQ=255)
                                                                                           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                                                      <  4:134391/149‑1 (MQ=255)
                                                                                                  gCCGGTTTCCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGATCGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtct                                               <  7:21262/149‑1 (MQ=255)
                                                                                                  gCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtct                                               >  3:62236/1‑149 (MQ=255)
                                                                                                   ccGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtctc                                              >  2:122200/1‑149 (MQ=255)
                                                                                                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACATCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc                                             >  2:387317/1‑149 (MQ=255)
                                                                                                     ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt                                            <  3:293393/149‑1 (MQ=255)
                                                                                                       ttACCAGCGTACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCGACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa                                          <  6:30905/149‑1 (MQ=255)
                                                                                                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa                                          <  1:122200/149‑1 (MQ=255)
                                                                                                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCCGGCTTCTCCTTa                                          <  2:167712/149‑1 (MQ=255)
                                                                                                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAATCAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc                                       >  2:65299/1‑149 (MQ=255)
                                                                                                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGCACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACg                                      >  7:72986/1‑149 (MQ=255)
                                                                                                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACg                                      >  4:150644/1‑149 (MQ=255)
                                                                                                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt                                     >  5:142795/1‑149 (MQ=255)
                                                                                                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt                                     >  7:42403/1‑149 (MQ=255)
                                                                                                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGCTTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtt                                    >  5:84052/1‑149 (MQ=255)
                                                                                                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                                                             <  5:190929/124‑1 (MQ=255)
                                                                                                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                                                             >  6:190929/1‑124 (MQ=255)
                                                                                                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtt                                    <  6:142795/149‑1 (MQ=255)
                                                                                                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGACATCGGCGGATTCCACATCAGCGTGGATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCTGGCTTCTCCTTAGACGtttt                                  <  1:68500/149‑1 (MQ=255)
                                                                                                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTaa                                >  2:318345/1‑146 (MQ=255)
                                                                                                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTaa                                <  1:318345/146‑1 (MQ=255)
                                                                                                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCACAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGc                        <  5:187143/149‑1 (MQ=255)
                                                                                                                           cGGCGAATATTGCAGCGGTGAACGTTGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCt                       <  2:262738/148‑1 (MQ=255)
                                                                                                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTgg                     <  5:270325/149‑1 (MQ=255)
                                                                                                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCGGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTgg                     <  6:109622/149‑1 (MQ=255)
                                                                                                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCCTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTATAGACGACGCTgg                     >  7:281488/1‑149 (MQ=255)
                                                                                                                                aaTATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTGCTCCTTAGACGTTTTAAGACGACGCTGGTGAt                 <  8:113943/149‑1 (MQ=255)
                                                                                                                                 atatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGAtt                >  8:351907/1‑149 (MQ=255)
                                                                                                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGc             <  3:334856/149‑1 (MQ=255)
                                                                                                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGCg            >  1:371801/1‑149 (MQ=255)
                                                                                                                                       cAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGCGt           <  3:322832/148‑1 (MQ=255)
                                                                                                                                             tGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGAGTTCAcgcg    >  3:336093/1‑149 (MQ=255)
                                                                                                                                               aaCGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGCGTTCACgcggc  >  3:395189/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGCGTCAACACCGAGGTACGCTCATGCACCGCTTCCACGATCCCTTTGACATTTACGCCTTCCGCCGCCGGTTTACCCGGACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTTAAGACGACGCTGGTGATTGGCGTTCACGCGGC  >  NC_000913/4391184‑4391475

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: