Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 814,470 A→G K574K (AAA→AAG uvrB → DNA repair; excision nuclease subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913814,4700AG97.4% 596.2 / 1.9 190K574K (AAA→AAGuvrBDNA repair; excision nuclease subunit B
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (98/87);  total (101/89)
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.46e-01

Only 100 of 198 total aligned reads displayed.

TAACGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAACAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACATTGC  >  NC_000913/814380‑814559
                                                                                          |                                                                                         
tAACGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGca                                                                                         >  2:1488227/1‑93 (MQ=255)
 aaCGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGcag                                                                                        >  1:936183/1‑93 (MQ=255)
 aaCGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGcag                                                                                        <  1:3062783/93‑1 (MQ=255)
   cGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAgcagca                                                                                      <  2:2440381/93‑1 (MQ=255)
   cGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAgcagca                                                                                      <  1:2345308/93‑1 (MQ=255)
     ttAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGa                                                                                    >  2:3217368/1‑93 (MQ=255)
     ttAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGa                                                                                    >  2:2737530/1‑93 (MQ=255)
      tAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGaa                                                                                   >  2:4592360/1‑93 (MQ=255)
      tAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGaa                                                                                   <  1:2527557/93‑1 (MQ=255)
       aaCGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAg                                                                                  >  1:4163519/1‑93 (MQ=255)
       aaCGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAg                                                                                  <  1:4681642/93‑1 (MQ=255)
          ggTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTAc                                                                               <  1:1633569/93‑1 (MQ=255)
            tAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACaa                                                                             <  1:2155794/93‑1 (MQ=255)
            tAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACaa                                                                             >  2:2927816/1‑93 (MQ=255)
             aaaGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAAc                                                                            <  2:3913222/93‑1 (MQ=255)
              aaGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACg                                                                           >  1:3732755/1‑93 (MQ=255)
               aGCGATTCTCTACGGCGATAAGCTCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGa                                                                          <  1:137428/93‑1 (MQ=255)
               aGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGa                                                                          >  2:4180223/1‑93 (MQ=255)
               aGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGa                                                                          <  1:4467204/93‑1 (MQ=255)
                tcGATTCTCTACGGCTATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAg                                                                         <  2:539749/92‑1 (MQ=255)
                gCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAg                                                                         >  1:4645252/1‑93 (MQ=255)
                   aTTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGaa                                                                      <  1:3359349/93‑1 (MQ=255)
                   aTTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGaa                                                                      >  2:4538183/1‑93 (MQ=255)
                   aTTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGaa                                                                      >  1:2546747/1‑93 (MQ=255)
                    ttCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAac                                                                     <  1:1473734/93‑1 (MQ=255)
                    ttCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAac                                                                     >  1:1200926/1‑93 (MQ=255)
                      ctctACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAacac                                                                   <  1:2284372/93‑1 (MQ=255)
                      ctctACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAacac                                                                   >  1:1825854/1‑93 (MQ=255)
                      ctctACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAacac                                                                   >  1:1322712/1‑93 (MQ=255)
                       tctACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACg                                                                  <  2:4366705/93‑1 (MQ=255)
                        ctACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACgg                                                                 <  1:1507790/93‑1 (MQ=255)
                         tACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGa                                                                <  2:2152646/93‑1 (MQ=255)
                          aCGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGaa                                                               >  1:2204039/1‑93 (MQ=255)
                           cGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAAt                                                              >  2:1799127/1‑93 (MQ=255)
                            ggCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAAtt                                                             >  2:4133590/1‑93 (MQ=255)
                            ggCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAAtt                                                             >  2:3858875/1‑93 (MQ=255)
                             gCGATAAGATCACACCATCAATGACGAAAGCGATTAGCGAAACCTATCGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTa                                                            <  1:1964011/93‑1 (MQ=255)
                               gATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTAcg                                                          >  1:3750717/1‑93 (MQ=255)
                               gATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAACAGCAGAAGTACAACGAGGACCACGGAATTAcc                                                          >  1:1552821/1‑92 (MQ=255)
                                aTAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTAcgc                                                         <  1:4605012/93‑1 (MQ=255)
                                aTAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTAcgc                                                         <  2:2535159/93‑1 (MQ=255)
                                     aTCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCaa                                                    <  1:2051492/93‑1 (MQ=255)
                                        aCCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGc                                                 >  1:2279461/1‑93 (MQ=255)
                                          cccATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCtt                                               >  2:1636370/1‑93 (MQ=255)
                                          cccATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCtt                                               >  2:1303079/1‑93 (MQ=255)
                                          cccATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAACAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCtt                                               >  2:1585873/1‑93 (MQ=255)
                                            cATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGa                                             >  1:4619619/1‑93 (MQ=255)
                                             aTCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGGACAACGAGGAACACGGAATTACGCCGCAAGGCTTGaa                                            >  2:1996610/1‑93 (MQ=255)
                                              tCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAAc                                           >  2:1360441/1‑93 (MQ=255)
                                                aaTGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACaa                                         >  2:1099088/1‑93 (MQ=255)
                                                aaTGGCGAAAGCGATTGGCGAAACCGAACGGCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACaa                                         <  1:1522890/93‑1 (MQ=255)
                                                aaTGCCGAAACCGATTGGCGAAACCGAACGTCGCCGGGAGAAGCAGCAGAAGGACAAGGAGGAACACGGAATTACGCCGCAAGGCTTGAACaa                                         >  2:32546/1‑93 (MQ=255)
                                                  aggCGAAAGCGATTGGCGAAGCCGAACGTCGCCGTGAGAAACAGCAGAAGTACAACGAGGACCACGGAATTACGCCGCAAGGCTTGAACAAGa                                       <  2:3015212/92‑1 (MQ=255)
                                                     cGAACGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAg                                    <  1:1164727/93‑1 (MQ=255)
                                                     cGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAg                                    <  1:3668131/93‑1 (MQ=255)
                                                     cGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAg                                    <  2:834553/93‑1 (MQ=255)
                                                       aaaGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTg                                  >  2:2790411/1‑93 (MQ=255)
                                                       aaaGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTg                                  >  2:3626951/1‑93 (MQ=255)
                                                        aaGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAggag                                 >  2:884561/1‑90 (MQ=255)
                                                        aaGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTgg                                 >  2:745167/1‑93 (MQ=255)
                                                         aGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGt                                <  2:2794461/93‑1 (MQ=255)
                                                         aGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGt                                >  2:4236228/1‑93 (MQ=255)
                                                         aGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGt                                <  1:4465952/93‑1 (MQ=255)
                                                         aGCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGt                                <  1:3987067/93‑1 (MQ=255)
                                                          gCGATTGGCGAAACGGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTc                               <  2:4447484/93‑1 (MQ=255)
                                                          gCGATTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTc                               <  2:1170102/93‑1 (MQ=255)
                                                             aTTGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGat                            >  1:376590/1‑93 (MQ=255)
                                                              ttGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGata                           <  2:2721565/93‑1 (MQ=255)
                                                               tGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGatat                          <  1:388472/93‑1 (MQ=255)
                                                               tGGCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGatat                          >  2:2349544/1‑93 (MQ=255)
                                                                ggCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATc                         <  1:2608387/93‑1 (MQ=255)
                                                                 gCGAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATcc                        <  1:3137021/93‑1 (MQ=255)
                                                                   gAAACCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTg                      <  1:875004/93‑1 (MQ=255)
                                                                    aaaCCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTgg                     <  1:1673665/93‑1 (MQ=255)
                                                                     aaCCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGgc                    <  1:4712581/93‑1 (MQ=255)
                                                                     aaCCGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGgc                    >  1:770109/1‑93 (MQ=255)
                                                                       ccGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGgcgc                  <  1:492228/93‑1 (MQ=255)
                                                                        cGAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCt                 <  1:2417898/93‑1 (MQ=255)
                                                                         gAACGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTg                >  2:336855/1‑93 (MQ=255)
                                                                          aaCGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTgg               >  1:1079287/1‑93 (MQ=255)
                                                                           aCGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTggg              <  2:4502213/93‑1 (MQ=255)
                                                                           aCGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTggg              >  1:1769712/1‑93 (MQ=255)
                                                                            cGTCGCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTgggg             <  1:4408220/93‑1 (MQ=255)
                                                                              tagCCGTGAGAAGCAGCAGAAGTACAACGAGGAACGCGGAATTACGCCGCAAGGTTTGAACAGGAAAGTGGTCGATATCCTGGCGCTGGGGCa           <  1:4721847/91‑1 (MQ=255)
                                                                                gCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGa         >  2:3231099/1‑93 (MQ=255)
                                                                                gCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGa         <  1:565263/93‑1 (MQ=255)
                                                                                gCCGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGa         >  1:1137478/1‑93 (MQ=255)
                                                                                 ccGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTTGCGCTGGCGCAGaa        >  1:1755173/1‑93 (MQ=255)
                                                                                 ccGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGaa        <  2:173757/93‑1 (MQ=255)
                                                                                 ccGTGAGAAACAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGaa        <  1:4700094/93‑1 (MQ=255)
                                                                                  cGTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAAc       >  2:3459278/1‑93 (MQ=255)
                                                                                   gTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACa      <  2:1268426/93‑1 (MQ=255)
                                                                                   gTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACa      >  1:3586810/1‑93 (MQ=255)
                                                                                   gTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACa      >  1:296417/1‑93 (MQ=255)
                                                                                   gTGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACa      >  2:481703/1‑93 (MQ=255)
                                                                                    tGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACAt     >  2:3766720/1‑93 (MQ=255)
                                                                                    tGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACAt     <  2:1650725/93‑1 (MQ=255)
                                                                                    tGAGAAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACAt     <  1:1079492/93‑1 (MQ=255)
                                                                                     gagaAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACAtt    >  1:1274551/1‑93 (MQ=255)
                                                                                       gaAGCAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACATTGc  >  1:2712377/1‑93 (MQ=255)
                                                                                          |                                                                                         
TAACGTTAACGGTAAAGCGATTCTCTACGGCGATAAGATCACCCCATCAATGGCGAAAGCGATTGGCGAAACCGAACGTCGCCGTGAGAAACAGCAGAAGTACAACGAGGAACACGGAATTACGCCGCAAGGCTTGAACAAGAAAGTGGTCGATATCCTGGCGCTGGGGCAGAACATTGC  >  NC_000913/814380‑814559

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