Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,966 T→A 5.4% intergenic (‑39/+18) ralR ← / ← recT Rac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,9660TA 5.4% 312.9 / 4.1 103intergenic (‑39/+18)ralR/recTRac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (40/57);  new base A (3/3);  total (43/60)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 103 total aligned reads displayed.

TTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATAT  >  NC_000913/1413819‑1414113
                                                                                                                                                   |                                                                                                                                                   
ttCGGGATTCGGATCCGTTACACTTCGCTCGGTAATTTCCTGATATGGCTTTCACCGTTACTGATTGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                                                                   <  1:865264/150‑1 (MQ=255)
ttCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                                                                   <  2:1627888/150‑1 (MQ=255)
 tCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTac                                                                                                                                                  <  2:708903/150‑1 (MQ=255)
  cGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTCCTGATGGACTACCACAAAATGGAGATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGAGGTGCATTacc                                                                                                                                                 >  2:265259/1‑149 (MQ=255)
    ggATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAcc                                                                                                                                               <  1:543434/150‑1 (MQ=255)
      aTTCGCATCCGTAACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGTCATGGTTTAACATTGTCATATCTCATAACTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGc                                                                                                                                             <  1:385122/150‑1 (MQ=255)
         cGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                                                                          <  1:957702/150‑1 (MQ=255)
         cGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                                                                          <  2:801955/150‑1 (MQ=255)
         cGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTCTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                                                                          <  1:44918/150‑1 (MQ=255)
           cATCCGTTACACTTCGCTCGGTAATATCCTGAAATGTCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                                                                        <  2:1378960/150‑4 (MQ=255)
           cATCCGTTACACTTCGCTCGGTAATATCCTGAAATCGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAGCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                                                                        <  2:304639/150‑4 (MQ=255)
             tCCGTTACACTTCGCTCTGTAATATCATGAAATGTCTTTAACCGTTACTGATGGACAACCACAAAATTGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTg                                                                                                                                      <  1:1165353/150‑1 (MQ=255)
             tCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCACATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTg                                                                                                                                      >  1:311486/1‑150 (MQ=255)
              ccGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGa                                                                                                                                     >  2:913936/1‑150 (MQ=255)
               cGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                                                                    >  2:1010195/1‑150 (MQ=255)
               cGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                                                                    >  2:855600/1‑150 (MQ=255)
                gTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                                                                   <  1:1218175/150‑1 (MQ=255)
                   acacTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattat                                                                                                                                <  2:569800/150‑1 (MQ=255)
                    cacTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                                                               >  1:1109749/1‑150 (MQ=255)
                    cacTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                                                               >  1:784715/1‑150 (MQ=255)
                     acTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTAGTc                                                                                                                              >  2:1577257/1‑150 (MQ=255)
                          gCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                                                                         <  2:1034381/150‑1 (MQ=255)
                          gCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTAAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                                                                         <  1:700433/150‑1 (MQ=255)
                             cGGTAATATCCTGAATTGGCTTTCACCGATACTGATGGACAACCACAAAATGGAAATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTTCTTTACACCGCAAGGCTGAATTATTCCTCTGAAt                                                                                                                      >  1:782838/1‑150 (MQ=255)
                             cGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAt                                                                                                                      >  2:664249/1‑150 (MQ=255)
                             cGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAt                                                                                                                      >  1:665742/1‑150 (MQ=255)
                              ggTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAtt                                                                                                                     <  2:912282/150‑1 (MQ=255)
                              ggTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAtt                                                                                                                     >  2:1071313/1‑150 (MQ=255)
                                  atatCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                                                                 <  1:1070488/150‑1 (MQ=255)
                                     tcatGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCACGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAtt                                                                                                              <  1:112587/147‑1 (MQ=255)
                                           aaTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                                                                        <  1:11514/150‑1 (MQ=255)
                                            aTGGCTTTCACCGTTACTGATTGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGt                                                                                                       <  1:1063746/150‑1 (MQ=255)
                                            aTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTTCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGt                                                                                                       <  1:1069303/150‑1 (MQ=255)
                                                   tCACCGTAACTGATGGACAACCAAAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcccc                                                                                                <  1:293257/150‑1 (MQ=255)
                                                    cACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCg                                                                                               <  2:843328/150‑1 (MQ=255)
                                                    cACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCg                                                                                               >  2:607332/1‑150 (MQ=255)
                                                       cGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                                                                            <  2:489562/150‑1 (MQ=255)
                                                             tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACa                                                                                      >  1:858345/1‑150 (MQ=255)
                                                               aTGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                                                                    <  1:1717168/150‑1 (MQ=255)
                                                               aTGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTAATCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                                                                    >  1:377735/1‑150 (MQ=255)
                                                                tGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                                                                   >  2:1675578/1‑150 (MQ=255)
                                                                 ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                                                                  <  1:789707/150‑1 (MQ=255)
                                                                 ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCAAGGTTAATACAGAgg                                                                                  >  2:1286568/1‑150 (MQ=255)
                                                                 ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATCTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCATCTGAATTATCGATTACACTGTATTCACCGGTTAATACAGAgg                                                                                  >  1:59573/1‑150 (MQ=255)
                                                                      aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                                                                             <  2:1326092/150‑1 (MQ=255)
                                                                        cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                                                                           >  2:450973/1‑150 (MQ=255)
                                                                        cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACAGACAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTACCAGGTTAATACAGAGGAATCTGc                                                                           >  1:865894/1‑150 (MQ=255)
                                                                         acaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCa                                                                          <  2:1259301/150‑1 (MQ=255)
                                                                          caAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAg                                                                         <  1:855600/150‑1 (MQ=255)
                                                                            aaaTGTACATGGTTTAACATTGTCAAATCTCATAATTTTTCTCATAAATAATATTTCTAGTTGGCGGTGCTTTACTCCGCCAGGCTTATTTATTCCTCTGAATTATCGATTCCTCTGTATTAACCGGGGAATACAGTGTAATcgataatt                                                                       <  1:277411/150‑8 (MQ=21)
                                                                               tgtaCATGGTTTAACATTGTCTTATCTAATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCg                                                                    <  2:311486/147‑1 (MQ=255)
                                                                                   cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                                                                <  1:1048967/150‑1 (MQ=255)
                                                                                   cATGGATTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                                                                <  2:1203618/150‑1 (MQ=255)
                                                                                     tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc                                                              <  1:469997/150‑1 (MQ=255)
                                                                                      ggTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCa                                                             >  1:2146/1‑150 (MQ=255)
                                                                                       gTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCATCTGAATTATCGATTACACTGTATTCCACGGTTAGTACAGAGGAATCTGCAGGATCGATTGTctg                                                            >  2:422074/1‑148 (MQ=255)
                                                                                          tAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCGGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAGCCAGAGCAATCTGCAGGATCGATTGTCAGtcg                                                         >  2:1336507/1‑148 (MQ=255)
                                                                                              aTTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTcc                                                     <  2:518950/150‑1 (MQ=255)
                                                                                                 gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTACTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt                                                  >  2:650375/1‑150 (MQ=255)
                                                                                                  tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTTCACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt                                                 <  2:695447/150‑1 (MQ=255)
                                                                                                  tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt                                                 <  2:1652031/150‑1 (MQ=255)
                                                                                                    ttatCTCATAATTTTTCTCATAAAAAATTTTTCAAGTTGGCGGTGCTTTACACCGCGAGGCTGAATTTTTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACATAGGAATCTGCAGGATCGATTGTCAGTGGTTGCTTTTCa                                               <  1:1341922/149‑1 (MQ=255)
                                                                                                    atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                               <  1:393354/150‑1 (MQ=255)
                                                                                                     tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTcc                                                     <  2:934155/143‑1 (MQ=255)
                                                                                                     tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTcc                                                     >  1:934155/1‑143 (MQ=255)
                                                                                                     tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCAt                                              <  1:1259419/150‑1 (MQ=255)
                                                                                                     tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCAt                                              <  1:1586717/150‑1 (MQ=255)
                                                                                                         tcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCAt                                          >  1:1082883/1‑150 (MQ=255)
                                                                                                          cATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCAtt                                         <  2:1731132/150‑1 (MQ=255)
                                                                                                            tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGa                                       >  2:1705527/1‑150 (MQ=255)
                                                                                                              aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATa                                     <  1:1078577/150‑1 (MQ=255)
                                                                                                              aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATa                                     <  1:913936/150‑1 (MQ=255)
                                                                                                               tttttCTCATAAAAAATATTTCAAGTTGGCTGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATAc                                    <  1:635953/150‑1 (MQ=255)
                                                                                                               tttttCTCATAAAAAATATTTCAAGTTGGCGGTTCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATAc                                    <  1:412708/150‑1 (MQ=255)
                                                                                                                 tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTg                                  <  1:167488/150‑1 (MQ=255)
                                                                                                                   tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCa                                <  2:698586/150‑1 (MQ=255)
                                                                                                                   tatcATAAAAAATATTTCAAGTTGGCTGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCa                                <  2:132036/148‑1 (MQ=255)
                                                                                                                      aaTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGc                             <  2:1102943/149‑1 (MQ=255)
                                                                                                                       aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGAACGCt                            >  1:717160/1‑150 (MQ=255)
                                                                                                                         aaaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTgg                          <  1:1071313/150‑1 (MQ=255)
                                                                                                                         aaaaaaTATTTCAAGTTGGCGGTGCAATACACCGCCAGGCTGAATTATTCCTCTGAAATATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTCGTACCTTATCATCCATTCATACTGCACTCTga                          >  1:1292183/1‑149 (MQ=255)
                                                                                                                           aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAt                        >  2:1269026/1‑150 (MQ=255)
                                                                                                                           aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAt                        <  2:1683111/150‑1 (MQ=255)
                                                                                                                              atatTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAtctc                     <  2:784715/150‑1 (MQ=255)
                                                                                                                                   tCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTg                <  1:739770/150‑1 (MQ=255)
                                                                                                                                    gatcttGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTGCCCGGTTAATACTGAGGAATCTGCAGGATCGATTGTCTGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTgt               <  1:1089513/146‑2 (MQ=255)
                                                                                                                                      gcctGGCGGTGTAATGCACCGCCAGGATGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATa             >  1:1392259/4‑150 (MQ=255)
                                                                                                                                       cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATAc            >  2:1653051/3‑150 (MQ=255)
                                                                                                                                       cctGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATAGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATAc            >  2:91501/3‑150 (MQ=255)
                                                                                                                                         ttgagttGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATGGATTACACTGTATTCCGCGGTTAATTCTGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCTGGCTGGATCTCAATTGATACgg          <  2:782838/144‑1 (MQ=255)
                                                                                                                                         tGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGa               >  2:1089514/1‑145 (MQ=255)
                                                                                                                                         tGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGa               <  1:1089514/145‑1 (MQ=255)
                                                                                                                                           gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGc        >  2:1277811/1‑150 (MQ=255)
                                                                                                                                              gTGCATTACACCGCCAGGCTGCATTATTCCTCTGAATTTTCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCAGCAGGATCGATTGTCAGTGGTTCCTGTTCATCCATTGATACTGCACGCTGGCTCTCAATTGATCCGGGCaaa     >  2:475814/1‑150 (MQ=255)
                                                                                                                                              gTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGCTCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCaaa     >  1:1500236/1‑150 (MQ=255)
                                                                                                                                               tGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGATTCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAat    >  2:911156/1‑150 (MQ=255)
                                                                                                                                               tGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAat    >  1:775249/1‑150 (MQ=255)
                                                                                                                                                gCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAata   >  2:1331592/1‑150 (MQ=255)
                                                                                                                                                 cATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAatat  >  1:1566616/1‑150 (MQ=255)
                                                                                                                                                 cATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAatat  >  1:1732346/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATAT  >  NC_000913/1413819‑1414113

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: