Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,208 T→A 5.5% intergenic (‑30/+134) rbfA ← / ← infB 30s ribosome binding factor/translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,2080TA 5.5% 422.3 / 7.1 142intergenic (‑30/+134)rbfA/infB30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base T (64/70);  new base A (4/4);  total (68/74)
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.19e-01

Only 100 of 145 total aligned reads displayed.

ACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAC  >  NC_000913/3313059‑3313355
                                                                                                                                                      |                                                                                                                                                   
 ccGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                    <  2:737736/150‑1 (MQ=255)
 ccGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                    <  2:1080834/150‑1 (MQ=255)
  cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTa                                                                                                                                                   >  2:1598702/1‑150 (MQ=255)
  cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTa                                                                                                                                                   >  2:475068/1‑150 (MQ=255)
  cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATACTTTCTCCATAATAAATTCTCCTGACAAAAAACGTGCTGTTa                                                                                                                                                   >  2:1275259/1‑150 (MQ=255)
  cGGAAAAGGTGGTCATCATGCCCAGGCGATGATCTTTAATTTCACGCTGCAGGATGAGAGCGATATC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTa                                                                                                                                                   <  1:301646/150‑1 (MQ=255)
     aaaCGGTGGTCATCATGCGCAGGCGAGTATCTTTAATTTCACGCTGCAGGATGAGAGCGATATC‑TTTTTGCATTTCCTGCGCTACTCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGcc                                                                                                                                                <  1:133151/150‑1 (MQ=255)
     aaaCGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGcc                                                                                                                                                <  1:543396/150‑1 (MQ=255)
     aaaCGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGcc                                                                                                                                                >  2:18188/1‑150 (MQ=255)
          gtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCTAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                           <  1:895602/150‑1 (MQ=255)
           tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtttt                                                                                                                                          <  2:1039608/150‑1 (MQ=255)
           tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtttt                                                                                                                                          <  1:1630342/150‑1 (MQ=255)
           tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtttt                                                                                                                                          <  2:1379991/150‑1 (MQ=255)
           tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt                                                                                                                                          <  2:1404079/150‑1 (MQ=255)
            ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                         <  2:1430070/150‑1 (MQ=255)
            ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                         <  1:1486215/150‑1 (MQ=255)
              tcatcaTGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaa                                                                                                                                       >  1:1015878/1‑150 (MQ=255)
               catcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaa                                                                                                                                      >  2:1539525/1‑150 (MQ=255)
                atcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                     >  1:1328406/1‑150 (MQ=255)
                atcatGCCAAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGACATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                     <  1:448578/150‑1 (MQ=255)
                 tcatTCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaat                                                                                                                                    <  2:320074/150‑1 (MQ=255)
                 tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCTCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaat                                                                                                                                    <  1:309151/150‑1 (MQ=255)
                  catGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaatt                                                                                                                                   >  2:41613/1‑150 (MQ=255)
                   ttGCCCAGGCGAGTATCTTTAATTTCACGCGGCCGGCTGAGAGAGATCTCTTTTTTGCATTTCCTTCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaatt                                                                                                                                   <  1:728693/149‑1 (MQ=255)
                   atGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaatta                                                                                                                                  >  1:1293782/1‑150 (MQ=255)
                     tcccAGGCGATGATCTTTAATTTCACGCTTCATGATGTTAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGTCCAAATTATTTCGCCATAATAAATTCTACTGACAAAAAATGGGCTGTTAGGCCCTTTTTAAaattaat                                                                                                                                <  2:1002666/149‑1 (MQ=255)
                      cccACGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGTTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaatt                                                                                                                               <  2:1607832/150‑1 (MQ=255)
                          ggCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAg                                                                                                                           >  2:702892/1‑150 (MQ=255)
                            cGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGt                                                                                                                         <  2:1515307/150‑1 (MQ=255)
                                  tCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAgg                                                                                                                   >  1:1531/1‑150 (MQ=255)
                                    tttAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAGTTCTCATGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGATGGGc                                                                                                                 >  2:1523032/1‑150 (MQ=255)
                                     ttAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCt                                                                                                                <  2:467588/150‑1 (MQ=255)
                                        aTTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATATTAAATTCTCCTGACAAAAAAGGTGCTGTTAGCCACTTTTTAATATTAATTTCAGGTGGATGGGCTGtt                                                                                                             >  1:1100418/1‑150 (MQ=255)
                                         tttCCCGCTGCAGGATGAGAGCGAGCTC‑TTTTTGCATTTCCTGCGCGACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTc                                                                                                            <  2:1152919/150‑1 (MQ=255)
                                         tttCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTc                                                                                                            <  1:1677331/150‑1 (MQ=255)
                                            cACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAACTTCTTTCGCCATAATAAATTCTCCAGACAAAAAAGGCGCTCTTAGCCCCTTTTTAAAATTAATTTCAGGTGCAAGGGCTGTTCACg                                                                                                         >  2:68963/1‑150 (MQ=255)
                                            cACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACg                                                                                                         <  2:1650430/150‑1 (MQ=255)
                                             aCGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGt                                                                                                        <  2:966991/150‑1 (MQ=255)
                                               gcggccggATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTg                                                                                                      <  2:940404/144‑1 (MQ=255)
                                                cTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGa                                                                                                     >  1:319397/1‑150 (MQ=255)
                                                cTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGa                                                                                                     >  2:1355519/1‑150 (MQ=255)
                                                 tGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGAc                                                                                                    <  2:913173/150‑1 (MQ=255)
                                                 tGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGAc                                                                                                    <  1:90989/150‑1 (MQ=255)
                                                      gATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGAt                                                                                               <  1:1436903/150‑1 (MQ=255)
                                                      gATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGAt                                                                                               <  1:421019/150‑1 (MQ=255)
                                                         gagaGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAg                                                                                            >  2:513316/1‑150 (MQ=255)
                                                          agagCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                           >  2:215962/1‑150 (MQ=255)
                                                          agagCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGa                                                                                           >  1:1638276/1‑150 (MQ=255)
                                                               gATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                      <  1:846550/150‑1 (MQ=255)
                                                               gATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                      >  1:359320/1‑150 (MQ=255)
                                                               gATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                      <  2:1061897/150‑1 (MQ=255)
                                                                aTCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGcc                                                                                     <  2:1293782/150‑1 (MQ=255)
                                                                aTCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGcc                                                                                     <  2:1057818/150‑1 (MQ=255)
                                                                 tctc‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCa                                                                                    <  1:747933/150‑1 (MQ=255)
                                                                  ctc‑tTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAg                                                                                   <  1:117211/150‑1 (MQ=255)
                                                                        tttGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGATGCGCCAGCGTca                                                                              >  1:1454400/1‑150 (MQ=255)
                                                                           gCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACAt                                                                           >  1:1399210/1‑150 (MQ=255)
                                                                              tttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAg                                                                        <  1:1539525/150‑1 (MQ=255)
                                                                                  cTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCaa                                                                    >  1:1680233/1‑150 (MQ=255)
                                                                                     cgcTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATcc                                                                 <  2:411589/150‑1 (MQ=255)
                                                                                      gcTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGTGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCa                                                                <  2:694716/150‑1 (MQ=255)
                                                                                       cTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCAt                                                               <  2:1712839/150‑1 (MQ=255)
                                                                                         aCGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAAATCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTGCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGc                                                             >  1:1521939/1‑150 (MQ=255)
                                                                                               tGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCAATGCCGGATg                                                       >  1:995602/1‑150 (MQ=255)
                                                                                                       aCCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaa                                               >  1:1724035/1‑150 (MQ=255)
                                                                                                        ccAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGCCCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaaa                                              >  1:1283835/1‑150 (MQ=255)
                                                                                                        ccAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaaa                                              >  1:1443194/1‑150 (MQ=255)
                                                                                                          aaaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACg                                            >  2:1029005/1‑150 (MQ=255)
                                                                                                          aaaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACg                                            <  1:524195/150‑1 (MQ=255)
                                                                                                            aTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGcc                                          <  1:1603354/150‑1 (MQ=255)
                                                                                                             ttCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCt                                         <  1:1355519/150‑1 (MQ=255)
                                                                                                              tCTTTAGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCtt                                        <  1:1049288/150‑1 (MQ=255)
                                                                                                                   cGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATccc                                   <  1:407794/150‑1 (MQ=255)
                                                                                                                    gCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCg                                  >  1:87505/1‑150 (MQ=255)
                                                                                                                    gCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCg                                  >  1:1433349/1‑150 (MQ=255)
                                                                                                                     ccATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGc                                 >  1:1499077/1‑150 (MQ=255)
                                                                                                                         aataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGTATGCAGCGTAAACGCCTTATCCCGCATgg                             >  1:185191/1‑150 (MQ=255)
                                                                                                                            aacttCTCCTTACAAACAAGTGGCGGTTAGCCGCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAc                          <  2:254144/147‑1 (MQ=38)
                                                                                                                            aaaTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAc                          >  2:893809/1‑150 (MQ=255)
                                                                                                                            aaaTTCTCCTGACAAAAAAGGGGCTGATAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAc                          >  1:1478815/1‑150 (MQ=255)
                                                                                                                             aaTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATAAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCGGGATAAGGCGTTTACGCTGCATCCGGCATGGattg                         <  1:1448799/150‑4 (MQ=37)
                                                                                                                              atactcCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCTCGTTGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGAGTCACATCAGGCAATCCATGAGGGATAAGGCGTTTACGCTGCATCCGGCATGGattgc                        <  1:289599/147‑5 (MQ=25)
                                                                                                                              aTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAcc                         >  2:1539344/1‑149 (MQ=255)
                                                                                                                              aTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAcc                         <  1:1539344/149‑1 (MQ=255)
                                                                                                                                 ctcCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTaa                     >  2:1662348/1‑150 (MQ=255)
                                                                                                                                      gACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAcc                >  1:1014556/1‑150 (MQ=255)
                                                                                                                                       tttaaaaaGGGGCTAACAGCCCGTTTTTAATATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTTAACGCCTTATCCCGCATGGAACACTAAAAAcat               >  1:135164/4‑148 (MQ=38)
                                                                                                                                        cAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCtt              <  2:32483/150‑1 (MQ=255)
                                                                                                                                        cAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCtt              >  2:180458/1‑150 (MQ=255)
                                                                                                                                         aaaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTa             <  1:215962/150‑1 (MQ=255)
                                                                                                                                          aaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa            >  2:117389/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa            >  2:244071/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa            >  2:1072936/1‑150 (MQ=255)
                                                                                                                                           aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg           <  1:1598702/150‑1 (MQ=255)
                                                                                                                                            aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGc          <  2:1627504/150‑1 (MQ=255)
                                                                                                                                            aaaGGGGCTAAAAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGc          <  1:445318/150‑1 (MQ=255)
                                                                                                                                                gggCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATg      >  1:945398/1‑150 (MQ=255)
                                                                                                                                                 ggCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATgg     >  1:1350860/1‑150 (MQ=255)
                                                                                                                                                   cTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTa   >  2:1227180/1‑150 (MQ=255)
                                                                                                                                                    tGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGACGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAc  <  2:1515948/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTC‑TTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAC  >  NC_000913/3313059‑3313355

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