Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,801,368 (T)7→6 coding (2000/2145 nt) nrdE → ribonucleoside‑diphosphate reductase 2 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,801,3620T.96.7% 526.7 / 9.4 122coding (1994/2145 nt)nrdEribonucleoside‑diphosphate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (48/70);  total (50/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 152 total aligned reads displayed.

GTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATC  >  NC_000913/2801270‑2801450
                                                                                            |                                                                                        
gTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCt                                                                                          >  1:4035623/1‑93 (MQ=255)
 tATCAGGACTCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                                                                         <  2:3190465/93‑1 (MQ=255)
 tATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                                                                         >  1:4260701/1‑93 (MQ=255)
 tATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                                                                         <  1:1768383/93‑1 (MQ=255)
  aTCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttt                                                                                        >  2:1929616/1‑93 (MQ=255)
  aTCAGGACGCTTACGAAATGGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttt                                                                                        <  2:1492368/93‑1 (MQ=255)
   tCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtttt                                                                                       <  1:3311847/93‑1 (MQ=255)
    cAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttt                                                                                      <  1:624629/93‑1 (MQ=255)
    cAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttt                                                                                      <  2:1880978/93‑1 (MQ=255)
     aGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtttttt                                                                                     >  1:1016904/1‑93 (MQ=255)
      ggACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttc                                                                                    >  2:1204646/1‑92 (MQ=255)
      ggACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttc                                                                                    <  1:3789477/93‑2 (MQ=255)
      ggACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttc                                                                                    >  2:1855168/1‑92 (MQ=255)
      ggACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttc                                                                                    >  1:3283175/1‑92 (MQ=255)
          gCTTACGAAATTGGCGCAGCAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTccccg                                                                                >  2:4570294/1‑92 (MQ=255)
           ctgaCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGa                                                                              <  1:206137/90‑1 (MQ=255)
           cTTACGAAATTGGCGCAGAAAAGATCATTGACACCTACGCGGAAGCGACTCGCCATGCCGATCAGGGGCTGTCGCTGGCGCTTTTTTCCCCGa                                                                              >  1:4359835/1‑93 (MQ=255)
           cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGa                                                                              >  2:4247093/1‑93 (MQ=255)
           cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGa                                                                              <  2:184422/93‑1 (MQ=255)
            ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGAt                                                                             >  1:2978671/1‑93 (MQ=255)
              aCGAAATTGGCGCAGAAAAGATCATCGACCCCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATAc                                                                           <  1:214682/93‑1 (MQ=255)
              aCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATAc                                                                           <  1:2604563/93‑1 (MQ=255)
                  aaTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGcc                                                                       <  1:4366774/93‑1 (MQ=255)
                  aaTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTAGCTGACGCTTTTTTCCACGATACCGcc                                                                       >  2:519822/1‑93 (MQ=255)
                    ttGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGccac                                                                     <  1:1460354/93‑1 (MQ=255)
                        cgcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACt                                                                 <  2:4488158/93‑1 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                >  2:1376141/1‑93 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                >  1:284381/1‑93 (MQ=255)
                           aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTcgc                                                              >  2:2011364/1‑93 (MQ=255)
                            gAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTcgcg                                                             >  2:2773763/1‑93 (MQ=255)
                            gAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTcgcg                                                             >  2:3561946/1‑93 (MQ=255)
                             aaaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGa                                                            <  2:453217/93‑1 (MQ=255)
                               aaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGata                                                          <  1:4314372/93‑1 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGatata                                                        >  2:1599788/1‑92 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                        >  1:4258396/1‑93 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                        <  1:1990766/93‑1 (MQ=255)
                                      tcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCTATACCGCCACCACTCGCGATATcaacaa                                                   >  2:1576115/1‑93 (MQ=255)
                                       cGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACaaa                                                  <  2:4656389/93‑1 (MQ=255)
                                       cGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACaaa                                                  <  2:1415190/93‑1 (MQ=255)
                                       cGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACaaa                                                  <  1:1306465/93‑1 (MQ=255)
                                        gACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAg                                                 >  1:3635843/1‑93 (MQ=255)
                                        gACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAg                                                 <  2:2501192/93‑1 (MQ=255)
                                        gACACCGACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAg                                                 <  1:1826110/93‑1 (MQ=255)
                                         acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgc                                                <  2:2834689/93‑1 (MQ=255)
                                          cacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAcgcg                                               >  2:1414861/1‑93 (MQ=255)
                                           acCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcgc                                              >  1:4265843/1‑93 (MQ=255)
                                           acCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcgc                                              <  1:1710844/93‑1 (MQ=255)
                                            ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCa                                             <  1:3488694/93‑1 (MQ=255)
                                            ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCa                                             >  1:1120391/1‑93 (MQ=255)
                                             cTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAg                                            >  1:559458/1‑93 (MQ=255)
                                              tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGa                                           <  2:1023976/93‑1 (MQ=255)
                                              tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGa                                           >  1:2024581/1‑93 (MQ=255)
                                              tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGa                                           <  1:4615224/93‑1 (MQ=255)
                                              tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGa                                           >  1:4505814/1‑93 (MQ=255)
                                              caCGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGa                                           <  1:137015/92‑1 (MQ=255)
                                               aCGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAt                                          <  1:194919/93‑1 (MQ=255)
                                                cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAtt                                         <  2:3983627/93‑1 (MQ=255)
                                                cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAtt                                         <  2:455280/93‑1 (MQ=255)
                                                cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAtt                                         >  1:3883166/1‑93 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        <  2:2944073/93‑1 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        >  2:4163921/1‑93 (MQ=255)
                                                   ggAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCTCTCGCGATATCAACAAAGCGCAGATTTAc                                      <  1:4126510/93‑1 (MQ=255)
                                                   ggAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTAc                                      <  2:1716707/93‑1 (MQ=255)
                                                      aGCGCCCCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                   <  2:2360078/93‑1 (MQ=255)
                                                      aGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                   <  1:1679695/93‑1 (MQ=255)
                                                       gCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCt                                  <  1:945002/93‑1 (MQ=255)
                                                        cGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTg                                 <  1:2983273/93‑1 (MQ=255)
                                                           cTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcg                              >  2:3043178/1‑93 (MQ=255)
                                                           cTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcg                              >  2:3629697/1‑93 (MQ=255)
                                                             cGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCa                            <  1:1240937/93‑1 (MQ=255)
                                                              gCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCaa                           <  1:4541780/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGTCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  2:807691/93‑1 (MQ=255)
                                                                cATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAgg                         <  2:2319432/93‑1 (MQ=255)
                                                                cATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAgg                         <  2:3122411/93‑1 (MQ=255)
                                                                 ctctCGATCAGGGGCTGTCGCTGCCGCTTTTTTCCCCGATTCCGCCACCACTCACGATATAAACAAAGCGCAGATTTACTCCTGACGCAAggg                        <  2:609283/90‑1 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        >  2:2700662/1‑93 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        <  2:2557745/93‑1 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        >  1:3547580/1‑93 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        >  2:3019905/1‑93 (MQ=255)
                                                                  tGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                       >  2:4440115/1‑93 (MQ=255)
                                                                  tGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                       >  2:2878002/1‑93 (MQ=255)
                                                                      gATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCa                   >  2:204225/1‑93 (MQ=255)
                                                                         cAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaaa                <  1:3704253/93‑1 (MQ=255)
                                                                          aGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAAc               <  2:2561789/93‑1 (MQ=255)
                                                                           ggggCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg              >  2:136786/1‑93 (MQ=255)
                                                                            gggCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGc             <  2:3273464/93‑1 (MQ=255)
                                                                               cTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct          <  2:2489014/93‑1 (MQ=255)
                                                                               cTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct          >  1:3284306/1‑93 (MQ=255)
                                                                               cTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct          <  1:4670664/93‑1 (MQ=255)
                                                                                tGTCGCTCACGCTTTTTTCCCCGAGACCGCCACCGCTCGCGATATCAACAAAGCGCAGATTTACGCGTGGCGCAAGGGTATCAAAACGCTCTa         <  2:4309934/93‑1 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        <  1:150176/93‑1 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        <  2:679548/93‑1 (MQ=255)
                                                                                  tCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt       >  1:1999360/1‑93 (MQ=255)
                                                                                   cGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTa      <  1:432297/93‑1 (MQ=255)
                                                                                   cGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTa      >  2:693934/1‑93 (MQ=255)
                                                                                    tcTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTAc     <  2:4204202/92‑1 (MQ=255)
                                                                                    gCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTAc     <  1:3506880/93‑1 (MQ=255)
                                                                                     cTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACa    >  1:2543384/1‑93 (MQ=255)
                                                                                       gACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc  >  1:3706557/1‑93 (MQ=255)
                                                                                       gACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc  <  1:4027392/93‑1 (MQ=255)
                                                                                            |                                                                                        
GTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATC  >  NC_000913/2801270‑2801450

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