Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,536,978 T→C T1213A (ACC→GCC)  narZ ← cryptic nitrate reductase 2 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,9780TC85.5% 597.1 / 79.4 234T1213A (ACC→GCC) narZcryptic nitrate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (24/10);  new base C (111/89);  total (135/99)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 246 total aligned reads displayed.

CTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAAT  >  NC_000913/1536886‑1537068
                                                                                            |                                                                                          
ctcctGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGt                                                                                            <  1:2695266/93‑1 (MQ=255)
  cctGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGgcgc                                                                                          >  2:475383/1‑93 (MQ=255)
  cctGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGgcgc                                                                                          >  1:3498842/1‑93 (MQ=255)
     gTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgt                                                                                       >  2:3233605/1‑93 (MQ=255)
     gTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgt                                                                                       <  1:6551765/93‑1 (MQ=255)
     gTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgt                                                                                       >  1:3851105/1‑93 (MQ=255)
      tACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgta                                                                                      <  1:3360810/93‑1 (MQ=255)
       aCCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCgtag                                                                                     >  2:6544453/1‑93 (MQ=255)
       aCCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgtag                                                                                     <  1:5644190/93‑1 (MQ=255)
       aCCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgtag                                                                                     >  1:2301251/1‑93 (MQ=255)
        ccTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCgtagt                                                                                    >  2:680620/1‑93 (MQ=255)
           gATCGCGACCTCCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGACCTCATCGCGGTTCGATCCGACGGTGCCGTagtagt                                                                                 <  1:5723056/93‑1 (MQ=255)
             tCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTa                                                                               >  1:423726/1‑93 (MQ=255)
               gcgACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTaaa                                                                             <  2:815524/93‑1 (MQ=255)
               gcgACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTaaa                                                                             <  1:5180416/93‑1 (MQ=255)
                cgACCTTCATCACCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACGCATCGCGATTCGATCCGACGGCGCCGTAGTAGTTAAAg                                                                            <  1:606991/93‑1 (MQ=255)
                  aCCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGcc                                                                          >  1:2762580/1‑93 (MQ=255)
                  aCCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGcc                                                                          >  2:2117195/1‑93 (MQ=255)
                   ccTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGccc                                                                         >  1:1370382/1‑93 (MQ=255)
                      tcatcaTCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAg                                                                      >  2:6357067/1‑93 (MQ=255)
                      tcatcaTCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAAATCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAg                                                                      >  1:1924100/1‑93 (MQ=255)
                       catcatCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAgg                                                                     >  2:3068549/1‑93 (MQ=255)
                        atcatcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGc                                                                    >  2:1809331/1‑93 (MQ=255)
                        atcatcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGc                                                                    <  2:3730812/93‑1 (MQ=255)
                        atcatcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGc                                                                    <  2:2082442/93‑1 (MQ=255)
                         tcatcCAGTCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGcc                                                                   >  1:6478159/1‑93 (MQ=255)
                         tcatcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATTAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGcc                                                                   >  1:6377597/1‑93 (MQ=255)
                          catcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCa                                                                  <  1:3950221/93‑1 (MQ=255)
                           atcCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAg                                                                 >  2:2639717/1‑93 (MQ=255)
                            tcCAGCCAGTTAACGTTCTTCATCATGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGc                                                                <  1:1759624/93‑1 (MQ=255)
                             ccagccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCt                                                               >  1:2277369/1‑93 (MQ=255)
                              cagccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTg                                                              <  2:3364692/93‑1 (MQ=255)
                                gccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTgcg                                                            >  2:5182352/1‑93 (MQ=255)
                                gccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTgcg                                                            <  1:6089000/93‑1 (MQ=255)
                                gccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTgcg                                                            >  1:1255307/1‑93 (MQ=255)
                                gccagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGACCCAGGCCAGCTgcg                                                            >  2:2329997/1‑93 (MQ=255)
                                  cagTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGcgcg                                                          >  1:4868553/1‑93 (MQ=255)
                                   agTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGt                                                         >  1:4835427/1‑93 (MQ=255)
                                   agTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGt                                                         <  2:2728116/93‑1 (MQ=255)
                                     ttAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTaa                                                       <  2:2498072/93‑1 (MQ=255)
                                      tAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAAc                                                      <  1:2796338/93‑1 (MQ=255)
                                        aCGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCg                                                    >  1:3647210/1‑93 (MQ=255)
                                          gTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGcc                                                  >  2:5312117/1‑93 (MQ=255)
                                           ttcttcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCa                                                 >  1:5550971/1‑93 (MQ=255)
                                            tcttcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCaa                                                >  2:6479346/1‑93 (MQ=255)
                                              ttcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATc                                              >  1:1529602/1‑93 (MQ=255)
                                               tcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCa                                             >  1:3248348/1‑93 (MQ=255)
                                               tcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCa                                             <  2:5267526/93‑1 (MQ=255)
                                               tcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCa                                             >  2:3929382/1‑93 (MQ=255)
                                               tcATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCa                                             >  1:5300164/1‑93 (MQ=255)
                                                 aTCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCata                                           <  1:1499721/93‑1 (MQ=255)
                                                  tCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCatat                                          <  2:6270264/93‑1 (MQ=255)
                                                  tCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCatat                                          <  1:6179155/93‑1 (MQ=255)
                                                  tCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCatat                                          >  2:1351767/1‑93 (MQ=255)
                                                  tCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCatat                                          <  1:4081254/93‑1 (MQ=255)
                                                   cTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATg                                         <  1:3404475/93‑1 (MQ=255)
                                                    ttGCGGATCATGATGCACTCATCGCGCTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGc                                        <  2:1579382/93‑1 (MQ=255)
                                                    ttGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAACCATAAGc                                        >  1:3741375/1‑93 (MQ=255)
                                                      gCGGATCATGATGAACTCATCGCGGTTCGTTCCGACGGCTCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGt                                      <  1:4444515/93‑1 (MQ=255)
                                                      gCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGt                                      <  1:1606746/93‑1 (MQ=255)
                                                       cGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGtt                                     <  2:5065158/93‑1 (MQ=255)
                                                       cGGATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGtt                                     >  2:4451695/1‑93 (MQ=255)
                                                        ggATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTg                                    <  1:505200/93‑1 (MQ=255)
                                                         gATCATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTgg                                   >  1:979666/1‑93 (MQ=255)
                                                            cATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGttt                                >  1:4351544/1‑93 (MQ=255)
                                                            cATGATGAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGttt                                <  2:5804700/93‑1 (MQ=255)
                                                             atgatgAACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTc                               <  2:5626751/93‑1 (MQ=255)
                                                              tgatgaACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCg                              <  2:4227447/93‑1 (MQ=255)
                                                               gatgaACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCgg                             >  2:2264243/1‑93 (MQ=255)
                                                                 tgaACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGc                           >  1:834552/1‑93 (MQ=255)
                                                                 tgaACTCATCGCGGTTCGATCCGACGGCGCCTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCa                          <  2:3426240/93‑1 (MQ=255)
                                                                  gaACTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCa                          <  2:282816/93‑1 (MQ=255)
                                                                   aaCTCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCaa                         >  2:3627843/1‑93 (MQ=255)
                                                                      tCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAcg                      >  2:4393222/1‑93 (MQ=255)
                                                                      tCATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAcg                      >  2:4601054/1‑93 (MQ=255)
                                                                       cATCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAcgc                     <  2:6574735/93‑1 (MQ=255)
                                                                        aTCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAcgcg                    >  2:2939609/1‑93 (MQ=255)
                                                                         tCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCgg                   <  2:4656510/93‑1 (MQ=255)
                                                                         tCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCgg                   >  1:1014398/1‑93 (MQ=255)
                                                                         tCGCGGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCgg                   <  2:57143/93‑1 (MQ=255)
                                                                           gcgGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGt                 <  2:6071734/93‑1 (MQ=255)
                                                                            cgGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTa                >  1:2883960/1‑93 (MQ=255)
                                                                            cgGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTa                >  1:588282/1‑93 (MQ=255)
                                                                            cgGTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTa                <  2:1832017/93‑1 (MQ=255)
                                                                             ggTTCGATCCGACGGTGCCGTGGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTaa               >  2:5426393/1‑93 (MQ=255)
                                                                             ggTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTaa               <  1:2064829/93‑1 (MQ=255)
                                                                              gTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAc              >  1:1838052/1‑93 (MQ=255)
                                                                              gTTCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAc              >  2:4875601/1‑93 (MQ=255)
                                                                               ttCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAcc             <  1:4783889/93‑1 (MQ=255)
                                                                               ttCGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAcc             <  1:1136676/93‑1 (MQ=255)
                                                                                 cGATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGa           <  2:6408840/93‑1 (MQ=255)
                                                                                  gATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAg          >  1:6287527/1‑93 (MQ=255)
                                                                                  gATCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAg          >  1:4612637/1‑93 (MQ=255)
                                                                                   aTCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGt         >  1:5470640/1‑93 (MQ=255)
                                                                                    tCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGtt        <  2:5408742/93‑1 (MQ=255)
                                                                                    tCCGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGtt        >  2:5873889/1‑93 (MQ=255)
                                                                                      cGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTAt      <  1:2452413/93‑1 (MQ=255)
                                                                                      cGACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTAt      >  1:5508067/1‑93 (MQ=255)
                                                                                       gACGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATg     <  1:360785/93‑1 (MQ=255)
                                                                                         cGGCGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGaa   >  2:3836536/1‑93 (MQ=255)
                                                                                            |                                                                                          
CTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAAT  >  NC_000913/1536886‑1537068

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