Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,142,505 G→C 20.1% R148R (CGC→CGG)  udk ← uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,142,5050GC20.1% 339.9 / 83.8 145R148R (CGC→CGG) udkuridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base G (56/59);  new base C (11/18);  total (67/78)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 157 total aligned reads displayed.

TCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATG  >  NC_000913/2142407‑2142604
                                                                                                  |                                                                                                   
tCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                    <  2:17970/100‑1 (MQ=255)
tCAATGAATTGCAGGAACATCGGGAGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggc                                                                                                    <  2:1472231/100‑1 (MQ=255)
 cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggcg                                                                                                   <  2:64328/100‑1 (MQ=255)
 cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggcg                                                                                                   <  1:236781/100‑1 (MQ=255)
 cAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggcg                                                                                                   <  1:3679377/100‑1 (MQ=255)
  aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATgcggcg                                                                                                   <  1:276676/99‑1 (MQ=255)
  aaTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGgcgc                                                                                                  <  1:1992879/100‑1 (MQ=255)
   aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCa                                                                                                 <  1:1431206/100‑1 (MQ=255)
   aTGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCa                                                                                                 >  1:1839244/1‑100 (MQ=255)
     gAATTGCAGGAACATCGGGCGCACGGATTTTTGATATTGTGCCATCACTGAATCCATTGCACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATc                                                                                               >  2:4013072/1‑100 (MQ=255)
      aaTTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCa                                                                                              <  1:1526569/100‑1 (MQ=255)
        ttGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAgg                                                                                            <  2:3329229/100‑1 (MQ=255)
         tGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                           <  1:527602/100‑1 (MQ=255)
         tGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGc                                                                                           <  1:2498471/100‑1 (MQ=255)
         tGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGc                                                                                           <  1:1707433/100‑1 (MQ=255)
          gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCa                                                                                          <  1:907330/100‑1 (MQ=255)
          gCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGcg                                                                                          <  1:3786654/100‑2 (MQ=255)
           cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGccg                                                                                         >  1:1108993/1‑98 (MQ=255)
           cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAg                                                                                         <  1:1038411/100‑1 (MQ=255)
           cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAg                                                                                         >  1:2388453/1‑100 (MQ=255)
           cAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGcgt                                                                                         <  2:3391654/100‑3 (MQ=255)
            aGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaagcgtg                                                                                        <  1:543327/100‑7 (MQ=255)
            aGGAACATCGGGCGCACGGTTTTTTGATATCGCGCCATCACTGAATCCATTGAACGCCCACGCGCGTTAACGTCACGCTTGATGCGCCGCATCAAGcg                                                                                          <  2:3028486/98‑2 (MQ=255)
             ggAACATCGGGCGCACGGTTTTTTTTTTTTTTTTCATCACTGAAGCCATTGAACGCACACTCTCGTTTACGGCACGCTTGATGCGGCGCATCACGCata                                                                                        >  2:3641742/1‑97 (MQ=255)
             ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGat                                                                                       >  1:2790018/1‑100 (MQ=255)
             ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGat                                                                                       >  2:4086895/1‑100 (MQ=255)
             ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaagcgtga                                                                                       <  1:1800164/100‑8 (MQ=255)
             ggAACATCGGGCGCACGGTTTTTTGATATTGCGCCAGCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCaagcgtga                                                                                       <  1:2940392/100‑8 (MQ=255)
             ggAACATCGAGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATCGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGCCGCATCAAGcgtg                                                                                        <  2:2639355/99‑4 (MQ=255)
              tccgatctggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGata                                                                                      <  2:3132905/92‑1 (MQ=255)
                aCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCGGATATc                                                                                    <  1:479222/100‑1 (MQ=255)
                aCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATc                                                                                    >  1:3225543/1‑100 (MQ=255)
                aCATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATc                                                                                    <  1:492308/100‑1 (MQ=255)
                 cATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGcgtgacgtt                                                                                   <  2:3222181/100‑9 (MQ=255)
                  aTCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCa                                                                                  >  2:481437/1‑100 (MQ=255)
                  aTCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCa                                                                                  <  1:804490/100‑1 (MQ=255)
                      ggCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGataagatcgga                                                                              >  1:3132919/1‑92 (MQ=255)
                        cgcACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgggg                                                                            <  1:889903/100‑1 (MQ=255)
                        cgcACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgggg                                                                            <  1:2629141/100‑1 (MQ=255)
                        cgcACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCgggg                                                                            <  1:3371556/100‑1 (MQ=255)
                         gcACGGTTTTTTGATATTGCGCCATCTCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGt                                                                           <  1:3352162/100‑1 (MQ=255)
                          cACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTa                                                                          >  1:3859669/1‑100 (MQ=255)
                           aCGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTAt                                                                         <  1:1053270/100‑1 (MQ=255)
                           aCGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTAt                                                                         >  1:2568245/1‑100 (MQ=255)
                            cGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATc                                                                        <  2:478227/100‑1 (MQ=255)
                                ttcttGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACg                                                                    >  1:575362/4‑100 (MQ=255)
                                    tGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGaaaa                                                                >  1:3412921/1‑100 (MQ=255)
                                     gATCTTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATGTCCAGCGGGGTATCAACGa                                                                   >  1:452773/1‑96 (MQ=255)
                                       tatTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATgg                                                             >  1:1318377/1‑100 (MQ=255)
                                        atTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                            >  2:2089908/1‑100 (MQ=255)
                                        atTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGga                                                            <  1:297432/100‑1 (MQ=255)
                                         ttGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgag                                                           <  1:3701419/100‑1 (MQ=255)
                                          tGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGgaga                                                          <  1:1207054/100‑1 (MQ=255)
                                           gcgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCCGGCAGATATCCAGCGGGGTATCAACGAAAATGgaga                                                          >  1:625889/1‑99 (MQ=255)
                                            cgcCATCCCCGAATCCATTGAACGCCCACGCTCGTTAACCTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaa                                                         <  2:2956889/99‑1 (MQ=255)
                                            cgcCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaag                                                        <  2:3393417/100‑1 (MQ=255)
                                               cATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagtta                                                     <  2:3859644/100‑1 (MQ=255)
                                                aTCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                    <  2:36725/100‑1 (MQ=255)
                                                aTCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaa                                                    >  2:3057787/1‑100 (MQ=255)
                                                  cACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGAGGTATCAACGAAAATGGAGaagttaagt                                                  >  1:2297521/1‑100 (MQ=255)
                                                   aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagtt                                                 >  2:84619/1‑100 (MQ=255)
                                                   aCTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGaagttaagtt                                                 <  1:2592918/100‑1 (MQ=255)
                                                     tGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCg                                               <  1:4106780/100‑1 (MQ=255)
                                                     tGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCg                                               <  1:4107292/100‑1 (MQ=255)
                                                      gAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGt                                              >  2:1065945/1‑100 (MQ=255)
                                                         tCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATAAACGAAAATGGAGAAGTTAAGTTCGTCAc                                           >  2:1401358/1‑100 (MQ=255)
                                                           cATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGc                                         >  2:1261286/1‑100 (MQ=255)
                                                               gAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAc                                     >  1:3525686/1‑100 (MQ=255)
                                                                    cccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCa                                >  1:3002965/1‑100 (MQ=255)
                                                                    cccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCa                                <  1:1134747/100‑1 (MQ=255)
                                                                     ccACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCAt                               >  1:1105168/1‑100 (MQ=255)
                                                                     ccACGCTCGTTAACGTCACGCTTGATGCGGCCCATAAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAAcgcg                                  <  2:3262139/97‑1 (MQ=255)
                                                                          cTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTc                          <  2:211374/100‑1 (MQ=255)
                                                                            cGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAg                        <  1:1117076/100‑1 (MQ=255)
                                                                            cGTTAACGTCACGCTTCATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAg                        >  2:2385258/1‑100 (MQ=255)
                                                                             gTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGc                       >  1:4173353/1‑100 (MQ=255)
                                                                                 aCGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGcaaca                   >  2:4041144/1‑100 (MQ=255)
                                                                                  cGTCACGCTTGATGCGGCGCATCAGGCCGCTATCCAGCGGGGTCTCACCGAAAACGGAGAAGTTCAGTTCGTCACGCAAACGCGCATCCGTCAGcaacac                  <  2:2500740/100‑2 (MQ=255)
                                                                                   gTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGGCACGCAAACGCGCATCCGTCAGCAACaaa                 >  2:3556429/1‑100 (MQ=255)
                                                                                   gTCACGCTTGATGCGGCGCAGCAGGCAGATATCCAGCGGGTGATCAACGAAAATGGAGATGTTAGGTGCGTCACGCAAGCGCGCAGCCg                            >  2:3555687/1‑89 (MQ=255)
                                                                                      aCGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATg              >  1:801047/1‑100 (MQ=255)
                                                                                       tgcctGATGCGCCGCATCAAGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGc             >  2:165592/5‑100 (MQ=255)
                                                                                       cGCTTGATGCGGCGCATCAGACAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGc             >  2:2556113/1‑100 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAGCGAAACTGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAACTGc             >  2:3357398/4‑99 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGcc            >  2:2360659/4‑100 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGcc            >  2:1673171/4‑100 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGCAGATATCCAGCGAGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCGACAAAATGcc            >  2:2748637/4‑100 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGCAGATATCCAGCGAGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCGACAAAATGcc            >  2:2748743/4‑100 (MQ=255)
                                                                                        gcctGATGCGCCGCATCAGGAAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGcc            >  2:2561206/4‑100 (MQ=255)
                                                                                          ttGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCtt          <  2:3908289/100‑1 (MQ=255)
                                                                                           tGATGCGGCGCATCAGGCAGCTATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTc         >  2:2561696/1‑100 (MQ=255)
                                                                                            gATGCGCCGCATCAGGCAGATATCCAGCGGGGTATTAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCATTCg        >  2:2932315/1‑100 (MQ=255)
                                                                                             aTGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCga       >  2:3483272/1‑100 (MQ=255)
                                                                                              tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCgag      <  1:2530286/100‑1 (MQ=255)
                                                                                              tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCgag      >  1:2851258/1‑100 (MQ=255)
                                                                                              tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCgag      >  1:2026082/1‑100 (MQ=255)
                                                                                              tGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCgag      <  1:3192033/100‑1 (MQ=255)
                                                                                               gcggcgCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCgaga     >  2:2298316/1‑100 (MQ=255)
                                                                                                 ggcgCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAat   <  1:3649589/100‑1 (MQ=255)
                                                                                                  gcgcATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAatg  <  2:801047/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCAATGAATTGCAGGAACATCGGGCGCACGGTTTTTTGATATTGCGCCATCACTGAATCCATTGAACGCCCACGCTCGTTAACGTCACGCTTGATGCGGCGCATCAGGCAGATATCCAGCGGGGTATCAACGAAAATGGAGAAGTTAAGTTCGTCACGCAAACGCGCATCCGTCAGCAACAAAATGCCTTCGAGAATG  >  NC_000913/2142407‑2142604

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: ‑