Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,896,050 A→C N416H (AAC→CAC)  pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,896,0500AC97.1% 588.7 / ‑1.1 173N416H (AAC→CAC) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base C (85/83);  total (86/87)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

CCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCC  >  NC_000913/1895954‑1896146
                                                                                                |                                                                                                
ccGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACg                                                                                               <  1:4555114/100‑1 (MQ=255)
ccGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACg                                                                                               >  1:5425359/1‑100 (MQ=255)
ccGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACg                                                                                               <  1:1691936/100‑1 (MQ=255)
      cGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACaaa                                                                                         <  2:708016/100‑1 (MQ=255)
       gACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAAt                                                                                        >  2:1107433/1‑100 (MQ=255)
         cGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAAttt                                                                                      <  1:839873/100‑1 (MQ=255)
           cAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtc                                                                                    <  1:1934291/100‑1 (MQ=255)
              tGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTGGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctgc                                                                                 <  2:1129283/100‑1 (MQ=255)
               gCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGACAAATTTtctgct                                                                                <  2:1712979/100‑1 (MQ=255)
                   ggTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctgctctgc                                                                            <  1:4623986/100‑1 (MQ=255)
                    gTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcg                                                                           <  1:2383486/100‑1 (MQ=255)
                    gTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcg                                                                           <  2:501791/100‑1 (MQ=255)
                     tGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgg                                                                          >  2:1317220/1‑100 (MQ=255)
                      gCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTTCTATCCGCACGCTGACTGCCGTTCACGGACAAATTTTCTGCTCTgcggg                                                                         <  2:5280876/100‑1 (MQ=255)
                       cGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggc                                                                        >  2:2284554/1‑100 (MQ=255)
                        ggCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggcg                                                                       >  2:3179077/1‑100 (MQ=255)
                         gcagcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggcgg                                                                      >  1:3268302/1‑100 (MQ=255)
                           agcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTg                                                                    <  2:2092782/100‑1 (MQ=255)
                            gcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTgg                                                                   <  2:5455943/100‑1 (MQ=255)
                            gcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTgg                                                                   <  2:4351772/100‑1 (MQ=255)
                             caTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGa                                                                  <  2:3498483/100‑1 (MQ=255)
                             caTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGa                                                                  <  1:3601882/100‑1 (MQ=255)
                               ttGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAt                                                                >  2:1300790/1‑100 (MQ=255)
                                tGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAtt                                                               >  2:325276/1‑100 (MQ=255)
                                tGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAtt                                                               <  1:485163/100‑1 (MQ=255)
                                 ggCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTg                                                              >  2:90587/1‑100 (MQ=255)
                                 ggCTATTTGAGCTTTGGCGGCAACATGGATGCCAGTATTACTATCAGCACGCTGACTCCCATTCACGGACAAATTTACTGCTCTGCGGGCAGTGGAATTg                                                              >  1:4844525/1‑100 (MQ=255)
                                   cTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTc                                                            >  2:399929/1‑100 (MQ=255)
                                    tATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTTTGCGGGCGGTGGAATTGTCg                                                           <  2:2236799/100‑1 (MQ=255)
                                    tATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCg                                                           >  1:1064759/1‑100 (MQ=255)
                                    tATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCg                                                           <  2:88252/100‑1 (MQ=255)
                                     aTTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGc                                                          >  1:85142/1‑100 (MQ=255)
                                      tttGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGcc                                                         <  1:1028848/100‑1 (MQ=255)
                                       ttGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCg                                                        >  1:990210/1‑100 (MQ=255)
                                       ttGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCg                                                        <  2:1417955/100‑1 (MQ=255)
                                       ttGAGCTTTTGCGGAAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCg                                                        <  2:570888/100‑1 (MQ=255)
                                          aGCTTTTGAGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATa                                                     >  2:5266531/1‑100 (MQ=255)
                                           gCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAg                                                    <  2:5352411/100‑1 (MQ=255)
                                            cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGTACAAATTTTCTGCCCTGCGGGCGGTGGAATTGTCGCCGATAGc                                                   >  2:5422280/1‑100 (MQ=255)
                                            cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGc                                                   >  1:4756562/1‑100 (MQ=255)
                                              tttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCa                                                 >  1:3324407/1‑100 (MQ=255)
                                              tttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCa                                                 <  1:2108369/100‑1 (MQ=255)
                                              tttGCGGCAACATGGATACCAGTACTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCa                                                 <  1:5012377/100‑1 (MQ=255)
                                                 gCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGga                                              <  2:4461957/100‑1 (MQ=255)
                                                  cGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAGATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaa                                             >  2:5146001/1‑100 (MQ=255)
                                                  cGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaa                                             <  2:2165167/100‑1 (MQ=255)
                                                  cGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaa                                             <  1:416089/100‑1 (MQ=255)
                                                  cGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaa                                             >  1:1465425/1‑100 (MQ=255)
                                                   ggCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaag                                            >  2:2884096/1‑100 (MQ=255)
                                                   ggCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaag                                            <  1:4321925/100‑1 (MQ=255)
                                                    gCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaaga                                           >  1:3621915/1‑100 (MQ=255)
                                                     cAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaagaa                                          <  2:4743947/100‑1 (MQ=255)
                                                      aaCATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGaagaag                                         <  1:2363724/100‑1 (MQ=255)
                                                      aaCATGGATACCAGTATTACTATCCGCACGCTGACGGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGatgaag                                         <  2:4640536/100‑1 (MQ=255)
                                                        cATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCg                                       >  1:4261891/1‑100 (MQ=255)
                                                          tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGTGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCggn                                     <  1:532923/100‑2 (MQ=255)
                                                          tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     >  2:34454/1‑100 (MQ=255)
                                                          tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     <  2:4823641/100‑1 (MQ=255)
                                                            gttACCAGTAATACTATCCGCACGGTGAATGCCATTCACGGACAAATTTTCTGCACGGCGCACGGGGGAATTGTCGCCGATGTCCAGTAAGGAGCGGAat                                   >  2:1909537/3‑100 (MQ=255)
                                                            gATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAat                                   <  1:4567540/100‑1 (MQ=255)
                                                              tACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAatat                                 <  1:3714940/100‑1 (MQ=255)
                                                                ccAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCa                               >  1:2154756/1‑100 (MQ=255)
                                                                ccAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCa                               >  1:934322/1‑100 (MQ=255)
                                                                 cAGTATTACTATCCGCACGCTGGCTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAg                              <  1:1370896/100‑1 (MQ=255)
                                                                 cAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAg                              >  1:490425/1‑100 (MQ=255)
                                                                  aGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAgg                             >  2:4883654/1‑100 (MQ=255)
                                                                  aGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAgg                             >  1:260559/1‑100 (MQ=255)
                                                                   gTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGa                            >  1:2220361/1‑100 (MQ=255)
                                                                    tATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGaa                           <  2:1194729/100‑1 (MQ=255)
                                                                       tACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACt                        >  1:2097193/1‑100 (MQ=255)
                                                                       tACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACt                        <  2:4692165/100‑1 (MQ=255)
                                                                        aCTATCCTCACGCTTACGGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtt                       >  2:2771496/1‑100 (MQ=255)
                                                                        aCTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtt                       >  1:391278/1‑100 (MQ=255)
                                                                         cTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACttt                      >  2:4185757/1‑100 (MQ=255)
                                                                          tATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtttt                     <  1:4576920/100‑1 (MQ=255)
                                                                          tATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtttt                     <  1:2657921/100‑1 (MQ=255)
                                                                             ccGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                  <  2:3381918/100‑1 (MQ=255)
                                                                              cGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATa                 <  2:2357085/100‑1 (MQ=255)
                                                                              cGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATa                 >  2:4822141/1‑100 (MQ=255)
                                                                              cGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATa                 >  1:478509/1‑100 (MQ=255)
                                                                               gCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATaa                >  2:25668/1‑100 (MQ=255)
                                                                               gCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATaa                >  2:4084085/1‑100 (MQ=255)
                                                                               gCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATaa                >  2:3287015/1‑100 (MQ=255)
                                                                                cACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATaaa               <  1:5274653/100‑1 (MQ=255)
                                                                                cACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATaaa               <  1:1346863/100‑1 (MQ=255)
                                                                                 aCGCTGACTGCCATTAACGGAGAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAg              <  2:3327332/100‑1 (MQ=255)
                                                                                  cGCTGACTGCCATTCACGGACCAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGGCAGGAAGAAGCGGAATATCGGGAACCTTTTGATAAAGt             <  2:3089041/100‑1 (MQ=255)
                                                                                   gCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGtt            <  2:1925460/100‑1 (MQ=255)
                                                                                   gCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGtt            <  2:4777455/100‑1 (MQ=255)
                                                                                   gCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGtt            <  1:2011211/100‑1 (MQ=255)
                                                                                     tGACCGTCATTCACGGACACAGAGTCTGCTATGCGGTCTGGCGAATTGTCGCTTATAGCTACGACCGATCGTAATATCAGGAAACTTTTGATAAAGTTaa          <  2:4174232/100‑1 (MQ=255)
                                                                                      gACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAAt         <  2:5335017/100‑1 (MQ=255)
                                                                                        cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGGAGAAGCGGAATAGCAGGAAACTGTGGATAAAGGTAAGCg       >  1:1648655/1‑100 (MQ=255)
                                                                                        cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAATAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCg       <  2:3660512/100‑1 (MQ=255)
                                                                                        cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCg       >  2:2571193/1‑100 (MQ=255)
                                                                                        cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCg       >  1:3931401/1‑100 (MQ=255)
                                                                                           acATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTan    <  1:4431281/99‑2 (MQ=255)
                                                                                             agtAACGGGCAAATTTTCTGCTGAGAGGGCGGTGGAATTGTTGCCGATAGCCAGGAAGAAGCGGAATATCCGGAAACTTTTGATAAAGTTAATCGTATcg  <  1:130134/98‑2 (MQ=255)
                                                                                             aTTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATcc  >  1:2933810/1‑100 (MQ=255)
                                                                                             aTTCACGGACAAATTTTCTGCTCTGCGGGCGGAGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATcc  >  1:2933878/1‑100 (MQ=255)
                                                                                                |                                                                                                
CCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCC  >  NC_000913/1895954‑1896146

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: