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,0500AC99.3% 489.5 / ‑4.2 144N416H (AAC→CAC) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (66/76);  total (67/77)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 144 total aligned reads displayed.

ACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCC  >  NC_000913/1895953‑1896146
                                                                                                 |                                                                                                
  cGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCTGCACGCTGACTGCCATTCACgg                                                                                              <  1:3412429/100‑1 (MQ=255)
  cGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACgg                                                                                              >  1:1587530/1‑100 (MQ=255)
    cAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGAc                                                                                            <  2:1404770/100‑1 (MQ=255)
    cAGCGACGCAATGCAAGGTGCGGAAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGAc                                                                                            <  2:188755/100‑1 (MQ=255)
     aGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACa                                                                                           >  1:3063156/1‑100 (MQ=255)
        gACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAAt                                                                                        >  1:1987972/1‑100 (MQ=255)
         aCGCACTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAAtt                                                                                       >  2:1603077/1‑100 (MQ=255)
            cAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtc                                                                                    <  2:506232/100‑1 (MQ=255)
            cAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtc                                                                                    >  1:2570133/1‑100 (MQ=255)
            cAATGCCTGGTGCGACAGCATTGGCTGTTTGAGCTTTTGCGGCAACAGGGATACGAGTATTACTATCCGCACGCTGACTGCCCTTCACGGACAAATTTtc                                                                                    >  2:3738087/1‑100 (MQ=255)
             aaTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtct                                                                                   <  1:2882080/100‑1 (MQ=255)
             aaTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtct                                                                                   <  2:2144012/100‑1 (MQ=255)
              aTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctg                                                                                  >  1:1719802/1‑100 (MQ=255)
              aTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctg                                                                                  >  2:124425/1‑100 (MQ=255)
              aTGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctg                                                                                  >  2:347427/1‑100 (MQ=255)
                 ccTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctgctc                                                                               >  2:2664341/1‑100 (MQ=255)
                   tGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTtctgctctg                                                                             >  2:2957209/1‑100 (MQ=255)
                     gTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcg                                                                           >  1:2442523/1‑100 (MQ=255)
                        cGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggc                                                                        <  1:352721/100‑1 (MQ=255)
                         ggCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggcg                                                                       >  1:2854379/1‑100 (MQ=255)
                         ggCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggcg                                                                       >  2:1812468/1‑100 (MQ=255)
                          gcggcaGTGGCTATTTGAGCGTGTGCGGTAAGATGGATACCAGTATTACTATCCGCCCGCTGACTGCGATTCACGGACAAATTTTCTGCTCTgcgggcgg                                                                      <  2:2468788/97‑1 (MQ=255)
                          gcagcaTTGGCTATTTGAGCTTTTGCGGCAACGTGGATACCAGTATTACTATCCGCAGGCTGACTGCCATTCATGGACAAATTTTCTGCTCTgcgggcgg                                                                      >  2:1444946/1‑100 (MQ=255)
                          gcagcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTgcgggcgg                                                                      >  1:49778/1‑100 (MQ=255)
                             gcaTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTgg                                                                   <  2:2774080/100‑1 (MQ=255)
                              caTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGa                                                                  <  2:3435370/100‑1 (MQ=255)
                              caTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGa                                                                  <  1:2492326/100‑1 (MQ=255)
                              caTTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGa                                                                  <  2:2534263/100‑1 (MQ=255)
                                ttGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAt                                                                <  2:1041295/100‑1 (MQ=255)
                                ttGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAt                                                                <  1:982290/100‑1 (MQ=255)
                                ctGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAAt                                                                <  2:1041372/99‑1 (MQ=255)
                                  ggCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTg                                                              >  1:722280/1‑100 (MQ=255)
                                    cTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTc                                                            <  1:2463807/100‑1 (MQ=255)
                                    cTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTc                                                            <  1:561821/100‑1 (MQ=255)
                                     tATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCg                                                           <  1:1284596/100‑1 (MQ=255)
                                      aTTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGc                                                          >  1:3031138/1‑100 (MQ=255)
                                        ttGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCg                                                        <  1:882034/100‑1 (MQ=255)
                                         tGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGa                                                       >  1:3383809/1‑100 (MQ=255)
                                          gAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGAt                                                      >  2:2587085/1‑100 (MQ=255)
                                            gCTTTTGCGGCAACATGGATACCAGTATTACTATCTGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAg                                                    <  2:3942048/100‑1 (MQ=255)
                                            gCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAg                                                    <  2:3566834/100‑1 (MQ=255)
                                             cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGc                                                   >  2:1735245/1‑100 (MQ=255)
                                              ttttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGcc                                                  >  1:1077576/1‑100 (MQ=255)
                                              ttttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGcc                                                  <  1:2283767/100‑1 (MQ=255)
                                              ttttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGcc                                                  >  1:1488744/1‑100 (MQ=255)
                                               tttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCa                                                 <  2:3308020/100‑1 (MQ=255)
                                                ttGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAg                                                <  2:2873049/100‑1 (MQ=255)
                                                  gCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGga                                              <  2:3194708/100‑1 (MQ=255)
                                                   cGGCGACGTGGGTACGAGGATTACTATCCTCACGCTGACTGCCATTCACGGATAAATTTGCGTCGCTGCGGGCGGTGGAATTGTCGCCGATAGTCAGtaa                                             >  2:2404670/1‑100 (MQ=255)
                                                   cGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaa                                             >  2:1629240/1‑100 (MQ=255)
                                                   cGGCAACATGGATACCAGTATTACCATCCCCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGTATTGTCGCCGATAGCCAGgaa                                             >  1:3805562/1‑100 (MQ=255)
                                                    ggCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGAGAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaag                                            >  2:3619615/1‑100 (MQ=255)
                                                    ggCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCAGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCAAGgaag                                            <  2:3069695/100‑1 (MQ=255)
                                                        aCATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGc                                        <  1:3281785/100‑1 (MQ=255)
                                                        aCATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGc                                        >  1:4069803/1‑100 (MQ=255)
                                                          aTGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCgg                                      <  2:509944/100‑1 (MQ=255)
                                                           tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     <  1:216707/100‑1 (MQ=255)
                                                           tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     >  1:3908355/1‑100 (MQ=255)
                                                           tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     <  2:323753/100‑1 (MQ=255)
                                                           tGGATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGa                                     >  1:2780396/1‑100 (MQ=255)
                                                            ggATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGaa                                    >  1:2168762/1‑100 (MQ=255)
                                                            ggATACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGaa                                    <  1:3228529/100‑1 (MQ=255)
                                                              ttACCAATATTAAGATGCGCACGCGGACTGACATTCACTGACAAATTGTCTGCTCCGGGGGCGCTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAata                                  <  2:3317711/99‑1 (MQ=255)
                                                              aTACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAata                                  <  1:428695/100‑1 (MQ=255)
                                                              aTACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAata                                  <  2:3106765/100‑1 (MQ=255)
                                                               tACCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAGGAAGCGGAatat                                 >  1:1203607/1‑100 (MQ=255)
                                                                aCCAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATc                                <  1:1745956/100‑1 (MQ=255)
                                                                 ccAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCa                               >  2:2888921/1‑100 (MQ=255)
                                                                 ccAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCa                               >  1:3673200/1‑100 (MQ=255)
                                                                 ccAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCa                               >  2:1020833/1‑100 (MQ=255)
                                                                  cAGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAg                              <  2:1340016/100‑1 (MQ=255)
                                                                  cAGTATTACTATCCGCACGCTGACTGCCACTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAg                              <  1:3726912/100‑1 (MQ=255)
                                                                   aGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAgg                             <  1:1656874/100‑1 (MQ=255)
                                                                   aGTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAgg                             <  1:2892394/100‑1 (MQ=255)
                                                                    gTATTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGa                            <  2:1517756/100‑1 (MQ=255)
                                                                      aTTACTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGaaa                          >  2:1765408/1‑100 (MQ=255)
                                                                         aCTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATCTCAGGAAACtt                       >  2:1093538/1‑100 (MQ=255)
                                                                         aCTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtt                       >  1:753037/1‑100 (MQ=255)
                                                                          cTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACttt                      <  1:624818/100‑1 (MQ=255)
                                                                          cTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACttt                      >  2:259709/1‑100 (MQ=255)
                                                                          cTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACttt                      >  1:3165055/1‑100 (MQ=255)
                                                                          cTATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACttt                      <  1:1042473/100‑1 (MQ=255)
                                                                           tATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtttt                     <  2:1959188/100‑1 (MQ=255)
                                                                           tATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtttt                     <  2:1922095/100‑1 (MQ=255)
                                                                           tATCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACtttt                     <  2:1170272/100‑1 (MQ=255)
                                                                            aTCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTg                    >  2:595885/1‑100 (MQ=255)
                                                                             tCCGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGa                   <  2:558395/100‑1 (MQ=255)
                                                                              ccGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                  <  1:2709811/100‑1 (MQ=255)
                                                                              ccGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                  <  1:3961810/100‑1 (MQ=255)
                                                                               cGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATa                 >  1:2201313/1‑100 (MQ=255)
                                                                               cGCACGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATa                 >  1:1368115/1‑100 (MQ=255)
                                                                                   cGCTGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGt             >  1:2523923/1‑100 (MQ=255)
                                                                                      tGACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTaa          >  2:463987/1‑100 (MQ=255)
                                                                                       gACTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAAt         <  2:1281013/100‑1 (MQ=255)
                                                                                        aCTGCCATTGACGGACAAGTTATCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCTGAATATCAGGAAGCTTTTGATAAAGTTAATc        <  2:1548263/100‑1 (MQ=255)
                                                                                         cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCg       >  1:1468109/1‑100 (MQ=255)
                                                                                         cTGCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCg       >  1:2100206/1‑100 (MQ=255)
                                                                                           gCCATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTa     >  1:3695355/1‑100 (MQ=255)
                                                                                            ccATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTAt    <  1:2705478/100‑1 (MQ=255)
                                                                                             cATTCACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATc   >  1:1658937/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGACAAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCC  >  NC_000913/1895953‑1896146

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