Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,783,402 A→G 53.9% V87A (GTC→GCC)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,783,4020AG53.9% 63.4 / 265.1 211V87A (GTC→GCC) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (58/39);  new base G (66/47);  total (125/86)
Fisher's exact test for biased strand distribution p-value = 8.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 212 total aligned reads displayed.

AACAGACGTCCGGTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATCGGGG  >  NC_000913/3783303‑3783500
                                                                                                   |                                                                                                  
aaCAGACGTCCGGTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGa                                                                                                    >  1:734309/1‑100 (MQ=255)
   aGACGTCCGGTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCAt                                                                                                 >  2:3143740/1‑100 (MQ=255)
         ccGGTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCCCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTgg                                                                                           >  2:3220576/1‑100 (MQ=255)
           ggTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGt                                                                                         <  2:2741349/100‑1 (MQ=255)
            gTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTa                                                                                        >  2:3309701/1‑100 (MQ=255)
              ttCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTacg                                                                                      >  1:3949227/1‑100 (MQ=255)
                 cGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTacgacg                                                                                   >  1:646066/1‑100 (MQ=255)
                 cGCTGCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTacgacg                                                                                   <  1:355777/100‑1 (MQ=255)
                  gCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACAAGGCTGGGTacgacga                                                                                  >  1:4068592/1‑100 (MQ=255)
                     tCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTAcgacgacga                                                                               >  1:3160721/1‑100 (MQ=255)
                      ccAGCCCTTTGGTCGCTCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACgag                                                                              >  2:2615217/1‑100 (MQ=255)
                       cAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACgaga                                                                             >  1:190344/1‑100 (MQ=255)
                       cAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACgaga                                                                             <  2:191390/100‑1 (MQ=255)
                        aGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGaa                                                                            >  2:874621/1‑100 (MQ=255)
                        aGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGCCgaga                                                                             >  2:2095996/1‑99 (MQ=255)
                         gCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAat                                                                           >  2:3118695/1‑100 (MQ=255)
                         gCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAat                                                                           <  2:132127/100‑1 (MQ=255)
                          cccTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAata                                                                          <  2:2204953/100‑1 (MQ=255)
                           ccTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAatat                                                                         >  1:2571171/1‑100 (MQ=255)
                            cTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTTATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATAtt                                                                        <  1:1011334/100‑1 (MQ=255)
                            cTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATAtt                                                                        >  1:1363255/1‑100 (MQ=255)
                               tGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACg                                                                     >  1:4094181/1‑100 (MQ=255)
                                gggcGCCCACACCAGTCGCGCATCAGGACGCATCAGTGGTTTAATCTGGAGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACgg                                                                    <  1:2359027/97‑1 (MQ=255)
                                 gTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACggc                                                                   >  1:1689544/1‑100 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACggctg                                                                 >  2:2010178/1‑100 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACggctg                                                                 >  1:3108468/1‑100 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACggct                                                                  >  2:1203849/1‑99 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACggctg                                                                 <  2:74858/100‑1 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACggctg                                                                 <  2:3017717/100‑1 (MQ=255)
                                   cGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACggctg                                                                 <  2:646060/100‑1 (MQ=255)
                                     cccACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACggctgg                                                                >  2:4101259/1‑99 (MQ=255)
                                       cacaCCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCCCCAAAGACATGGCTGGGTACGACGACGAGAATATTACGgctggctg                                                             >  2:2999355/1‑100 (MQ=255)
                                         cacCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGcc                                                           >  1:2042189/1‑100 (MQ=255)
                                           ccAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAg                                                         >  2:3908265/1‑100 (MQ=255)
                                              gACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAgcgc                                                      <  2:1397657/100‑1 (MQ=255)
                                                  cgcATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTg                                                  <  2:11337/100‑1 (MQ=255)
                                                  cgcATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTg                                                  >  1:1247540/1‑100 (MQ=255)
                                                   gcATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTgg                                                 <  2:3618443/100‑1 (MQ=255)
                                                    caacAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCCAAGGCATGGCTGGGTACGACGACGAGAATATTCCGGCTGGCTGCCAGCGCAGTgg                                                 <  1:2956139/96‑1 (MQ=255)
                                                     aTCAGGACGCGTCAGTGGTTTAAGCTGGCGCAGCACTTCACCAAGGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCg                                               <  2:1234239/100‑1 (MQ=255)
                                                      tCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCga                                              >  1:2256868/1‑100 (MQ=255)
                                                      tCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCga                                              <  1:1294248/100‑1 (MQ=255)
                                                       cAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCGTGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCgag                                             <  1:1502345/100‑1 (MQ=255)
                                                       cAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCgag                                             <  1:1855704/100‑1 (MQ=255)
                                                        aGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCgaga                                            >  2:351906/1‑100 (MQ=255)
                                                        aGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCgaga                                            >  2:571226/1‑100 (MQ=255)
                                                         ggACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGAt                                           >  1:1135912/1‑100 (MQ=255)
                                                          gACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATc                                          <  1:2061304/100‑1 (MQ=255)
                                                            cGCATCAGTGGTTTAATCTGGCGCAGCACGTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGc                                        <  2:542664/100‑1 (MQ=255)
                                                            cGCATCAGTGGTTTAACCGGGCGCAGCACTACACCAAAGGCATGGCTGAGTCCGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCg                                               >  2:513856/1‑93 (MQ=255)
                                                             gCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCt                                       >  2:3832572/1‑100 (MQ=255)
                                                              cATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCtt                                      >  2:2239083/1‑100 (MQ=255)
                                                              cATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCtt                                      >  1:434810/1‑100 (MQ=255)
                                                              cATCAGAGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGCCTGGGTACGACGACGAACATATCACGGCCAGCTGCCAGCGCAGTGGCAAGATcac                                        >  2:1271796/1‑96 (MQ=255)
                                                               aTCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCttt                                     <  2:3680644/100‑1 (MQ=255)
                                                                 cAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCa                                   >  1:3284628/1‑100 (MQ=255)
                                                                 cAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCa                                   >  1:1113015/1‑100 (MQ=255)
                                                                  aGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCaa                                  >  1:1353292/1‑100 (MQ=255)
                                                                  aGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCaa                                  <  1:2814372/100‑1 (MQ=255)
                                                                  aGTGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCaa                                  >  2:3982150/1‑100 (MQ=255)
                                                                   gTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAg                                 >  1:2090170/1‑100 (MQ=255)
                                                                    tGGTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGa                                <  2:2647902/100‑1 (MQ=255)
                                                                     ggTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAAAATATTACGGCTGGCTGCAAGCGCAGTGGCGAGATCGCTTTCAAGAt                               >  1:3078950/1‑100 (MQ=255)
                                                                     ggTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGAt                               >  1:2926273/1‑100 (MQ=255)
                                                                      gTTTAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATg                              >  2:1637342/1‑100 (MQ=255)
                                                                      gTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGAt                               >  2:1800728/1‑99 (MQ=255)
                                                                      gTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATg                              >  1:951300/1‑100 (MQ=255)
                                                                      gTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATg                              >  2:3766598/1‑100 (MQ=255)
                                                                       tttAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATgg                             <  1:2183043/100‑1 (MQ=255)
                                                                       tttAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATgg                             <  1:548365/100‑1 (MQ=255)
                                                                         tAATCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTAAGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCAAGATAGCTTTAAAGATGGAg                           >  1:342965/1‑100 (MQ=255)
                                                                         tAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAg                           <  1:3816183/100‑1 (MQ=255)
                                                                          aaTCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGc                          >  1:3721966/1‑100 (MQ=255)
                                                                          aaTCTGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGc                          <  1:3513465/100‑1 (MQ=255)
                                                                           aTCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCg                         >  1:2665227/1‑100 (MQ=255)
                                                                           aTCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCg                         >  2:3073466/1‑100 (MQ=255)
                                                                           aTCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCg                         <  1:3676448/100‑1 (MQ=255)
                                                                            tCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGt                        >  2:1336695/1‑100 (MQ=255)
                                                                              tGGCGCAGCACTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTAt                      <  1:3908291/100‑1 (MQ=255)
                                                                              tGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTAt                      >  1:3139564/1‑100 (MQ=255)
                                                                               ggCGCAGCACTTCACCAAAGACTTGACTGGGTACGACGACAAGAATATAACGGCAGGCCGCCAGCGAAGTGACAAGATCCATTTGAAGATGGacgtcatc                     >  2:3636089/1‑93 (MQ=255)
                                                                                  gcagcaCTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCggg                  >  1:499830/1‑100 (MQ=255)
                                                                                  gcagcaCTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCggg                  <  2:1668608/100‑1 (MQ=255)
                                                                                     gcaCTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGaaa               >  1:3010712/1‑100 (MQ=255)
                                                                                     gcaCTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGaaa               <  1:1939037/100‑1 (MQ=255)
                                                                                       aCTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAg             >  1:374972/1‑100 (MQ=255)
                                                                                       aCTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAg             >  2:2944249/1‑100 (MQ=255)
                                                                                       aCTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGGCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAg             <  2:826096/100‑1 (MQ=255)
                                                                                       aCTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAg             >  2:2415428/1‑100 (MQ=255)
                                                                                        cTTCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTTCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAgg            >  1:3854243/1‑100 (MQ=255)
                                                                                         ttCACCAAAGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGc           >  1:881902/1‑100 (MQ=255)
                                                                                          tCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGca          >  2:1990588/1‑100 (MQ=255)
                                                                                           cACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGca          <  1:3643735/99‑1 (MQ=255)
                                                                                            aCCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGcaca        <  2:3813797/100‑1 (MQ=255)
                                                                                             ccAAAGACATGGCTGGGTACGACGACGAGACCATTACGGCTGGCTCGCAGCGCAGTGCCGAGAACGCTATCAAGCAGGAGCGGACCGGGAAAAGGcac         >  2:44847/1‑98 (MQ=255)
                                                                                             ccAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTACCGGGAAAAGGcaca        >  2:742394/1‑99 (MQ=255)
                                                                                              cAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATc      >  2:3301026/1‑100 (MQ=255)
                                                                                               aaaGGCATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATCg     >  1:2822163/1‑100 (MQ=255)
                                                                                               aaaGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATCg     >  2:3438148/1‑100 (MQ=255)
                                                                                                  gACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATCggg   <  1:3244852/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
AACAGACGTCCGGTTTCCGCTTCCAGCCCTTTGGTCGCCCACACCAGACGCGCATCAGGACGCATCAGTGGTTTAATCTGGCGCAGCACTTCACCAAAGACATGGCTGGGTACGACGACGAGAATATTACGGCTGGCTGCCAGCGCAGTGGCGAGATCGCTTTCAAGATGGAGCGTATCGGGAAAAGGCACATCGGGG  >  NC_000913/3783303‑3783500

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: