Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,050,226 C→T 5.1% S180S (TCG→TCA 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,2260CT 5.1% 443.5 / 6.1 135S180S (TCG→TCAyihAcell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (60/68);  new base T (2/5);  total (62/73)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 137 total aligned reads displayed.

TATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050077‑4050373
                                                                                                                                                     |                                                                                                                                                   
 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)
    gAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCAAGCCTTGTTTCTTCAACgaaa                                                                                                                                                 >  2:302438/1‑149 (MQ=255)
     aaTTTACGGGCCGGATACGCCACATCCGACACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaaga                                                                                                                                                <  1:129047/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                                                                                                                                         >  2:565693/1‑150 (MQ=255)
            gggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGtt                                                                                                                                         <  1:951017/150‑1 (MQ=255)
            gggCCGGATACGCCACATCCGGCACAAGAATTTAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaagg                                                                                                                                         <  1:1105095/150‑6 (MQ=255)
             ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaaggc                                                                                                                                        <  2:1257990/150‑7 (MQ=255)
             ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGttt                                                                                                                                        <  1:1530718/150‑1 (MQ=255)
               ccGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa                                                                                                                                      <  2:1201521/150‑1 (MQ=255)
               ccGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa                                                                                                                                      >  2:254046/1‑150 (MQ=255)
                     aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 >  2:946532/1‑133 (MQ=255)
                      cGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTTTCCACGCCTTGTTTCTTCGTTGAAGaaac                                                                                                                                             <  2:820133/136‑2 (MQ=255)
                      cGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGca                                                                                                                               >  1:1501374/1‑150 (MQ=255)
                        ccACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcaca                                                                                                                             <  1:1244593/150‑1 (MQ=255)
                           caTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCa                                                                                                                          <  2:1410520/150‑1 (MQ=255)
                            aTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATccc                                                                                                                         >  2:903865/1‑148 (MQ=255)
                                 gCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTa                                                                                                                    <  1:39768/150‑1 (MQ=255)
                                 gCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTa                                                                                                                    <  2:1539053/150‑1 (MQ=255)
                                   acaAGCATTAAGGCAAGAAAATTCTTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTaaa                                                                                                                  <  1:148434/150‑1 (MQ=255)
                                      aGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGc                                                                                                               >  1:1622122/1‑150 (MQ=255)
                                         aTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAg                                                                                                            <  1:1438505/150‑1 (MQ=255)
                                         aTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAg                                                                                                            <  1:1669811/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)
                                                        ttattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAc                                                                                             <  2:1696912/150‑1 (MQ=255)
                                                        ttattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAc                                                                                             <  2:1445459/150‑1 (MQ=255)
                                                         tattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAcc                                                                                            <  1:732854/150‑1 (MQ=255)
                                                           ttCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCat                                                                                          <  1:1580566/150‑1 (MQ=255)
                                                             cGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCatat                                                                                        <  1:1403017/150‑1 (MQ=255)
                                                             cGCCGTCCTGCGTTTCTTCGACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGGTAAACGCCAGTCCAGCTTCACGCACCatat                                                                                        >  2:491417/1‑150 (MQ=255)
                                                                cGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCa                                                                                     <  1:847306/150‑1 (MQ=255)
                                                                 gTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCaa                                                                                    <  1:911171/150‑1 (MQ=255)
                                                                  tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                   >  1:499010/1‑150 (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)
                                                                           tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          <  2:1674262/150‑1 (MQ=255)
                                                                           tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          >  1:361768/1‑150 (MQ=255)
                                                                           tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          >  1:368009/1‑150 (MQ=255)
                                                                           tcttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTa                                                                          >  1:156181/1‑150 (MQ=255)
                                                                               ctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTCCGTg                                                                      >  2:310180/1‑150 (MQ=255)
                                                                                         cATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGc                                                            >  1:502847/1‑150 (MQ=255)
                                                                                           tctcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCTCGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCTTTAAACGCCAGTACAGCATCACGCGCCATATTCAATTGCTATTTAAGTGCGCAGCTTgacc                                                          >  1:815972/1‑147 (MQ=255)
                                                                                            ctcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                         >  2:189044/1‑150 (MQ=255)
                                                                                              cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                       <  2:1459495/150‑1 (MQ=255)
                                                                                              cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCAGATTCAATTGAGCTTTACGTGCGGCGCTTGCCAGtt                                                       >  1:208607/1‑150 (MQ=255)
                                                                                                cTAAACCAGGTATCCAGATGCTGCCGCAGCATGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCTCTTGCCAGTTTg                                                     <  2:1092628/150‑1 (MQ=255)
                                                                                                  aaaCCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTc                                                   <  2:2351/150‑1 (MQ=255)
                                                                                                          tATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAACCGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTgg                                           <  1:1231025/150‑1 (MQ=255)
                                                                                                          tATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTgg                                           <  1:935762/150‑1 (MQ=255)
                                                                                                             ccAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTca                                        <  2:1004424/150‑1 (MQ=255)
                                                                                                             ccAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCACTTTGGTca                                        >  1:1207439/1‑150 (MQ=255)
                                                                                                              cAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcag                                       >  2:1423/1‑150 (MQ=255)
                                                                                                               aGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcagc                                      <  2:620471/150‑1 (MQ=255)
                                                                                                                gTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCTCATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTTAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTTGTcagca                                     <  1:67599/150‑1 (MQ=255)
                                                                                                                   tCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAgcagca                                  <  1:565693/150‑1 (MQ=255)
                                                                                                                    cTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCAc                                 <  1:683201/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)
                                                                                                                         gtaGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTATGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGaa                            >  2:1707047/3‑150 (MQ=255)
                                                                                                                         gCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGaa                            >  1:517748/1‑150 (MQ=255)
                                                                                                                          cAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAAc                           <  1:61855/150‑1 (MQ=255)
                                                                                                                           aGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACg                          <  2:156181/150‑1 (MQ=255)
                                                                                                                           aGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACg                          <  2:361768/150‑1 (MQ=255)
                                                                                                                           aGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACg                          <  2:368009/150‑1 (MQ=255)
                                                                                                                             cTTGTCCAAGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACAATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTCGTCAGCAGCACCAGAACGGc                        >  2:1196616/1‑150 (MQ=255)
                                                                                                                              ttGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCCGTACAGCTTCACGCACCATATTCAATTGCGCTTTTCGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCg                       <  2:966222/150‑1 (MQ=255)
                                                                                                                                gTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCCCCGTTAAACGCCAGTACAGCTTCACGCCCCATATTCAATTGCGCTTTACGTGCGCCGCTTGACAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGag                     >  2:1102214/1‑149 (MQ=255)
                                                                                                                                gTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCCCCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTCCGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGCACGGCGaa                     >  1:299445/1‑149 (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)
                                                                                                                                     cGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGc                <  1:982695/150‑1 (MQ=255)
                                                                                                                                       ccTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTg              >  1:564822/1‑150 (MQ=255)
                                                                                                                                         cgttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACAAGAACGGCGATATTGCTGTc            >  2:1611311/4‑150 (MQ=255)
                                                                                                                                          tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAATCGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           <  1:254046/150‑1 (MQ=255)
                                                                                                                                          tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           <  2:129993/150‑1 (MQ=255)
                                                                                                                                           gTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTa          >  1:198657/1‑150 (MQ=255)
                                                                                                                                            tttCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTAc         <  2:251155/150‑1 (MQ=255)
                                                                                                                                                ttcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGcc     <  1:228842/150‑1 (MQ=255)
                                                                                                                                                  cAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCa   >  2:613954/1‑150 (MQ=255)
                                                                                                                                                   aaCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAc  <  2:499010/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050077‑4050373

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