Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,671 A→T 9.0% intergenic (+22/‑314) ygaQ → / → csiD uncharacterized protein; putative enzyme/carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,6710AT 9.0% 334.9 / 15.2 121intergenic (+22/‑314)ygaQ/csiDuncharacterized protein; putative enzyme/carbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (57/53);  new base T (5/6);  total (62/59)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 121 total aligned reads displayed.

CGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGC  >  NC_000913/2788524‑2788819
                                                                                                                                                   |                                                                                                                                                    
cGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATa                                                                                                                                                    >  1:789625/1‑150 (MQ=255)
   aaaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaaa                                                                                                                                                 >  2:629455/1‑150 (MQ=255)
   aaaaTTATACCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGACCCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaaa                                                                                                                                                 >  1:590879/1‑150 (MQ=255)
     aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACg                                                                                                                                               <  1:181693/150‑1 (MQ=255)
     aaTTAACCCACCTGCTAAAGAATATACGCTTAGAGGCGGAACAATAAATATCAATATTCCTGGTCGTAGCGCTATTATTTTAGCGAAAACAGGAGAACCACAGAACAATCTCTATATAAAATTTAATGATAAATATGCCCCGATAAAAcc                                                                                                                                               >  1:987624/1‑149 (MQ=255)
         aTCCCCCCTGGTATAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                           <  1:1469264/150‑1 (MQ=255)
           ccccccTGGTAAATTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAg                                                                                                                                         <  2:1272098/150‑1 (MQ=255)
           ccccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAg                                                                                                                                         >  1:579751/1‑150 (MQ=255)
            cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGa                                                                                                                                        <  1:602113/150‑1 (MQ=255)
            cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                        >  2:254699/1‑150 (MQ=255)
             ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGat                                                                                                                                       <  2:658525/150‑1 (MQ=255)
             ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGat                                                                                                                                       <  2:460666/150‑1 (MQ=255)
              cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAATTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGata                                                                                                                                      <  1:466186/150‑1 (MQ=255)
              cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAATTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGata                                                                                                                                      >  2:466186/1‑150 (MQ=255)
                  ggTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAACGCCGAACTATCTCTATTAATAATTTAATGTTATATCTGCCACGATAAAACGGGGCAGATAatat                                                                                                                                  >  2:1226992/1‑150 (MQ=255)
                    tAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCATGGGCGTAGCGCTCTTCTGTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgt                                                                                                                                >  1:721158/1‑150 (MQ=255)
                    tAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATAATCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGatataacat                                                                                                                                <  1:1345235/150‑7 (MQ=255)
                      aaGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgttt                                                                                                                              <  2:1167830/150‑1 (MQ=255)
                        gTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttag                                                                                                                            <  1:714092/150‑1 (MQ=255)
                         tATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagt                                                                                                                           >  2:495484/1‑150 (MQ=255)
                          aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagtt                                                                                                                          <  2:25182/150‑1 (MQ=255)
                          aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCACGATAAAACGTTGCAGATAATAAgtttagta                                                                                                                          >  2:623409/1‑149 (MQ=255)
                            ttACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagttta                                                                                                                        >  2:736391/1‑150 (MQ=255)
                            ttACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagttta                                                                                                                        <  1:195205/150‑1 (MQ=255)
                              aCGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACt                                                                                                                      <  2:1350715/150‑1 (MQ=255)
                               cGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAACAAGGAGACCCGCCGAACTATCGCTTTTTATAATTTATTGTTATAAATGCCCCTATAAAACGGGGAAGATAATATGTTTAGTTTAATa                                                                                                                     >  1:541361/1‑150 (MQ=255)
                                 cTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAACCGGGGCAGATAATATGTTTAGTTTACTAAc                                                                                                                   >  2:272453/1‑150 (MQ=255)
                                  ttAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACg                                                                                                                  <  2:1673882/150‑1 (MQ=255)
                                   tAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACgg                                                                                                                 <  2:263803/150‑1 (MQ=255)
                                     gagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTc                                                                                                               >  2:1546738/1‑150 (MQ=255)
                                     gagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTc                                                                                                               <  1:1072574/150‑1 (MQ=255)
                                       ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                             >  1:173888/1‑150 (MQ=255)
                                        gCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATTATATGTTTAGTTTACTAACGGTCAtt                                                                                                            >  1:273342/1‑150 (MQ=255)
                                             aCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCa                                                                                                       <  2:1282087/150‑1 (MQ=255)
                                              ctctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAg                                                                                                      >  1:1181923/1‑150 (MQ=255)
                                               tctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGt                                                                                                     >  1:1275364/1‑150 (MQ=255)
                                               tctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGt                                                                                                     >  1:812420/1‑150 (MQ=255)
                                                ctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTg                                                                                                    <  2:812420/150‑1 (MQ=255)
                                                  aaaTATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGaa                                                                                                  >  2:1213494/1‑150 (MQ=255)
                                                  aaaTATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGaa                                                                                                  <  1:1213494/150‑1 (MQ=255)
                                                  aaaTATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGAAGATAATATGTTTAGTTTACTAACGTTCATTTTGCAGTGaa                                                                                                  >  1:154018/1‑150 (MQ=255)
                                                     tatCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGcc                                                                                               >  1:1003451/1‑150 (MQ=255)
                                                     tatCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGcc                                                                                               >  1:74576/1‑150 (MQ=255)
                                                     tatCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGcc                                                                                               >  2:1719554/1‑150 (MQ=255)
                                                       tCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAt                                                                                             >  2:1647041/1‑150 (MQ=255)
                                                        cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAATCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                            >  2:113059/1‑150 (MQ=255)
                                                        cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                            >  1:462376/1‑150 (MQ=255)
                                                          atatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTa                                                                                          >  2:259349/1‑150 (MQ=255)
                                                             ttCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTg                                                                                       <  2:356992/150‑1 (MQ=255)
                                                               ccTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtt                                                                                     >  2:1450107/1‑150 (MQ=255)
                                                                cTGGGCGTAGCGCTCTTCTTTTAGGGAAGACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttt                                                                                    <  2:63523/150‑1 (MQ=255)
                                                                cTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttt                                                                                    >  1:371003/1‑150 (MQ=255)
                                                                   ggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttttt                                                                                 <  2:1411804/150‑1 (MQ=255)
                                                                      gTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATc                                                                              <  1:733275/150‑1 (MQ=255)
                                                                          cgcTCTTCTTTTAGGGTAAACAGGAGAACCGCCGAACTAGCTCTATTTATAATGTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATGTAGTGTTTTTTATCGAcc                                                                          >  2:355920/1‑150 (MQ=255)
                                                                           gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                         <  2:789625/150‑1 (MQ=255)
                                                                           gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                         <  2:1541181/150‑1 (MQ=255)
                                                                            ctctTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAACCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCCTTTACTGTTTTTTATCGACCAg                                                                        >  1:1494677/1‑150 (MQ=255)
                                                                              cttcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGAt                                                                      <  2:1567416/150‑1 (MQ=255)
                                                                                tcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCGGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATaa                                                                    >  2:1087808/1‑150 (MQ=255)
                                                                                  ttttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATc                                                                  <  2:1717758/150‑1 (MQ=255)
                                                                                  ttttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGGCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATc                                                                  >  2:852722/1‑150 (MQ=255)
                                                                                     tAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGt                                                               <  1:1100018/150‑1 (MQ=255)
                                                                                      aGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                              >  1:401169/1‑150 (MQ=255)
                                                                                       gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                             >  2:1147459/1‑150 (MQ=255)
                                                                                         gAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtctc                                                           >  2:1334349/1‑150 (MQ=255)
                                                                                          aaaaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTctct                                                          <  2:661569/150‑1 (MQ=255)
                                                                                           aaaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTCTTATCGACCAGATAATCTGTTCTCTa                                                         <  2:517289/150‑1 (MQ=255)
                                                                                           aaaCAGGAGAACCGCCGAACTATCTCTATTTATAATATAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGttatata                                                         >  1:1035269/1‑145 (MQ=255)
                                                                                              cAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGATCTCTAAtt                                                      >  2:509032/1‑149 (MQ=255)
                                                                                                ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGtt                                                    >  1:177026/1‑150 (MQ=255)
                                                                                                ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGtt                                                    <  1:259349/150‑1 (MQ=255)
                                                                                                 gagaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                   >  2:932833/1‑150 (MQ=255)
                                                                                                   gaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAAc                                                 <  1:1719554/150‑1 (MQ=255)
                                                                                                     aCCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGGGAAGCCATTTACTGTTTTTTATCGACCAGAGAATCTGTTCTCTTATGTTAACta                                               >  2:890065/1‑149 (MQ=255)
                                                                                                     aCCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGAAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTc                                               >  2:1128681/1‑150 (MQ=255)
                                                                                                      ccgccgAACTATCTCTATTTATTATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTcc                                              >  2:949140/1‑150 (MQ=255)
                                                                                                      ccgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACCGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTcc                                              >  1:954903/1‑150 (MQ=255)
                                                                                                       cgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                             >  1:935841/1‑150 (MQ=255)
                                                                                                         cctaaCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccccc                                           <  2:401169/147‑1 (MQ=255)
                                                                                                             aCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAc                                       >  2:439243/1‑150 (MQ=255)
                                                                                                              cTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                      <  2:173888/150‑1 (MQ=255)
                                                                                                              atATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                      <  1:479080/149‑1 (MQ=255)
                                                                                                                 tctcTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGt                                   <  1:1396831/150‑1 (MQ=255)
                                                                                                                         tttaAATTAATGTTATATCTGCACCGATAAAACGGGGCAGAAAATATGTTTAGTTTACTAAAGGTCATTTTGCATTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                           <  1:272453/148‑1 (MQ=255)
                                                                                                                         tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                           <  2:1162134/150‑1 (MQ=255)
                                                                                                                         tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                           >  2:1576495/1‑150 (MQ=255)
                                                                                                                                tAATGTTATTTCTGACCCGATAAATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGTCCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCa                    <  1:1582406/150‑1 (MQ=255)
                                                                                                                                  aTGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAtt                  >  2:908645/1‑150 (MQ=255)
                                                                                                                                  aTGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAtt                  <  1:1546738/150‑1 (MQ=255)
                                                                                                                                    gTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTc                <  1:932833/150‑1 (MQ=255)
                                                                                                                                      tataTCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCt              >  1:271762/1‑150 (MQ=255)
                                                                                                                                       atatCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg             <  2:736734/150‑1 (MQ=255)
                                                                                                                                        tatCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGt            >  1:674011/1‑150 (MQ=255)
                                                                                                                                        tatCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCGATCCTGTTGCTTTAGTTATTCTTTTCCTGt            >  1:736708/1‑150 (MQ=255)
                                                                                                                                         atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc           <  1:64032/150‑1 (MQ=255)
                                                                                                                                          tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct          >  2:1302830/1‑150 (MQ=255)
                                                                                                                                               cccGATAAAACGGGGCGGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACtt     <  1:1087808/150‑1 (MQ=255)
                                                                                                                                               cccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACtt     >  1:696685/1‑150 (MQ=255)
                                                                                                                                                  gATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGc  <  2:1594130/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGC  >  NC_000913/2788524‑2788819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: