Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,271,825 T→A 5.6% intergenic (‑55/+42) garK ← / ← garR glycerate kinase I/tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,271,8250TA 5.6% 330.1 / 7.7 90intergenic (‑55/+42)garK/garRglycerate kinase I/tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base A (3/2);  total (46/44)
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.54e-01

AATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCC  >  NC_000913/3271681‑3271969
                                                                                                                                                |                                                                                                                                                
aaTTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCa                                                                                                                                              <  3:570107/149‑1 (MQ=255)
  tttCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAgg                                                                                                                                            >  2:1336624/1‑149 (MQ=255)
        gAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAtt                                                                                                                                                 <  3:673917/138‑1 (MQ=255)
        gAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAtt                                                                                                                                                 >  4:673917/1‑138 (MQ=255)
        gAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGa                                                                                                                                      <  1:1336624/149‑1 (MQ=255)
         aaaTCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGa                                                                                                                                      >  4:323034/1‑148 (MQ=255)
         aaaTCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGa                                                                                                                                      <  3:323034/148‑1 (MQ=255)
               tttttCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAgg                                                                                                                                            <  2:73418/136‑1 (MQ=255)
               tttttCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAgg                                                                                                                                            >  1:73418/1‑136 (MQ=255)
                         gCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAAt                                                                                                                               >  3:1314596/1‑139 (MQ=255)
                         gCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAAt                                                                                                                               <  4:1314596/139‑1 (MQ=255)
                              cgcAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGc                                                                                                                                           >  3:626877/1‑122 (MQ=255)
                              cgcAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGc                                                                                                                                           <  4:626877/122‑1 (MQ=255)
                                   ccTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGt                                                                                                           >  3:874516/1‑149 (MQ=255)
                                         tGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGccc                                                                                                                                         >  4:1329342/1‑113 (MQ=255)
                                         tGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGccc                                                                                                                                         <  3:1329342/113‑1 (MQ=255)
                                          ggCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAAc                                                                                                    <  4:874516/149‑1 (MQ=255)
                                                     tttctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCaa                                                                                                                                <  1:471367/110‑1 (MQ=255)
                                                     tttctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCaa                                                                                                                                >  2:471367/1‑110 (MQ=255)
                                                     tttctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAATCCAGGc                                                                                                                                           <  1:40122/99‑1 (MQ=255)
                                                     tttctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGAATCCAGGc                                                                                                                                           >  2:40122/1‑99 (MQ=255)
                                                       tctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATa                                                                                                                              <  4:1420335/110‑1 (MQ=255)
                                                       tctttATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATa                                                                                                                              >  3:1420335/1‑110 (MQ=255)
                                                             ttagagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCGCATGACGCCGGAGGGCGTTTCTCTTATTCGGCATGGATTCCAGGCCCGGAtt                                                                                                                                    <  4:1242015/96‑1 (MQ=255)
                                                             taAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGAtt                                                                                                                                    >  3:1242015/1‑98 (MQ=255)
                                                               agagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTaa                                                                                                     >  1:450174/1‑127 (MQ=255)
                                                               agagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTaa                                                                                                     <  2:450174/127‑1 (MQ=255)
                                                                gagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGc                                                                                                                                  <  2:108209/97‑1 (MQ=255)
                                                                gagTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGc                                                                                                                                  >  1:108209/1‑97 (MQ=255)
                                                                            aTTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGAtt                                                                                                                                    >  1:602635/1‑83 (MQ=255)
                                                                            aTTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGAtt                                                                                                                                    <  2:602635/83‑1 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                  >  1:1345635/1‑110 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                  <  2:1345635/110‑1 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGcc                                                           >  3:1419154/1‑149 (MQ=255)
                                                                                   ttttCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGcc                                                           >  4:1050465/1‑149 (MQ=255)
                                                                                        aTACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                  <  4:885058/105‑1 (MQ=255)
                                                                                        aTACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGa                                                                                                  >  3:885058/1‑105 (MQ=255)
                                                                                        aTACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTc                                                                                   <  4:374965/120‑1 (MQ=255)
                                                                                        aTACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTc                                                                                   >  3:374965/1‑120 (MQ=255)
                                                                                         tACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAgg                                                        >  2:33748/1‑146 (MQ=255)
                                                                                         tACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAgg                                                        <  1:33748/146‑1 (MQ=255)
                                                                                            cTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTa                                                  >  1:145769/1‑149 (MQ=255)
                                                                                             tATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTAt                                                 >  1:24359/1‑149 (MQ=255)
                                                                                                gCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGAt                                              <  2:145769/149‑1 (MQ=255)
                                                                                                 ccTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgta                                                                     >  2:1192630/1‑125 (MQ=255)
                                                                                                 ccTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgta                                                                     <  1:1192630/125‑1 (MQ=255)
                                                                                                  cTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgta                                                                     <  4:894342/124‑1 (MQ=255)
                                                                                                  cTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgta                                                                     >  3:894342/1‑124 (MQ=255)
                                                                                                   tGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTc                                                                        >  1:1442461/1‑120 (MQ=255)
                                                                                                   tGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTc                                                                        <  2:1442461/120‑1 (MQ=255)
                                                                                                         cacaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGgcgc                                                    <  2:36384/134‑1 (MQ=255)
                                                                                                         cacaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGgcgc                                                    >  1:36384/1‑134 (MQ=255)
                                                                                                         cacaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                     >  1:467996/1‑149 (MQ=255)
                                                                                                          acaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGa                                                                                        >  1:1264034/1‑97 (MQ=255)
                                                                                                          acaTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGa                                                                                        <  2:1264034/97‑1 (MQ=255)
                                                                                                            aTGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAAGACGCCATCAGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTa                                  >  1:1107498/1‑149 (MQ=255)
                                                                                                              gACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCAtt                                                                                                       >  4:294541/1‑78 (MQ=255)
                                                                                                              gACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCAtt                                                                                                       <  3:294541/78‑1 (MQ=255)
                                                                                                              gACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTaaa                                <  2:467996/149‑1 (MQ=255)
                                                                                                               aCGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAAc                               <  3:1050465/149‑1 (MQ=255)
                                                                                                                 gCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGtt                                                                            <  3:624833/102‑1 (MQ=255)
                                                                                                                 gCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGtt                                                                            >  4:624833/1‑102 (MQ=255)
                                                                                                                  ccGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCAt                            <  2:24359/149‑1 (MQ=255)
                                                                                                                    ggAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                     <  3:9469/138‑1 (MQ=255)
                                                                                                                    ggAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                     >  4:9469/1‑138 (MQ=255)
                                                                                                                         gCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATc                           <  1:359918/143‑1 (MQ=255)
                                                                                                                         gCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATc                           >  2:359918/1‑143 (MQ=255)
                                                                                                                           gTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGtt                                      <  1:1009582/130‑1 (MQ=255)
                                                                                                                           gTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGtt                                      >  2:1009582/1‑130 (MQ=255)
                                                                                                                           gTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCa                   >  1:949682/1‑149 (MQ=255)
                                                                                                                            tttCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCAGGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAg                  >  2:105753/1‑149 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGcc                                                                                                                          <  1:1150248/41‑1 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTc                                                                        >  4:1313567/1‑91 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTc                                                                        <  3:1313567/91‑1 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGcc              <  1:105753/149‑1 (MQ=255)
                                                                                                                                tctTATTCGGCCTGGATACCAGGCCCGGATTGCAATACGcc                                                                                                                          >  2:1150248/1‑41 (MQ=255)
                                                                                                                                 ctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                     <  2:1398678/125‑1 (MQ=255)
                                                                                                                                 ctTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTc                                     >  1:1398678/1‑125 (MQ=255)
                                                                                                                                    aTTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCAc                                                              <  2:833664/97‑1 (MQ=255)
                                                                                                                                    aTTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCAc                                                              >  1:833664/1‑97 (MQ=255)
                                                                                                                                     ttCGGCCTGGATTCCATGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGcat         <  2:1107498/149‑1 (MQ=255)
                                                                                                                                     ttCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGcat         >  1:467774/1‑149 (MQ=255)
                                                                                                                                       gggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgtagta                                                                  <  2:1444681/89‑1 (MQ=255)
                                                                                                                                       gggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCgtagta                                                                  >  1:1444681/2‑90 (MQ=255)
                                                                                                                                       gggCCTGGAATCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGcatca       >  2:859795/2‑149 (MQ=255)
                                                                                                                                         gCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCatcatc     <  2:949682/149‑1 (MQ=255)
                                                                                                                                           cTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCt                          >  1:1062220/1‑126 (MQ=255)
                                                                                                                                           cTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCt                          <  2:1062220/126‑1 (MQ=255)
                                                                                                                                               aTTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAtct    >  4:241100/1‑144 (MQ=255)
                                                                                                                                               aTTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAtct    <  1:557379/144‑1 (MQ=255)
                                                                                                                                               aTTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAtct    <  3:241100/144‑1 (MQ=255)
                                                                                                                                               aTTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCAtct    >  2:557379/1‑144 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTaaa                                <  4:648215/115‑1 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTaaa                                >  3:648215/1‑115 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTcc  <  1:885483/145‑1 (MQ=255)
                                                                                                                                                ttCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTcc  >  2:885483/1‑145 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTATGCCTGTTACCACATGACGCCGGAGGGCGTTTCTCTTATTCGGCCTGGATTCCAGGCCCGGATTGCAATACGCCATCCGGGCACGACGTCATTAACGAGTAACTTCGACTTTCGCCAGTTTTTCGTAGTAGCACGCCAGGGCGCTATGATCCGCCGTTCCTAAACCATCTGCTCGCAGTGCCTGCATCATCTCC  >  NC_000913/3271681‑3271969

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