Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,537,324 T→A 10.5% intergenic (+36/‑61) greB → / → yhgF transcript cleavage factor/putative transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,3240TA10.5% 384.2 / 26.4 143intergenic (+36/‑61)greB/yhgFtranscript cleavage factor/putative transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base T (60/68);  new base A (8/7);  total (68/75)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 149 total aligned reads displayed.

CCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGAC  >  NC_000913/3537176‑3537471
                                                                                                                                                    |                                                                                                                                                   
ccccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAtttt                                                                                                                                                    <  2:152640/150‑1 (MQ=255)
 cccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAttttt                                                                                                                                                   >  1:265289/1‑150 (MQ=255)
 cccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAttttt                                                                                                                                                   >  1:1405731/1‑150 (MQ=255)
 cccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAttttt                                                                                                                                                   >  1:1328812/1‑150 (MQ=255)
    gATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgcc                                                                                                                                                <  1:811570/150‑1 (MQ=255)
      tGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccag                                                                                                                                              <  2:963454/150‑1 (MQ=255)
        gCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccagcc                                                                                                                                            >  1:431017/1‑150 (MQ=255)
        gCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccagcc                                                                                                                                            >  1:1685262/1‑150 (MQ=255)
        gCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGGCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagcc                                                                                                                                            <  2:1247583/150‑1 (MQ=255)
           cgcgCATTGCTGAAAAAAGATGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatc                                                                                                                                         <  2:697613/150‑3 (MQ=255)
              gcATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatcggc                                                                                                                                      <  2:758358/150‑6 (MQ=255)
                    cTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAAc                                                                                                                                <  1:76775/150‑1 (MQ=255)
                    cTGAAAAAAGAAGTCGGCGATCCGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAAc                                                                                                                                <  1:845457/150‑1 (MQ=255)
                     tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACt                                                                                                                               >  2:1522218/1‑150 (MQ=255)
                      gAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCAGTATAACTa                                                                                                                              >  2:1239336/1‑150 (MQ=255)
                       aaaaaaGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCCGCCAGTCCGTATAACTAt                                                                                                                             >  1:987909/1‑150 (MQ=255)
                         aaaaGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGGGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATcc                                                                                                                           >  1:86629/1‑150 (MQ=255)
                                  cGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGATGGCATTTTTTCCATCCAGTCCGTATAAATATCCCCTGAtttt                                                                                                                  >  2:1247348/1‑150 (MQ=255)
                                     cGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGa                                                                                                               <  1:1522218/150‑1 (MQ=255)
                                      gATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTTATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGAt                                                                                                              >  2:291842/1‑150 (MQ=255)
                                      gATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTTATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGAt                                                                                                              >  2:944237/1‑150 (MQ=255)
                                        tCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATcc                                                                                                            <  2:1027106/150‑1 (MQ=255)
                                         cTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCg                                                                                                           >  2:1137842/1‑150 (MQ=255)
                                          tGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGa                                                                                                          >  2:78174/1‑150 (MQ=255)
                                           ggCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGaa                                                                                                         >  2:31910/1‑150 (MQ=255)
                                             cGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatcg                                                                                                                                        >  1:105891/1‑114 (MQ=255)
                                             cGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAAAAATgccagccatcg                                                                                                                                        <  2:105891/117‑4 (MQ=255)
                                              ggtggtGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAg                                                                                                      <  2:1640029/150‑1 (MQ=255)
                                               gtggtgAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGa                                                                                                     >  2:1197803/1‑150 (MQ=255)
                                                    gAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAAc                                                                                                >  2:714242/1‑150 (MQ=255)
                                                     aaTACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACt                                                                                               >  2:1017100/1‑150 (MQ=255)
                                                     aaTACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACt                                                                                               <  2:1462899/150‑1 (MQ=255)
                                                           ccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAcc                                                                                         <  1:463528/150‑1 (MQ=255)
                                                                 ggggAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatgatg                                                                                   <  2:405272/150‑1 (MQ=255)
                                                                  gggAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCAtgatga                                                                                  <  1:1022533/150‑1 (MQ=255)
                                                                     aaGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCGAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGatgaatg                                                                               <  2:1561621/150‑1 (MQ=255)
                                                                      aGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGatgaatga                                                                              >  1:515344/1‑150 (MQ=255)
                                                                        cGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTTATCCGAAAAGATGAACTCAAACCATGATGAATGAtt                                                                            >  1:865533/1‑150 (MQ=255)
                                                                          aGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCg                                                                          <  2:476679/150‑1 (MQ=255)
                                                                          aGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCg                                                                          <  2:1104082/150‑1 (MQ=255)
                                                                           gCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGt                                                                         >  1:1215253/1‑150 (MQ=255)
                                                                             tGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTc                                                                       <  2:1328812/150‑1 (MQ=255)
                                                                              gttATGTTAATTCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCt                                                                      <  1:778372/148‑1 (MQ=255)
                                                                                 aTGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGcc                                                                   <  2:520757/150‑1 (MQ=255)
                                                                                        tGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCattat                                                            >  2:135708/1‑150 (MQ=255)
                                                                                          cTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGc                                                                 >  1:560321/1‑143 (MQ=255)
                                                                                          cTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGc                                                                 <  2:560321/143‑1 (MQ=255)
                                                                                           tATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatga                                                                                     <  2:599537/122‑1 (MQ=255)
                                                                                           tATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatga                                                                                     >  1:599537/1‑122 (MQ=255)
                                                                                            aTCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAGGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCg                                                        <  1:703325/150‑1 (MQ=255)
                                                                                            aTCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCg                                                        <  1:959327/150‑1 (MQ=255)
                                                                                             tCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCgg                                                       <  1:841909/150‑1 (MQ=255)
                                                                                                aGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTg                                                    <  2:618305/150‑1 (MQ=255)
                                                                                                    ggTGAAACCGTTAGGAAGATTCATAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTTTCGCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                                <  2:231836/149‑1 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAAt                                                >  1:1032235/1‑150 (MQ=255)
                                                                                                    cGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTACGTATAACTATCCCAAGATTTTTTAAACGAAAAGATGAACTCAAACCATGATGAAAGATTACTTCTGCCCCATTATTGCGGGTGAAAt                                                >  2:735503/1‑150 (MQ=255)
                                                                                                      tGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTTTCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTTAATTTc                                              >  2:907488/1‑150 (MQ=255)
                                                                                                      tGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGAATCGTTCTGCCGCATTATTGCGGGTGAAAtta                                              >  1:1057642/1‑149 (MQ=255)
                                                                                                       gAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCa                                             >  2:646311/1‑150 (MQ=255)
                                                                                                       gAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCa                                             >  1:16723/1‑150 (MQ=255)
                                                                                                       gAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCa                                             <  2:265289/150‑1 (MQ=255)
                                                                                                         aaCCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAgg                                           <  2:1298150/150‑1 (MQ=255)
                                                                                                         aaCCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAgg                                           >  2:689578/1‑150 (MQ=255)
                                                                                                         aaCCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAgg                                           <  2:431017/150‑1 (MQ=255)
                                                                                                           ccGTAAGTAAGAGTCTTAACCTCCTGCCGATGGATGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcg                                         <  1:1238162/150‑1 (MQ=255)
                                                                                                           ccGTAAGGAAGAGTCTTAACCTCCTTCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcg                                         <  1:1344853/150‑1 (MQ=255)
                                                                                                           ccGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATATCTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcg                                         <  2:233289/150‑1 (MQ=255)
                                                                                                             gTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAACGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcgcg                                       <  2:987946/150‑1 (MQ=255)
                                                                                                              tAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcgcgc                                      <  1:291842/150‑1 (MQ=255)
                                                                                                              tAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGgcgcgc                                      <  1:944237/150‑1 (MQ=255)
                                                                                                                    agagTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGTTTTTTGATCCGATAAGATGTACTCAAACCATGATGATTGATTCGTTCTGCCGCATTATTGCGGGTGAAATTGAGGCGCGCCCGGat                                <  1:817760/150‑2 (MQ=255)
                                                                                                                      agTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACa                              >  1:1231353/1‑150 (MQ=255)
                                                                                                                       gTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAg                             <  1:649616/150‑1 (MQ=255)
                                                                                                                       gTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCAGCATTATTGCTGGTGAAATTCAGGCGCGCCCGGAACAg                             >  2:1142650/1‑150 (MQ=255)
                                                                                                                         cTTAACCTCCTGCCGATTGCTGGCATTTTTGCCAGCCAGTCCGTATAACAATCCCCTGATTTTTGATACGAAAAGATGAACTCAAACCCTGATGAATGATTCGTTCATCCGCATTATTGCGGGTTAAATTCAGGCGCGCCCGGAACAGGt                           >  1:1489123/1‑150 (MQ=255)
                                                                                                                         cTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTTATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGt                           >  1:481967/1‑150 (MQ=255)
                                                                                                                         cTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGt                           <  1:1017100/150‑1 (MQ=255)
                                                                                                                           tAACCTCCTGCCTATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTg                         <  2:1484302/150‑1 (MQ=255)
                                                                                                                             aCCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGAc                       >  2:766319/1‑150 (MQ=255)
                                                                                                                              cctcctGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACg                      >  2:397261/1‑150 (MQ=255)
                                                                                                                                  ctGCCGATGGCTGGCATTTTTGCCAGCCAGTCCATATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGc                  >  2:584744/1‑150 (MQ=255)
                                                                                                                                   tGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGGATAACTATCCCCGGAGTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGGTCTGCCGCAGTATTGAGGGTGCAATTCAGTCGCGCACGGAACCGGTTGACGCTGcc                 >  1:340955/1‑150 (MQ=255)
                                                                                                                                    gCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGCTCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCg                >  2:436474/1‑150 (MQ=255)
                                                                                                                                    gCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCg                >  2:1234308/1‑150 (MQ=255)
                                                                                                                                      cGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGtt              <  2:865533/150‑1 (MQ=255)
                                                                                                                                      cGATGGCGGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGtt              <  2:1405731/150‑1 (MQ=255)
                                                                                                                                       actggctggcAAAAATGCCAGCCAGTCCGTATACCGATCCCCTGATTTTTGATCCGAAAAGATGAACTCCAACCATGATGAATGAGTCGTTCTGCCGCAGTATTGCGGGTGAAATTCCGGCGCGCCCGGAACAGGTTGACGCTGCCGTTc             >  1:1584103/3‑150 (MQ=255)
                                                                                                                                       actggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGATCTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTc             >  1:1577689/3‑150 (MQ=255)
                                                                                                                                       actggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTc             >  2:1417319/3‑150 (MQ=255)
                                                                                                                                        ctggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCg            >  2:136448/2‑150 (MQ=255)
                                                                                                                                        ctggctggcAAAAATGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTAGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGAGCCCGGAACAGGTTGACGCTGCCGTTCg            >  2:1677586/2‑150 (MQ=255)
                                                                                                                                         tggctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGc           <  1:135708/150‑1 (MQ=255)
                                                                                                                                           gctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGTATGATTCGTTCTCCCGCATTATTGCGGGTGAAATTCAGGCGCGCGCGGAACAGGTTGACGCTGCCGTTCGccg         >  2:1546111/1‑149 (MQ=255)
                                                                                                                                           gctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGa                                                   >  2:297339/1‑108 (MQ=255)
                                                                                                                                           gctggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGa                                                   <  1:297339/108‑1 (MQ=255)
                                                                                                                                              gtcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCt      <  1:982438/148‑1 (MQ=255)
                                                                                                                                              ggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGACCTCAAACCCTGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCt      <  1:1197803/150‑1 (MQ=255)
                                                                                                                                              ggcATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCt      >  2:65397/1‑150 (MQ=255)
                                                                                                                                                cATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTg    <  2:515344/150‑1 (MQ=255)
                                                                                                                                                cATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTg    >  2:476788/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGGTGAAATTCAGGCGCGCCCGGAACAGGTTGACGCTGCCGTTCGCCTGCTTGAC  >  NC_000913/3537176‑3537471

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