Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 15.2% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC15.2% 382.3 / 64.9 163intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (56/82);  new base C (12/13);  total (68/95)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 166 total aligned reads displayed.

ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060165‑4060363
                                                                                                   |                                                                                                    
aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                      >  2:957823/1‑100 (MQ=255)
aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                      >  2:1556208/1‑100 (MQ=255)
 ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                     >  1:1917853/1‑100 (MQ=255)
  tCGTCACCGTCACCCGCCGCACACCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                    <  2:710307/100‑1 (MQ=255)
    gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                  <  1:3811021/100‑1 (MQ=255)
     tgaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 >  2:3541681/3‑100 (MQ=255)
     tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGcn                                                                                                 <  2:142715/100‑2 (MQ=255)
     tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 >  1:3711615/1‑100 (MQ=255)
       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                               <  2:1944462/100‑1 (MQ=255)
       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAACCTGCCATCGCTGGCAgn                                                                                               <  1:3148375/99‑2 (MQ=255)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                             <  2:3087067/100‑1 (MQ=255)
            cACCTGACGGAAAACGATCGGCGCCGGGCGAACCGCGCACCGAAAGACGATTAATGTCGGTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttt                                                                                         >  1:2222994/1‑100 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCAGCGATGGCAGG‑ttttttt                                                                                       <  1:3139015/100‑1 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttttta                                                                                       <  2:3336210/100‑2 (MQ=255)
               cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttat                                                                                      <  2:1239318/100‑3 (MQ=255)
               cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttat                                                                                      <  2:276309/100‑3 (MQ=255)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatg                                                                                     <  1:2777578/100‑4 (MQ=255)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatg                                                                                     <  2:549097/100‑4 (MQ=38)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                     <  2:3609588/100‑1 (MQ=255)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatga                                                                                    <  1:289579/100‑5 (MQ=38)
                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                                   <  1:565867/100‑1 (MQ=255)
                   cGGAAAACGATCGTCGCTGCGTAAAAAGCGCACCGAAAGAGGCTGAATTTCATTTTTTAGTAATAAAAAAGCTGCGAGCGATGGCAGG‑atttatttaga                                                                                   >  1:3432466/1‑96 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                  >  2:2817057/1‑100 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                  >  2:3032854/1‑100 (MQ=255)
                      aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGt                                                                               <  2:3192365/100‑1 (MQ=255)
                      aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGt                                                                               <  2:3506676/100‑1 (MQ=255)
                       aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                              >  2:1165537/1‑100 (MQ=255)
                            aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAg                                                                         >  2:2256427/1‑100 (MQ=255)
                            aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTGTTTGCCTGTCGTTAg                                                                         <  1:777012/100‑1 (MQ=255)
                             tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgc                                                                        >  2:2389951/1‑100 (MQ=255)
                              cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcg                                                                       <  2:2868967/100‑1 (MQ=255)
                              cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGAGATGGCAGG‑TTTTTTTTTGACTGTCGTGAgcg                                                                       >  1:535184/1‑100 (MQ=255)
                               gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcgc                                                                      <  2:3592854/100‑1 (MQ=255)
                               atcgCCGCGCCAACCGCGCACCGAAAGACGATTAGTTTCGTTGTTTATTCATATAAAACCTGCTAGCGAGGGCAGT‑TTATTTTTTGctt                                                                                 >  1:2612619/2‑87 (MQ=255)
                                tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAgcgc                                                                      <  2:2078510/100‑1 (MQ=255)
                                 cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcan                                                                    <  2:174644/100‑2 (MQ=255)
                                  gtcgagcCAAGCGAGCACAGTAAGACGATTAATTTCGTTTTTTAGTTATTAAAAACCTGCTAGCTATGGTAGG‑TTTTTTTTTGTCGGTCGTGa                                                                          <  1:1393634/88‑1 (MQ=255)
                                  gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                   >  2:3186014/1‑100 (MQ=255)
                                    cgcgCCAACCGCGCACCGAAATACTATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 <  1:2802683/100‑1 (MQ=255)
                                    cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 <  2:1352901/100‑1 (MQ=255)
                                    cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 <  1:2626827/100‑1 (MQ=255)
                                       gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                              >  1:1767165/1‑100 (MQ=255)
                                        ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                             <  2:2623779/100‑1 (MQ=255)
                                        ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                             <  1:2256416/100‑1 (MQ=255)
                                          aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                           <  2:694118/100‑1 (MQ=255)
                                          aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                           <  2:3065539/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  1:3240022/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  2:2499746/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  2:2367936/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         >  1:2737438/1‑100 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         >  1:424070/1‑100 (MQ=255)
                                                gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                     <  2:683268/100‑1 (MQ=255)
                                                   ccGAAAGACGATTAATTTCGTTTTTTTTTTTTTTTTAACCTGCCAGCGATGGCAGG‑TTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTCTCTcc                                                  >  1:1464742/1‑100 (MQ=255)
                                                   ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                  >  2:807410/1‑100 (MQ=255)
                                                    cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                 <  2:2098877/100‑1 (MQ=255)
                                                       aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                              <  2:2960332/100‑1 (MQ=255)
                                                         gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                            >  2:1775961/1‑100 (MQ=255)
                                                         gACGATTAATTTCGTTGTTTAGTCATAAAAAACGTGCCGGCGATGGCAGG‑TTTTTTGTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                            <  1:3512559/100‑1 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        <  2:3210389/100‑1 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                         >  1:2413282/1‑99 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        <  2:2132758/100‑1 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        <  2:2464370/100‑1 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        <  1:1956478/100‑1 (MQ=255)
                                                                 aTTTCGTTTTTTTGTTATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACGCCAATGCGAGTATCTCCCTTATAAGCCtgt                                     >  1:2576394/1‑99 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGGCTGTCGTGAGCGGACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    >  2:3171528/1‑100 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    <  2:557441/100‑1 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGTAGG‑TTTTTTTTTGACTGCCGCGAGAGACAAGCAATGCGAGTATCTCACatatcaatgtgcga                                   >  1:3621430/1‑90 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   >  2:1363074/1‑100 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   <  1:1539502/100‑1 (MQ=255)
                                                                       ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                              <  2:1015915/100‑1 (MQ=255)
                                                                       tcttttAGTCATAAAAAACTTGCCAGGGATGGTATG‑TTTTTTTTTGCCTGTCTTGAGCGCACACCAATGCGAGGATCGCCCTTAGAAGTCTGTGGTTa                                <  1:29785/96‑1 (MQ=255)
                                                                         ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCTTGTCGTGAGCGCCCACCAATGCGAGTATCTCCCTCAGATGTCTGTGGTTACTTc                            >  1:3701283/1‑100 (MQ=255)
                                                                         ttccgatcttAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                            <  2:3562120/91‑1 (MQ=255)
                                                                           gatctCATAAAAAACGTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                         <  2:3762394/96‑1 (MQ=255)
                                                                            tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                         >  2:456596/1‑100 (MQ=255)
                                                                            tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                         <  2:2889045/100‑1 (MQ=255)
                                                                               caaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                      >  2:589784/5‑100 (MQ=38)
                                                                                cATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagatc                    >  1:3762370/1‑98 (MQ=255)
                                                                                aaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATTTCCCTTATAAGTCTGTGGTTAAGTCAATa                        >  1:3574161/4‑97 (MQ=37)
                                                                                aaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                     >  2:1405323/4‑100 (MQ=38)
                                                                                 aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                    >  2:1721013/3‑100 (MQ=255)
                                                                                 aTAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                     <  2:1874495/99‑1 (MQ=255)
                                                                                 aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                    >  2:3726260/1‑100 (MQ=255)
                                                                                  tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCagatcggaa                   >  1:3562099/1‑92 (MQ=255)
                                                                                  tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                   >  2:2684275/1‑100 (MQ=255)
                                                                                   acaaaaCCTGCCAACGATGGCAGG‑TTTTTTTCTGCCTGTCGTAAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                  >  2:2573297/3‑100 (MQ=255)
                                                                                   aaaaaaCTTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCGCACCAATGCGTGTAGCTCCCTTATGAGTCGGTGGTTACGTCAATAGAGAGc                  <  1:199960/100‑1 (MQ=255)
                                                                                   aaaaaaCCTGCCAGCGATGGCAGT‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCAAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                   <  1:348156/99‑1 (MQ=255)
                                                                                     aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTTCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTTTGTGGTTACTTCAATAGAGAGAtt                >  1:1164229/1‑100 (MQ=255)
                                                                                      aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt               >  2:1510873/1‑100 (MQ=255)
                                                                                       aaCCTGCCATCGCTGGCAGG‑TTTTTTTTTCCTTGTCGTCAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAgcgc                  >  2:3243849/1‑96 (MQ=38)
                                                                                       aaCCTGCCAGCGATGGCAGG‑TTTTTTTGTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATCAGTCTGTGGTTACGTCAATAGAGAGCTTTg              <  1:379359/100‑1 (MQ=255)
                                                                                         ccTGCCAGCGATGGCAGG‑TTTTTTTTTTTTTTTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc            >  1:3099851/1‑100 (MQ=255)
                                                                                          cTGCCAGCGAGGGCAGG‑TTTTTTTTTGCCTGTCTTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg              <  1:1072555/97‑1 (MQ=255)
                                                                                           tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac          <  1:879995/100‑1 (MQ=255)
                                                                                             ccAGCGATGGCAGGTTTTTTTTTTTTTTTTCTTGATGGCATTTCAATGCTAGTATCTCGCTTATAATATTGTTTTTACGTCAATAAAGAACTTTTTcaca         >  1:665690/1‑100 (MQ=255)
                                                                                              cAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt       <  2:812671/100‑1 (MQ=255)
                                                                                              cAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt       <  1:1115972/100‑1 (MQ=255)
                                                                                               aGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt      >  2:747982/1‑100 (MQ=255)
                                                                                                 cGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTTTCTCCCTTTTAAGTCTGTGGTTACGACAATAAAGAGCTTTGTCACAATTAt    >  1:1294125/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060165‑4060363

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