Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,892,358 C→T 15.8% G293D (GGT→GAT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,3580CT15.8% 326.5 / 41.9 132G293D (GGT→GAT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base C (53/58);  new base T (8/13);  total (61/71)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 139 total aligned reads displayed.

AGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGG  >  NC_000913/1892259‑1892457
                                                                                                   |                                                                                                   
aGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTAc                                                                                                     <  1:951976/100‑1 (MQ=255)
   gACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                  <  2:2404416/100‑1 (MQ=255)
   gACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                  <  1:2999505/100‑1 (MQ=255)
      aGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                                               <  1:2579517/100‑1 (MQ=255)
       gTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATa                                                                                              >  2:130315/1‑100 (MQ=255)
            gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                                                         >  2:3306057/1‑100 (MQ=255)
            gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                                                         >  2:3259583/1‑100 (MQ=255)
             ccACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                                                        <  1:1446781/100‑1 (MQ=255)
             ccACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                                                        <  2:3245570/100‑1 (MQ=255)
             ccACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCACGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                                                        >  2:1173362/1‑100 (MQ=255)
              cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcg                                                                                       <  1:904302/100‑3 (MQ=255)
              cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcg                                                                                       <  2:4143452/100‑3 (MQ=255)
              cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                                                       <  1:3793393/100‑1 (MQ=255)
               aCGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgt                                                                                      <  2:1236162/100‑4 (MQ=255)
               aCGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCt                                                                                      <  1:4109118/100‑1 (MQ=255)
                cGTCATAACTACGTTCCAAGCTGTCATCGAGCAATAACAATGCCCGCTCAATCTCCGGATTAGCTAACAGCGCACGCAGGTTATCGCGGTAACCTgcgtg                                                                                     <  2:1984345/100‑5 (MQ=255)
                cGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtg                                                                                     <  2:4006056/100‑5 (MQ=255)
                cGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtg                                                                                     <  1:1312266/100‑5 (MQ=255)
                  tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgag                                                                                   <  1:750153/100‑7 (MQ=39)
                  tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgag                                                                                   <  1:1030740/100‑7 (MQ=39)
                  tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgag                                                                                   <  2:3046963/100‑7 (MQ=39)
                  tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                                                   <  2:1084559/100‑1 (MQ=255)
                   cATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgagc                                                                                  <  1:3246031/100‑8 (MQ=38)
                    aTAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCg                                                                                 >  2:2606938/1‑100 (MQ=255)
                     tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGa                                                                                >  2:2719157/1‑100 (MQ=255)
                     tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGa                                                                                >  1:2045636/1‑100 (MQ=255)
                      aaCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGa                                                                                >  2:1685772/1‑99 (MQ=255)
                      aaCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAg                                                                               >  1:3968020/1‑100 (MQ=255)
                        cAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCCCGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                                                              >  1:1116817/1‑99 (MQ=255)
                        cAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGtt                                                                             >  2:2327867/1‑100 (MQ=255)
                          aaGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGc                                                                           <  2:157659/100‑1 (MQ=255)
                            gTTCCAGGGTGTCAACGAGCAGTAAACATGCCCGCTGAATTTGCGGATTAGCTAACAGATCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTg                                                                            >  2:4111636/1‑97 (MQ=255)
                             ttCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                                                        <  2:1109525/100‑1 (MQ=255)
                              tCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                                                       >  2:58631/1‑100 (MQ=255)
                                cAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGa                                                                     >  2:1404409/1‑100 (MQ=255)
                                    gtgtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCGCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAAt                                                                 <  2:2326425/100‑1 (MQ=255)
                                      gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                                               <  2:2593745/100‑1 (MQ=255)
                                      gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                                               >  1:3382518/1‑100 (MQ=255)
                                      gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                                               <  2:491129/100‑1 (MQ=255)
                                      gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTACAGACGAAAATcc                                                               >  2:2274027/1‑100 (MQ=255)
                                       tCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCt                                                              <  2:467142/100‑1 (MQ=255)
                                        cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGCGTTGCAGACGAAAATCCt                                                              >  2:895602/1‑99 (MQ=255)
                                        cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                             >  2:3684835/1‑100 (MQ=255)
                                        cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                             >  1:2401417/1‑100 (MQ=255)
                                         aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                            <  1:3099388/100‑1 (MQ=255)
                                         aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGATCACGCAGGTTACCGCGATAAACTGGCTCACCGAGTTGCAGACGAAAATc                                                                >  2:2162250/1‑96 (MQ=255)
                                          tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcg                                                           <  2:424089/100‑1 (MQ=255)
                                          tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcg                                                           <  1:2151907/100‑1 (MQ=255)
                                           ctaGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcgc                                                          <  2:3814043/98‑1 (MQ=255)
                                           cGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgcgc                                                          >  1:1730133/1‑100 (MQ=255)
                                               cAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActc                                                      >  1:4059373/1‑100 (MQ=255)
                                                aGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACTGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActct                                                     >  2:1484191/1‑100 (MQ=255)
                                                 gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                                                    >  1:2456574/1‑100 (MQ=255)
                                                  tccgatctcccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTgg                                                   <  1:3405695/92‑1 (MQ=255)
                                                  tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTgg                                                   >  2:1313228/1‑100 (MQ=255)
                                                    aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGc                                                 <  2:211982/100‑1 (MQ=255)
                                                       aTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGCTCCCCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGa                                             >  2:1291960/1‑100 (MQ=255)
                                                          cccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTggagatcgga                                           >  2:3405679/1‑92 (MQ=255)
                                                           ccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGa                                          <  1:3622933/100‑1 (MQ=255)
                                                             gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATc                                        <  1:1664379/100‑1 (MQ=255)
                                                             gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATc                                        <  1:4069938/100‑1 (MQ=255)
                                                              cTGAAATTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                       >  2:2490818/1‑100 (MQ=255)
                                                                  aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcgc                                   >  1:792099/1‑100 (MQ=255)
                                                                   tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcgc                                   <  1:1251110/99‑1 (MQ=255)
                                                                   tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCa                                  >  1:275937/1‑100 (MQ=255)
                                                                    ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAc                                 <  1:3234394/100‑1 (MQ=255)
                                                                       cGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                              <  1:302730/100‑1 (MQ=255)
                                                                        ggATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACtttt                             >  2:729891/1‑100 (MQ=255)
                                                                         gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACtttgg                            <  1:2845931/100‑3 (MQ=255)
                                                                         gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTg                            <  1:2727006/100‑1 (MQ=255)
                                                                          aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGt                           <  1:1199432/100‑1 (MQ=255)
                                                                             aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGTGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                        >  1:3112246/1‑100 (MQ=255)
                                                                             aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                        <  2:795658/100‑1 (MQ=255)
                                                                                tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActgc                     >  2:259789/1‑100 (MQ=255)
                                                                                tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCAATCTGGGCAAGACGATCAGCGCACTTTTGTAActgc                     >  1:543391/1‑100 (MQ=255)
                                                                                   cAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACTATCAGCGCACTTTTGTAACTGCTgg                  <  1:259789/100‑1 (MQ=255)
                                                                                   cAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTgg                  <  1:58631/100‑1 (MQ=255)
                                                                                     gCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgt                >  1:1514074/1‑100 (MQ=255)
                                                                                     gCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgt                <  2:3463227/100‑1 (MQ=255)
                                                                                      cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAGTTTTGTAACTGCTGGgt                >  1:3663356/1‑99 (MQ=255)
                                                                                       tgagccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtgt              >  2:2405664/6‑100 (MQ=255)
                                                                                        cACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGgtg               >  1:806932/1‑98 (MQ=255)
                                                                                        cACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTc             >  1:2597914/1‑100 (MQ=255)
                                                                                         agccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTTCCCACTCCGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTc             >  2:1066325/4‑99 (MQ=255)
                                                                                         aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCt            >  1:3095925/1‑100 (MQ=255)
                                                                                         aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCt            <  1:1186373/100‑1 (MQ=255)
                                                                                          gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtt           >  2:2527518/3‑100 (MQ=255)
                                                                                          gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtt           >  1:1397258/3‑100 (MQ=255)
                                                                                          gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCt            >  2:3418141/3‑99 (MQ=255)
                                                                                          cgaaGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTCCGCACTCTGGGCAAGACGAACAGCGCACTTTTGGAAcggctg                   >  2:4055829/4‑92 (MQ=255)
                                                                                          cGCAAGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtt           <  1:2234299/100‑1 (MQ=255)
                                                                                           gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt          <  1:3684856/100‑1 (MQ=255)
                                                                                           gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt          <  2:659585/100‑1 (MQ=255)
                                                                                            cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtttt         >  2:3987993/1‑100 (MQ=255)
                                                                                             aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCACGACGATCAGCGCACTTTTGTAACTTCTGGGTGTCtttt         <  2:2588085/99‑1 (MQ=255)
                                                                                               gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAAt      <  2:506912/100‑1 (MQ=255)
                                                                                                 tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTc    <  1:1095604/100‑1 (MQ=255)
                                                                                                   ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTTCGCACTCTGCGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCg   >  2:620111/1‑99 (MQ=255)
                                                                                                   ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg  <  1:1464252/100‑1 (MQ=255)
                                                                                                   ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg  <  1:3306072/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGG  >  NC_000913/1892259‑1892457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: