Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,068 +C intergenic (‑178/‑242) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,0681.C99.5% 621.2 / ‑4.0 192intergenic (‑178/‑242)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (44/147);  total (44/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 203 total aligned reads displayed.

GGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTT  >  NC_000913/4476975‑4477154
                                                                                               |                                                                                      
ggTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacg                                                                                         <  1:4250748/93‑1 (MQ=255)
 gtgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                                                                        >  1:3776881/1‑92 (MQ=255)
  tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgca                                                                                       >  1:152797/1‑91 (MQ=255)
    tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgcatt                                                                                     >  2:2900490/1‑89 (MQ=255)
      gCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTAc                                                                                   >  1:1580455/1‑93 (MQ=255)
       ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                  >  2:3239100/1‑93 (MQ=255)
       ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                  <  1:5114762/93‑1 (MQ=255)
       ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                  <  2:2575818/93‑1 (MQ=255)
       ccTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                                                                  >  2:2792370/1‑93 (MQ=255)
        cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                                 >  2:2432201/1‑93 (MQ=255)
        cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                                 <  2:4157878/93‑1 (MQ=255)
        cTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                                                                 >  2:188473/1‑93 (MQ=255)
         tttATCAGCGGTGAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAgg                                                                                <  2:5301336/93‑1 (MQ=255)
         tttATCAGAGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGTATTATCCTGCAGCGTCAGGCGACGCATTACAgg                                                                                <  1:413970/93‑1 (MQ=255)
          ttATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                                                               <  2:3089735/93‑1 (MQ=255)
            aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCTGGCGACGCATTACAGGTGa                                                                             <  2:6159145/93‑1 (MQ=255)
            aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                                             >  2:6292715/1‑93 (MQ=255)
            aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                                             <  1:31865/93‑1 (MQ=255)
            aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                                                             <  1:703647/93‑1 (MQ=255)
             tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAt                                                                            <  2:4436400/93‑1 (MQ=255)
              cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                                                                           <  1:2059320/93‑1 (MQ=255)
                gCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                                                         <  2:2986572/93‑1 (MQ=255)
                gCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                                                                         <  2:3054321/93‑1 (MQ=255)
                 cGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcg                                                                        <  2:3731005/93‑1 (MQ=255)
                  ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                                                                       <  2:5757658/93‑1 (MQ=255)
                   gTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggg                                                                      >  2:5606024/1‑92 (MQ=255)
                       gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                                  >  1:2157042/1‑93 (MQ=255)
                       gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                                  <  2:899533/93‑1 (MQ=255)
                       gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                                  >  2:1310206/1‑93 (MQ=255)
                       gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                                                                  <  2:5951397/93‑1 (MQ=255)
                       aaCCGTATTCGGCGGATACCTGACGAATCACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGACTGCGGCGCaa                                                                  <  2:1118274/92‑1 (MQ=255)
                        aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAt                                                                 <  1:5047627/93‑1 (MQ=255)
                         ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                                                                <  1:5493575/93‑1 (MQ=255)
                          cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                                                               >  1:1717395/1‑93 (MQ=255)
                          cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                                                               >  1:2104947/1‑93 (MQ=255)
                          cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                                                               >  1:5750268/1‑93 (MQ=255)
                            tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                                             <  2:2144539/93‑1 (MQ=255)
                            tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                                             <  2:608636/93‑1 (MQ=255)
                            tATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTa                                                             <  2:1624233/93‑1 (MQ=255)
                             aTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAt                                                            <  1:1907114/93‑1 (MQ=255)
                              ttCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt                                                           <  1:4762182/93‑1 (MQ=255)
                                cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                         <  1:3483369/93‑1 (MQ=255)
                                cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa                                                         <  1:6009643/93‑1 (MQ=255)
                                 ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                                                        <  2:4322078/93‑1 (MQ=255)
                                 ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt                                                        <  1:3081613/93‑1 (MQ=255)
                                  gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa                                                       <  1:4413784/93‑1 (MQ=255)
                                  gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATa                                                       >  1:2571866/1‑93 (MQ=255)
                                    ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt                                                     <  1:3722131/93‑1 (MQ=255)
                                       tACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                                  >  1:3897640/1‑93 (MQ=255)
                                       tACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa                                                  <  1:4966460/93‑1 (MQ=255)
                                        aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                                                 <  2:3925070/93‑1 (MQ=255)
                                        aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa                                                 <  1:4325991/93‑1 (MQ=255)
                                          cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc                                               >  1:465465/1‑93 (MQ=255)
                                           tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                                              <  1:5697839/93‑1 (MQ=255)
                                           tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc                                              >  2:1721213/1‑93 (MQ=255)
                                            gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCa                                             >  2:4494519/1‑93 (MQ=255)
                                             aCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaa                                            <  2:5160311/93‑1 (MQ=255)
                                              cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaat                                           <  2:2202394/93‑1 (MQ=255)
                                               gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaata                                          <  1:3269176/93‑1 (MQ=255)
                                                aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa                                         >  1:482732/1‑93 (MQ=255)
                                                aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataa                                         >  1:3950003/1‑93 (MQ=255)
                                                 aTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaa                                        <  1:5362864/93‑1 (MQ=255)
                                                  tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaa                                       <  2:4078352/93‑1 (MQ=255)
                                                   gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaa                                      <  1:6254042/93‑1 (MQ=255)
                                                     cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaata                                    <  2:5633335/93‑1 (MQ=255)
                                                       cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaataaa                                  <  2:4831810/93‑1 (MQ=255)
                                                        gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAg                                 >  2:1857359/1‑93 (MQ=255)
                                                        gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAg                                 >  1:1983465/1‑93 (MQ=255)
                                                        gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAg                                 <  2:4528107/93‑1 (MQ=255)
                                                         ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGt                                <  2:795837/93‑1 (MQ=255)
                                                          gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACTCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa                               <  2:4702995/93‑1 (MQ=255)
                                                          gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa                               <  1:5590860/93‑1 (MQ=255)
                                                            gATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAt                             <  1:6238869/93‑1 (MQ=255)
                                                               tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                           <  2:5064800/89‑1 (MQ=255)
                                                               tagcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                           <  1:973781/89‑1 (MQ=255)
                                                                agcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtat                          <  1:1173210/90‑1 (MQ=255)
                                                                 gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata                         >  1:263734/3‑93 (MQ=255)
                                                                 gcaGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATtata                         <  1:2653444/91‑1 (MQ=255)
                                                                  caGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATAc                        >  2:2246161/2‑93 (MQ=255)
                                                                  caGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATAc                        <  1:5248079/92‑1 (MQ=255)
                                                                   aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa                       <  1:3369392/93‑1 (MQ=255)
                                                                   aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa                       <  2:2068974/93‑1 (MQ=255)
                                                                   aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACa                       <  2:4600850/93‑1 (MQ=255)
                                                                    ggATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAt                      <  1:828594/93‑1 (MQ=255)
                                                                     gATTATCGTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAtt                     <  1:3028591/93‑1 (MQ=255)
                                                                     gATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAtt                     <  2:3696536/93‑1 (MQ=255)
                                                                     gATAATCCGGCAGCTTCAGGCGACGCATAACAGGTGATTGCGGCGCAAGGTAATTCATATTTAAACCCAAAAAAAATAAAGTAATTATACAtt                     <  2:3945661/93‑1 (MQ=255)
                                                                        tATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGt                  <  2:1535616/93‑1 (MQ=255)
                                                                         aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGtt                 >  1:534630/1‑93 (MQ=255)
                                                                           ccTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTa               <  2:4126394/93‑1 (MQ=255)
                                                                            cTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTaa              >  2:1935712/1‑93 (MQ=255)
                                                                             tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAAt             <  2:3491053/93‑1 (MQ=255)
                                                                               cAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATAc           <  1:6290870/93‑1 (MQ=255)
                                                                                  cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc        >  2:4408169/1‑93 (MQ=255)
                                                                                  cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc        <  1:3869159/93‑1 (MQ=255)
                                                                                  cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc        <  1:2267678/93‑1 (MQ=255)
                                                                                   gTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAct       >  2:2123231/1‑93 (MQ=255)
                                                                                    tCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCActc      <  1:1442650/93‑1 (MQ=255)
                                                                                       ggCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCActctct   <  1:106575/93‑1 (MQ=255)
                                                                                       ggCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCActctct   <  1:3631949/93‑1 (MQ=255)
                                                                                               |                                                                                      
GGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTT  >  NC_000913/4476975‑4477154

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