Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,006,859 A→C S200R (AGC→CGC)  ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,006,8590AC98.2% 790.9 / 0.3 228S200R (AGC→CGC) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (110/113);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 228 total aligned reads displayed.

CAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGAC  >  NC_000913/3006763‑3006954
                                                                                                |                                                                                               
 aGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAGAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCt                                                                                              >  2:5090383/1‑99 (MQ=255)
 aGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCt                                                                                              >  2:428178/1‑99 (MQ=255)
  gTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTa                                                                                             <  2:2293113/99‑1 (MQ=255)
  gTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTa                                                                                             <  1:4629470/99‑1 (MQ=255)
  gTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTa                                                                                             <  2:1201068/99‑1 (MQ=255)
   tCAATGGCTGACCTCGCGTGGTTACGTGCACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTAt                                                                                            >  1:3548637/1‑99 (MQ=255)
   tCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTAt                                                                                            >  1:1202409/1‑99 (MQ=255)
   tCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTAt                                                                                            >  2:738027/1‑99 (MQ=255)
   tCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTAt                                                                                            <  1:325680/99‑1 (MQ=255)
    cAATGGCTGACCCCGCGTGGTTACGTGAACACCTTCGTTCACTGGAACACCTGAAAAGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCGATg                                                                                           <  2:3781534/99‑1 (MQ=255)
     aaTGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTAGACCTACTCTCCACGCTATgg                                                                                          >  1:3537518/1‑99 (MQ=255)
         gCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCaaa                                                                                      >  1:729227/1‑99 (MQ=255)
         gCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCaaa                                                                                      <  2:3593263/99‑1 (MQ=255)
            gACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCg                                                                                   <  1:2589560/99‑1 (MQ=255)
              ccTCGCGTGGTTACGTGAACACTTTGGTTCTCTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGct                                                                                 <  2:3409386/99‑1 (MQ=255)
               cTCGCGTGGTTACGTGAACGCTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctc                                                                                >  1:4510653/1‑99 (MQ=255)
               cTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctc                                                                                <  1:1055355/99‑1 (MQ=255)
                tCGCGTTGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                               <  1:2376760/99‑1 (MQ=255)
                tCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                               <  2:2605247/99‑1 (MQ=255)
                tCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                               >  2:2643090/1‑99 (MQ=255)
                tCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                               >  2:2940697/1‑99 (MQ=255)
                tCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctct                                                                               >  1:1336196/1‑99 (MQ=255)
                tCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGAGCCTACTCTCCAAGCTATGGCAAACCGctct                                                                               <  2:4643476/99‑1 (MQ=255)
                 cgcgTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctctc                                                                              <  2:2708002/99‑1 (MQ=255)
                  gcgTGGTTACGTGAACACTTTGGTTCATTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctctct                                                                             <  1:1787222/99‑1 (MQ=255)
                  gcgTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctctct                                                                             <  2:414865/99‑1 (MQ=255)
                  gcgTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGctctct                                                                             >  1:3501380/1‑99 (MQ=255)
                   cgTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTg                                                                            <  2:4834798/99‑1 (MQ=255)
                   cgTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTg                                                                            <  2:2020971/99‑1 (MQ=255)
                        ttACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACca                                                                       <  2:1374025/99‑1 (MQ=255)
                        ttACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACca                                                                       <  2:3863284/99‑1 (MQ=255)
                            gTGAACACTTTGGTTCACTGGAAACCCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                   <  1:3316001/99‑1 (MQ=255)
                            gTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                   >  2:838265/1‑99 (MQ=255)
                            gTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTAATCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                   >  1:2200221/1‑99 (MQ=255)
                            gTGAACACTTTGGCTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAg                                                                   >  2:3948650/1‑99 (MQ=255)
                              gAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGc                                                                 >  1:343256/1‑99 (MQ=255)
                               aaCACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGca                                                                <  1:4813113/99‑1 (MQ=255)
                               aaCACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGAACTGGGCCTACTATACACGCTATGGCACACCGCTCTCTGTACCACAAGGca                                                                >  2:1277936/1‑99 (MQ=255)
                                  acTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGcatca                                                             <  1:3158632/99‑1 (MQ=255)
                                  acTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGcatca                                                             >  2:3153322/1‑99 (MQ=255)
                                   cTTTGGTTCATTCGCAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAAGcatcat                                                            <  2:1802/99‑1 (MQ=255)
                                    tttGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCatcatc                                                           >  2:462023/1‑99 (MQ=255)
                                    tttGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCatcatc                                                           >  2:3535739/1‑99 (MQ=255)
                                     ttGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCg                                                          <  1:2482701/99‑1 (MQ=255)
                                       ggTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGTCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGt                                                        <  1:4078940/99‑1 (MQ=255)
                                        gTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGACAAACCGCTCTCTGTACCACAAGGCATCATCGGTc                                                       >  1:2965804/1‑99 (MQ=255)
                                         ttCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCt                                                      <  2:1899554/99‑1 (MQ=255)
                                         ttCACTGGAAAACCTGAAAGGTAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtg                                                     >  2:1524382/1‑99 (MQ=255)
                                            aCTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtgat                                                   >  1:2095919/1‑99 (MQ=255)
                                            aCTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtgat                                                   <  1:4882274/99‑1 (MQ=255)
                                            aCTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtgat                                                   <  1:818438/99‑1 (MQ=255)
                                              tGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtgatga                                                 >  2:3715129/1‑99 (MQ=255)
                                              tGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCtgatga                                                 >  1:2927122/1‑99 (MQ=255)
                                              tGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGAAAAAAGCTCTATGTACCACAAGGCATAATCGGTCtgatga                                                 >  1:3925669/1‑99 (MQ=255)
                                               ggAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGAc                                                <  1:4782749/99‑1 (MQ=255)
                                                gAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACt                                               >  1:1004683/1‑99 (MQ=255)
                                                 aaaaCCTGAAAGGTAAAAACATCGCCATGAACTGGGCATACTCTACACGCTATAGCAAACCGCTCTCTGTACAACAAGGCATCATCGGTATGATGAccc                                              <  2:1033367/99‑3 (MQ=255)
                                                   aaCCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGc                                            <  1:5230474/99‑1 (MQ=255)
                                                   aaCCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGc                                            <  1:1667820/99‑1 (MQ=255)
                                                    aCCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCt                                           <  2:3340910/99‑1 (MQ=255)
                                                    aCCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCt                                           <  2:4483313/99‑1 (MQ=255)
                                                      cTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTc                                         >  1:961093/1‑99 (MQ=255)
                                                      cTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTc                                         >  2:4049832/1‑99 (MQ=255)
                                                      cTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTc                                         >  1:1615220/1‑99 (MQ=255)
                                                       tGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCg                                        >  1:1232045/1‑99 (MQ=255)
                                                         aaaGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGt                                      >  2:1602302/1‑99 (MQ=255)
                                                          aaGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGtg                                     >  2:2447013/1‑98 (MQ=255)
                                                            ggTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATg                                   <  1:4895525/99‑1 (MQ=255)
                                                            ggTAAAAAAATCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATg                                   <  1:4291563/99‑1 (MQ=255)
                                                                aaaaaaTCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATg                               >  1:1091516/1‑99 (MQ=255)
                                                                    aaTCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCAc                           >  2:3009821/1‑99 (MQ=255)
                                                                     aTCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCAcc                          <  2:2258504/99‑1 (MQ=255)
                                                                      tCGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCAccc                         <  2:2501855/99‑1 (MQ=255)
                                                                       ctccATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCTCAAGGCATCATCGGTCTGATTACTCGCTTCGGTATGGATGTCACCCt                        <  2:5145269/97‑1 (MQ=255)
                                                                       cGCCATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCt                        <  2:1581349/99‑1 (MQ=255)
                                                                         ccATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTgg                      >  1:1445906/1‑99 (MQ=255)
                                                                         ccATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATATCACCCTgg                      >  1:2178185/1‑99 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTGACCCTGGc                     <  1:1565230/99‑1 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGc                     >  1:885658/1‑99 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGc                     >  1:1128259/1‑99 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGc                     <  1:1846718/99‑1 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGc                     <  1:1922515/99‑1 (MQ=255)
                                                                          cATGACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGc                     <  1:2254464/99‑1 (MQ=255)
                                                                           aTGAGCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTATACCACAAGGCATCATCGGGCTGATGACTCGCCTCGGTATGGATTGCACCCGGGcc                    >  2:2341723/1‑99 (MQ=255)
                                                                             gACCTGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCa                  >  1:4208322/1‑99 (MQ=255)
                                                                                 tGGGCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCg              >  2:4837684/1‑99 (MQ=255)
                                                                                  gggCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCgg             >  2:3469674/1‑99 (MQ=255)
                                                                                    gCCTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATTACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGaa           <  1:106906/99‑1 (MQ=255)
                                                                                     ccTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAg          <  1:3789472/99‑1 (MQ=255)
                                                                                     ccTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAg          >  1:3367558/1‑99 (MQ=255)
                                                                                      cTACTCTCCACGCTATGGCAAGCCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAgg         >  2:1978937/1‑99 (MQ=255)
                                                                                      cTACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAgg         <  1:4008095/99‑1 (MQ=255)
                                                                                       tACTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGc        <  2:938170/99‑1 (MQ=255)
                                                                                        aCTCTCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCt       >  1:3064314/1‑99 (MQ=255)
                                                                                         ctctCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACTCTGGCCCATCCGGAAGGCTa      >  1:3951543/1‑99 (MQ=255)
                                                                                         ctctCCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTa      <  2:1081915/99‑1 (MQ=255)
                                                                                          tctcCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTAc     <  1:555110/99‑1 (MQ=255)
                                                                                            tcCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGa   >  2:2862369/1‑99 (MQ=255)
                                                                                            tcCACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGa   <  2:3279163/99‑1 (MQ=255)
                                                                                             ccACGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGAc  >  2:1219559/1‑99 (MQ=255)
                                                                                                |                                                                                               
CAGTCAATGGCTGACCTCGCGTGGTTACGTGAACACTTTGGTTCACTGGAAAACCTGAAAGGTAAAAAAATCGCCATGACCTGGGCCTACTCTCCAAGCTATGGCAAACCGCTCTCTGTACCACAAGGCATCATCGGTCTGATGACTCGCTTCGGTATGGATGTCACCCTGGCCCATCCGGAAGGCTACGAC  >  NC_000913/3006763‑3006954

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