Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,050,225 A→T 5.1% L181M (TTG→ATG) ‡ yihA ← cell division GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,050,2250AT 5.1% 396.5 / 7.7 136L181M (TTG→ATG) ‡yihAcell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (61/68);  new base T (2/5);  total (63/73)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 138 total aligned reads displayed.

CTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050076‑4050373
                                                                                                                                                     |                                                                                                                                                    
 tATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAAc                                                                                                                                                     <  1:69983/150‑1 (MQ=255)
 tATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAAc                                                                                                                                                     >  2:228842/1‑150 (MQ=255)
  aTTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACg                                                                                                                                                    <  1:840993/150‑1 (MQ=255)
  aTTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACg                                                                                                                                                    >  2:683201/1‑150 (MQ=255)
     gAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 >  1:1410520/1‑150 (MQ=255)
     gAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 <  1:1268698/150‑1 (MQ=255)
       ttttACGGGCCGGATACGCCACATCCGGCACTAGCTTTAAGGCAAGAATATTATTCGCCGTCCTTCGTTTCTTCTACAGGCTGCATCTCGCTTAACCAGGTATCCAGTTTCTTCCGCAGCTTGTCCACGCCTTGTTTCTTCAACtaagaa                                                                                                                                               <  1:12669/149‑1 (MQ=255)
       aTTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaagaa                                                                                                                                               >  1:1561985/1‑150 (MQ=255)
        tttACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGaaa                                                                                                                                              >  2:1403017/1‑150 (MQ=255)
          tACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAAc                                                                                                                                            <  1:1103503/150‑1 (MQ=255)
           aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACg                                                                                                                                           <  2:550257/150‑1 (MQ=255)
           aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGCATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaa                                                                                                                                           <  1:835154/150‑4 (MQ=255)
            cGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaag                                                                                                                                          <  2:981448/150‑5 (MQ=255)
            cGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCGGCATCTCGCTAAACCAGGTATACAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAATCGt                                                                                                                                          >  2:190867/1‑150 (MQ=255)
             gggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGtt                                                                                                                                         <  1:951017/150‑1 (MQ=255)
             gggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGtt                                                                                                                                         >  2:565693/1‑150 (MQ=255)
             gggCCGGATACGCCACATCCGGCACAAGAATTTAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaagg                                                                                                                                         <  1:1105095/150‑6 (MQ=255)
              ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGttt                                                                                                                                        <  1:1530718/150‑1 (MQ=255)
                ccGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa                                                                                                                                      >  2:254046/1‑150 (MQ=255)
                     tACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTg                                                                                                                                 <  1:1517599/150‑1 (MQ=255)
                     tACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTg                                                                                                                                 <  2:1720609/150‑1 (MQ=255)
                      aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 >  2:946532/1‑133 (MQ=255)
                      aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 <  1:946532/133‑1 (MQ=255)
                      aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                >  1:54797/1‑150 (MQ=255)
                       cGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTTTCCACGCCTTGTTTCTTCGTTGAAGaaac                                                                                                                                             <  2:820133/136‑2 (MQ=255)
                       cGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaac                                                                                                                                             >  1:820133/1‑135 (MQ=255)
                        gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcac                                                                                                                              >  1:129993/1‑150 (MQ=255)
                         ccACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcaca                                                                                                                             <  1:1244593/150‑1 (MQ=255)
                            caTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCa                                                                                                                          >  2:389536/1‑150 (MQ=255)
                             aTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATccc                                                                                                                         >  2:903865/1‑148 (MQ=255)
                                  gCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTa                                                                                                                    <  1:39768/150‑1 (MQ=255)
                                    acaAGCATTAAGGCAAGAAAATTCTTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTaaa                                                                                                                  <  1:148434/150‑1 (MQ=255)
                                     caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACACCACCGTTAAAc                                                                                                                 >  1:501245/1‑150 (MQ=255)
                                     aaaGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAAc                                                                                                                 <  2:1064307/149‑1 (MQ=255)
                                       aGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGc                                                                                                               >  1:1622122/1‑150 (MQ=255)
                                          aTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAg                                                                                                            <  1:389536/150‑1 (MQ=255)
                                          aTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAg                                                                                                            <  1:1438505/150‑1 (MQ=255)
                                              aGGCATGAAAATTCATCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTGCACGGCTTGTTTCTTCAACGAAGAAATCGTTTCAACCTGCACATCACCGTTTAACGCCAGTACa                                                                                                        <  2:794695/150‑1 (MQ=255)
                                                 cAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCt                                                                                                     <  2:1590719/150‑1 (MQ=255)
                                                  aaGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCtt                                                                                                    <  1:1625205/150‑1 (MQ=255)
                                                      aaaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACg                                                                                                <  1:467328/150‑1 (MQ=255)
                                                         ttattCGCCGTCCTGCGTTTCTTCTACAGGCTGCCTCTCGCTAAACCAGGTATCCAGTTTCTGCCTCAGCTTGTCAACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGGTAAACGCCAGTACAGCTTCACGCAc                                                                                             >  1:251155/1‑150 (MQ=255)
                                                          tattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAcc                                                                                            <  1:732854/150‑1 (MQ=255)
                                                            ttCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCat                                                                                          <  1:1580566/150‑1 (MQ=255)
                                                              cGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCatat                                                                                        <  1:1403017/150‑1 (MQ=255)
                                                                  gTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCaa                                                                                    <  1:911171/150‑1 (MQ=255)
                                                                   tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                   >  1:1459495/1‑150 (MQ=255)
                                                                    ccTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAtt                                                                                  <  2:1215343/150‑1 (MQ=255)
                                                                        cGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTgcgc                                                                              >  1:627068/1‑150 (MQ=255)
                                                                         gTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCt                                                                             <  1:1543969/150‑1 (MQ=255)
                                                                         gTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCt                                                                             <  1:1716615/150‑1 (MQ=255)
                                                                            tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          >  1:156181/1‑150 (MQ=255)
                                                                            tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          <  2:1674262/150‑1 (MQ=255)
                                                                            tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          >  1:361768/1‑150 (MQ=255)
                                                                                ctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTCCGTg                                                                      >  2:310180/1‑150 (MQ=255)
                                                                                ctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATAACCGTTAAACGCCAGTACAGCTTCCCGCACCATATTCAATTGAGCTTTACGTg                                                                      >  2:281302/1‑150 (MQ=255)
                                                                                  aCAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTgcg                                                                    <  2:1666626/150‑1 (MQ=255)
                                                                                        tGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCtt                                                              >  2:1414546/1‑150 (MQ=255)
                                                                                         gCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTTTTTCATCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTg                                                             <  1:302438/150‑1 (MQ=255)
                                                                                         gCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTg                                                             <  2:1222874/150‑1 (MQ=255)
                                                                                          cATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGc                                                            >  1:502847/1‑150 (MQ=255)
                                                                                             ctcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                         <  2:1468944/150‑1 (MQ=255)
                                                                                               cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                       <  2:1459495/150‑1 (MQ=255)
                                                                                               cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCAGATTCAATTGAGCTTTACGTGCGGCGCTTGCCAGtt                                                       >  1:208607/1‑150 (MQ=255)
                                                                                                 cTAAACCAGGTATCCAGATGCTGCCGCAGCATGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCTCTTGCCAGTTTg                                                     <  2:1092628/150‑1 (MQ=255)
                                                                                                        aGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCtt                                              >  1:624445/1‑150 (MQ=255)
                                                                                                         ggTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                             <  2:1303405/150‑1 (MQ=255)
                                                                                                          gTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTg                                            <  2:551489/150‑1 (MQ=255)
                                                                                                           tATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAACCGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTgg                                           <  1:1231025/150‑1 (MQ=255)
                                                                                                              ccAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTca                                        <  2:1004424/150‑1 (MQ=255)
                                                                                                              ccAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCACTTTGGTca                                        >  1:1207439/1‑150 (MQ=255)
                                                                                                               cAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcag                                       >  2:61855/1‑150 (MQ=255)
                                                                                                               cAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcag                                       >  2:1423/1‑150 (MQ=255)
                                                                                                                aGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCCCATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcagc                                      >  2:1150734/1‑150 (MQ=255)
                                                                                                                 gTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCTCATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTTAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTTGTcagca                                     <  1:67599/150‑1 (MQ=255)
                                                                                                                    tCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAgcagca                                  <  1:565693/150‑1 (MQ=255)
                                                                                                                      tGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCAcc                                <  1:1225579/150‑1 (MQ=255)
                                                                                                                       gCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCa                               >  1:1347012/1‑150 (MQ=255)
                                                                                                                        ccGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAg                              >  1:936139/1‑150 (MQ=255)
                                                                                                                         cGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCGCGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGa                             >  2:1491711/1‑150 (MQ=255)
                                                                                                                         cGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGa                             >  1:360916/1‑150 (MQ=255)
                                                                                                                         cGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGa                             >  2:340751/1‑150 (MQ=255)
                                                                                                                          gtaGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTATGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGaa                            >  2:1707047/3‑150 (MQ=255)
                                                                                                                            aGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACg                          <  2:361768/150‑1 (MQ=255)
                                                                                                                            aGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACg                          <  2:156181/150‑1 (MQ=255)
                                                                                                                              cTTGTCCAAGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACAATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTCGTCAGCAGCACCAGAACGGc                        >  2:1196616/1‑150 (MQ=255)
                                                                                                                               ttGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCCGTACAGCTTCACGCACCATATTCAATTGCGCTTTTCGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCg                       <  2:966222/150‑1 (MQ=255)
                                                                                                                                  tCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGata                    >  2:181867/1‑150 (MQ=255)
                                                                                                                                   ccACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGatat                   >  2:1236234/1‑150 (MQ=255)
                                                                                                                                     aCGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTg                 <  2:386158/150‑1 (MQ=255)
                                                                                                                                      cGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGc                <  1:1364996/150‑1 (MQ=255)
                                                                                                                                        ccTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTg              >  1:564822/1‑150 (MQ=255)
                                                                                                                                           tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAATCGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           <  1:254046/150‑1 (MQ=255)
                                                                                                                                           tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           <  2:129993/150‑1 (MQ=255)
                                                                                                                                           tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           <  2:1561985/150‑1 (MQ=255)
                                                                                                                                            gTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTa          >  1:198657/1‑150 (MQ=255)
                                                                                                                                             tttCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTAc         <  2:251155/150‑1 (MQ=255)
                                                                                                                                                  tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCTCCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc    <  1:903865/150‑1 (MQ=255)
                                                                                                                                                   cAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCa   >  2:613954/1‑150 (MQ=255)
                                                                                                                                                    aaCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAc  <  2:499010/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050076‑4050373

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