Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,940 T→A 7.6% F61I (TTC→ATC)  dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9400TA 7.6% 395.8 / 14.1 143F61I (TTC→ATC) dinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base T (70/62);  new base A (4/7);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 144 total aligned reads displayed.

TAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTG  >  NC_000913/3817791‑3818087
                                                                                                                                                     |                                                                                                                                                   
      tGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCgagg                                                                                                                                               >  1:242752/1‑150 (MQ=255)
      tGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCgagg                                                                                                                                               <  1:1688875/150‑1 (MQ=255)
           aaTGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATc                                                                                                                                          <  2:1135764/150‑1 (MQ=255)
           aaTGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATAGTCACTTTATGCCTGTATTAACGATCGCCAAAGCAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGAta                                                                                                                                          >  1:757534/1‑149 (MQ=255)
                 aaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                    >  1:1329915/1‑150 (MQ=255)
                 aaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                    <  1:1584926/150‑1 (MQ=255)
                  aCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACATTTCGAGGAGATCCTCGatat                                                                                                                                  <  2:196736/150‑1 (MQ=255)
                   cGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                                                                  <  2:291201/150‑1 (MQ=255)
                   cGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGGCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                                                                  <  1:421/150‑1 (MQ=255)
                    ggAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATg                                                                                                                                 >  1:1251206/1‑150 (MQ=255)
                    ggAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCGCGATATg                                                                                                                                 >  1:1257910/1‑150 (MQ=255)
                     gAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATgg                                                                                                                                <  2:1647832/150‑1 (MQ=255)
                     gAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACATTGGGCACACAATTGTTGACCATATCGAGGATCTCCTCGAAATGGt                                                                                                                               <  2:49954/150‑1 (MQ=255)
                      aGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGt                                                                                                                               >  1:415978/1‑150 (MQ=255)
                      aGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGt                                                                                                                               <  2:210889/150‑1 (MQ=255)
                        cttaaCAATGGTCAGCAAGACATCTTGGGAAACTAGTGGGTTATTCAGAGTATCGTCAGTTTCTACCTGTATTAAGGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACTATGGAGGACCATATCGAGTATCGCCTCGAATTGGTCa                                                                                                                             <  2:777677/147‑1 (MQ=255)
                        cAGAACAATGGGCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTTTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGTATCTCCTCGAAATGGTCa                                                                                                                             <  1:1371661/150‑1 (MQ=255)
                         aGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCat                                                                                                                            <  2:643240/150‑2 (MQ=255)
                          gAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaaa                                                                                                                           >  2:1159448/1‑150 (MQ=255)
                          caaCAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCatc                                                                                                                           <  1:775407/149‑3 (MQ=255)
                            aCAATGGTCAGCAAGACAACTTGGGAAACTACTGTGTTATTCAGAGTATCGTCACTATAAACCTGTATTACCGCGCGCCAAAGCAGCCTGTGAGAACACTGGTCACACAATTGATGACCATTTCTAGGAGATCCTCGATATGGTCAAAAt                                                                                                                         >  1:266493/1‑150 (MQ=255)
                             cAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAAtt                                                                                                                        >  1:635317/1‑150 (MQ=255)
                              aaTGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                       <  1:667318/150‑1 (MQ=255)
                               attgTCAGAAAGACAACTTTGGAAACTACTGGGTTATTCATAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                      <  1:651250/147‑1 (MQ=255)
                               aTGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                      >  2:721343/1‑150 (MQ=255)
                                 ggTCAGCAAGACAAGTTGGGAAACTACGGGGTTATTCAGAGTATCGTCACTTTATACCTGTCTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                    <  2:1041660/150‑1 (MQ=255)
                                 ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                    >  1:399628/1‑150 (MQ=255)
                                 ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                    >  1:1122669/1‑150 (MQ=255)
                                 ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCCTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                    <  1:722699/150‑1 (MQ=255)
                                   tCAGCAAGACAACTTGGTAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTCTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGAGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCa                                                                                                                  >  1:288748/1‑150 (MQ=255)
                                   tCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCa                                                                                                                  <  1:99098/150‑1 (MQ=255)
                                    cAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaa                                                                                                                 >  1:1636690/1‑150 (MQ=255)
                                     aGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGCAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaaa                                                                                                                >  2:100827/1‑150 (MQ=255)
                                      gCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                               >  2:726804/1‑150 (MQ=255)
                                      gCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                               <  1:770233/150‑1 (MQ=255)
                                       cAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATAAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATg                                                                                                              <  2:1186000/150‑1 (MQ=255)
                                        aaGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGc                                                                                                             >  2:355892/1‑150 (MQ=255)
                                          gACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCa                                                                                                           >  2:555763/1‑150 (MQ=255)
                                              aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAag                                                                                                       >  1:1484273/1‑150 (MQ=255)
                                               cTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                      <  1:314676/150‑1 (MQ=255)
                                                ttGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAagag                                                                                                     <  2:1719904/150‑1 (MQ=255)
                                                  gggAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCGAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                   >  2:1702876/1‑150 (MQ=255)
                                                         tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAg                                                                                            <  1:739870/150‑1 (MQ=255)
                                                         tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAg                                                                                            >  2:1556840/1‑150 (MQ=255)
                                                          actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                 <  1:408633/144‑1 (MQ=255)
                                                          actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                 >  2:408633/1‑144 (MQ=255)
                                                          actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAAAGAGCGCCAAAGCAGCCTGTGAGAACAGTGGACACACAATTAATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATGGGCTCAAATGCCAAAAGAACATTAAAAgc                                                                                           >  2:283866/1‑149 (MQ=255)
                                                             gggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAGATGCCAAAAGAGCATTAAAAGACAt                                                                                        <  1:1078416/150‑1 (MQ=255)
                                                                ttATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGt                                                                                     <  2:1442097/150‑1 (MQ=255)
                                                                ttATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGt                                                                                     >  2:616363/1‑150 (MQ=255)
                                                                 tATTCAGAGTATCGTCACTTTATACCTGTTTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTa                                                                                    >  2:1531562/1‑150 (MQ=255)
                                                                   ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAct                                                                                  >  1:1455583/1‑150 (MQ=255)
                                                                   ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAct                                                                                  >  1:449642/1‑150 (MQ=255)
                                                                     cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAatct                                                                                >  1:381453/1‑150 (MQ=255)
                                                                       gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTcc                                                                              >  1:1007769/1‑150 (MQ=255)
                                                                        agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCAAAAATAGCATTAAAAGACATCGTACTCTccc                                                                             >  2:871366/1‑150 (MQ=255)
                                                                             cGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTAt                                                                        >  1:545888/1‑150 (MQ=255)
                                                                                cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGTGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                     <  2:974108/150‑1 (MQ=255)
                                                                                cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                     <  2:1586847/150‑1 (MQ=255)
                                                                                        aCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCACCCGCTATGCCTGTTACCt                                                             >  1:1257483/1‑150 (MQ=255)
                                                                                          cTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                           >  2:832429/1‑150 (MQ=255)
                                                                                              aTTAACGCGCGACAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagt                                                       <  2:1198719/150‑1 (MQ=255)
                                                                                               ttAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta                                                      <  2:36487/150‑1 (MQ=255)
                                                                                                tAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTAc                                                     <  1:1058630/150‑1 (MQ=255)
                                                                                                 aaCGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACa                                                    <  2:1036773/150‑1 (MQ=255)
                                                                                                   cgcgcgCCAAAGAAGCCTGCGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                  >  1:472100/1‑150 (MQ=255)
                                                                                                     cgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAAc                                                <  2:1288473/150‑1 (MQ=255)
                                                                                                       cgcCAAAGAAGCCTGAGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGCTCCGCGATATGGTCAAAATTGGCTCAACTGCCAAAAGAGCATTAAAAGACATCGTACTCACCCGCTATGCCTGTTACCTGGTAGTACAAAACgg                                              >  2:265924/1‑150 (MQ=255)
                                                                                                        gcaaaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATACTCGATATGGTCAAAATAGGCTCAAATGCCAAAAGAGCAATAAAAGACATCGTACTCTCCCCATAAGCCTGATACCAGGAAGTACAAAAcccc                                             >  2:63507/4‑147 (MQ=255)
                                                                                                        gcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                             <  2:415978/150‑1 (MQ=255)
                                                                                                        gcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                             <  1:614188/150‑1 (MQ=255)
                                                                                                        gcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                             >  1:1003762/1‑150 (MQ=255)
                                                                                                           aaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAc                                          >  1:758937/1‑150 (MQ=255)
                                                                                                                    tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGTAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                 >  2:1595426/1‑150 (MQ=255)
                                                                                                                    tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                 >  2:1011448/1‑150 (MQ=255)
                                                                                                                       gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                              >  2:1627336/1‑150 (MQ=255)
                                                                                                                         gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGt                            <  1:1406267/150‑1 (MQ=255)
                                                                                                                         caaCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGt                            <  1:917278/149‑1 (MQ=255)
                                                                                                                           aCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCa                          <  2:1122669/150‑1 (MQ=255)
                                                                                                                           aCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCa                          >  2:505885/1‑150 (MQ=255)
                                                                                                                            cAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCAt                         >  2:738044/1‑150 (MQ=255)
                                                                                                                             aGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCAtt                        <  2:1492155/150‑1 (MQ=255)
                                                                                                                               tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGc                      <  1:1556840/150‑1 (MQ=255)
                                                                                                                               tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGc                      <  2:1329915/150‑1 (MQ=255)
                                                                                                                                 gTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCgg                    <  2:1007769/150‑1 (MQ=255)
                                                                                                                                  tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggc                   >  1:874600/1‑150 (MQ=255)
                                                                                                                                   cacacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggca                  <  1:100827/150‑1 (MQ=255)
                                                                                                                                    acacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcag                 <  1:721343/150‑1 (MQ=255)
                                                                                                                                       caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggc              >  2:64425/1‑150 (MQ=255)
                                                                                                                                        aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggca             <  1:155817/150‑1 (MQ=255)
                                                                                                                                        aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAATCTCCCGCTATGCCGTTTACCTGGTAGTACAAAACGACGACACTGAGAAACCGGTCATTGCggcagggca             >  2:1568991/1‑150 (MQ=255)
                                                                                                                                        aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGGACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggca             <  1:871366/150‑1 (MQ=255)
                                                                                                                                              tgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTa       >  1:551521/1‑150 (MQ=255)
                                                                                                                                               gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAt      >  2:5662/1‑150 (MQ=255)
                                                                                                                                                aCCATTTCGAGGAGATCCTCGATATGGTAAAAATTGGCTCAAAAGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAACTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt     >  2:47666/1‑150 (MQ=255)
                                                                                                                                                  cATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtttt   >  1:1316906/1‑150 (MQ=255)
                                                                                                                                                   aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGTTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg  <  2:989919/150‑1 (MQ=255)
                                                                                                                                                   aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg  >  1:559291/1‑150 (MQ=255)
                                                                                                                                                   aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg  >  2:918380/1‑150 (MQ=255)
                                                                                                                                                   aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGACAGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg  >  1:1104703/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTG  >  NC_000913/3817791‑3818087

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: