Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,050,226 C→T 8.5% 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 8.5% 351.9 / 14.0 118S180S (TCG→TCAyihAcell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (53/54);  new base T (4/6);  total (58/60)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

Only 100 of 119 total aligned reads displayed.

ATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCACTC  >  NC_000913/4050078‑4050375
                                                                                                                                                    |                                                                                                                                                     
aTTGAATTTACGGGCCGGATACGCCACATCCGGCACCAGCATTAAGGCAAGAAAATAATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCAAGTTTCTGCCGCAGCTTGTCCACGACTTGTTTCTTCACCg                                                                                                                                                      >  2:1570693/1‑150 (MQ=255)
 ttGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACga                                                                                                                                                     >  2:530985/1‑150 (MQ=255)
    aaTTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaaga                                                                                                                                                  >  1:1341740/1‑150 (MQ=255)
      tttACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGaaa                                                                                                                                                >  1:1565206/1‑150 (MQ=255)
        tACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaaca                                                                                                                                              <  1:1122511/150‑3 (MQ=255)
         aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTACACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAAAGAAGAAAACg                                                                                                                                             >  2:881350/1‑150 (MQ=255)
         aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaa                                                                                                                                             <  1:913784/150‑4 (MQ=255)
         aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaa                                                                                                                                             <  1:1187580/150‑4 (MQ=255)
           gggCCGGATACGCCACCTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCGGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaagg                                                                                                                                           <  1:1320517/150‑6 (MQ=255)
           gggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacaagg                                                                                                                                           <  2:521403/150‑6 (MQ=255)
            ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGttt                                                                                                                                          >  2:1016933/1‑150 (MQ=255)
            ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGttt                                                                                                                                          >  1:1305569/1‑150 (MQ=255)
            ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGAATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAACCtttt                                                                                                                                          >  2:250891/1‑150 (MQ=255)
             gCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCGACGCCTTGTTTCTTCAACGAAGAAAACGtttg                                                                                                                                         >  1:476122/1‑149 (MQ=255)
             gCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTc                                                                                                                                         <  1:928811/150‑1 (MQ=255)
                 gATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAAcc                                                                                                                                     <  2:1019867/150‑1 (MQ=255)
                    aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                  <  1:469906/150‑1 (MQ=255)
                      gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcac                                                                                                                                <  1:310672/150‑1 (MQ=255)
                      gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcac                                                                                                                                >  2:1103517/1‑150 (MQ=255)
                           aTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCAc                                                                                                                           <  2:4533/150‑1 (MQ=255)
                                 cacaAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGCAAACGTTTCACCCTGCCCAGCCCCGTTaa                                                                                                                     >  2:1146924/1‑150 (MQ=255)
                                 cacaAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTaa                                                                                                                     <  2:1059021/150‑1 (MQ=255)
                                   caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaaca                                                                                                                                              <  2:1118179/123‑3 (MQ=255)
                                   caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaaca                                                                                                                                              >  1:1118179/1‑121 (MQ=255)
                                   caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAAc                                                                                                                   <  2:1166100/150‑1 (MQ=255)
                                   caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAAc                                                                                                                   >  1:1155820/1‑150 (MQ=255)
                                       cATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAACCGTTTCAACCTGCACGTCACCGATAAACGCCa                                                                                                               >  2:288830/1‑150 (MQ=255)
                                       cATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCATTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCa                                                                                                               >  1:151869/1‑150 (MQ=255)
                                       aaTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCa                                                                                                               <  2:1269382/149‑1 (MQ=255)
                                             ggCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAg                                                                                                         >  1:361924/1‑150 (MQ=255)
                                             ggCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAg                                                                                                         >  1:1255421/1‑150 (MQ=255)
                                              gCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGTCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGc                                                                                                        <  1:1048360/150‑1 (MQ=255)
                                               cAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCt                                                                                                       >  1:1344742/1‑150 (MQ=255)
                                                  gAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACCCCAGTACAGCTTcc                                                                                                    >  2:138318/1‑149 (MQ=255)
                                                   aaaaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCAc                                                                                                   <  1:1726999/150‑1 (MQ=255)
                                                     aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                                 <  1:920690/150‑1 (MQ=255)
                                                     aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                                 <  2:1303541/150‑1 (MQ=255)
                                                     aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                                 <  2:1255421/150‑1 (MQ=255)
                                                        tattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCCCGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCCCCGTTAAACGCCAGTACAGCTTCACGCAcc                                                                                              >  2:631169/1‑150 (MQ=255)
                                                        tattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCACCGCCTTGTTTCTTCACCGGAGAAATCGTTTAAACCTGCACCTCACCGGTAAACGCCAGTACAGCTTCAAGCAcc                                                                                              >  1:863199/1‑150 (MQ=255)
                                                        tattCGCCGTCCTGCGTTTCTTCTACAGGCCGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAcc                                                                                              <  2:1710808/150‑1 (MQ=255)
                                                           tCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCata                                                                                           <  1:1632423/150‑1 (MQ=255)
                                                           tCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCata                                                                                           <  1:1584866/150‑1 (MQ=255)
                                                            cGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACTTCACCGTTAAACGCCATTACAGCTTCACGCACCatat                                                                                          >  2:877582/1‑150 (MQ=255)
                                                             gCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTTCCGCTTCACGCCCCATAtt                                                                                         >  1:834732/1‑150 (MQ=255)
                                                                 tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                     <  2:1589493/150‑1 (MQ=255)
                                                                 tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCAGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                     <  1:130452/150‑1 (MQ=255)
                                                                   cTGCGTTTCTTCTACAGGCTGCATCTCGCTAATCCAGGTATCCAGTTTCTTCCGCAGCTTGCCCACGCCTTGTTTCTTCAACGAAGATAACGTTTCAACCTGCACATCACCCTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTg                                                                                   <  1:1179499/150‑1 (MQ=255)
                                                                     gCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATCTTCAATTgcg                                                                                 >  1:1638190/1‑150 (MQ=255)
                                                                       gTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCt                                                                               >  1:660679/1‑150 (MQ=255)
                                                                           cttctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTAc                                                                           >  2:754048/1‑150 (MQ=255)
                                                                            ttctACAGGCTGCATATCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATANTCAATTGCGCTTTACg                                                                          <  1:1570693/150‑1 (MQ=255)
                                                                              ctACAGGCTGCATCTCGCTAAACCAGGTAACCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTg                                                                        >  1:1034968/1‑150 (MQ=255)
                                                                              caaCAGGCTGCATCTTGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTg                                                                        <  2:254914/148‑1 (MQ=255)
                                                                                 cAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGTGCTTTACGTgcgc                                                                     <  1:223989/150‑1 (MQ=255)
                                                                                   ggCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACATGCACATCACCGTTATACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGcgccc                                                                   >  2:1529683/1‑149 (MQ=255)
                                                                                         aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGcc                                                             <  2:618017/150‑1 (MQ=255)
                                                                                         aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGcc                                                             >  1:390986/1‑150 (MQ=255)
                                                                                         aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGACAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGcc                                                             >  1:1226349/1‑150 (MQ=255)
                                                                                           ctcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTTCGATTTACTTGCGCCGCTTGCCAg                                                           >  2:287790/1‑150 (MQ=255)
                                                                                           ctcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                           <  1:1403914/150‑1 (MQ=255)
                                                                                           ctcGCTAAACCAGGTATCCAGGTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                           <  2:1148598/150‑1 (MQ=255)
                                                                                             cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                         <  1:502005/150‑1 (MQ=255)
                                                                                             cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                         >  2:1350622/1‑150 (MQ=255)
                                                                                             cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                         <  2:1341740/150‑1 (MQ=255)
                                                                                             cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                         <  1:367034/150‑1 (MQ=255)
                                                                                             cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGAGCCGCTTGCCAGtt                                                         >  2:1603554/1‑150 (MQ=255)
                                                                                               cTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTg                                                       <  1:695788/150‑1 (MQ=255)
                                                                                                tAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGt                                                      <  1:1016933/150‑1 (MQ=255)
                                                                                                   aCCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCg                                                   >  1:765959/1‑150 (MQ=255)
                                                                                                    ccAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTTAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGc                                                  >  1:1559650/1‑150 (MQ=255)
                                                                                                      aGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTGTCCACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCtt                                                >  2:1098563/1‑150 (MQ=255)
                                                                                                       ggTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                               >  1:523652/1‑150 (MQ=255)
                                                                                                       ggTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                               <  2:1550766/150‑1 (MQ=255)
                                                                                                         tATCCAGGTTCTGCCGCAGCTTGTCCACGCCTTGGTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTgg                                             <  2:89139/150‑1 (MQ=255)
                                                                                                           tCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCACCCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGATTTggaa                                           >  2:862798/1‑148 (MQ=255)
                                                                                                            ccAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAATACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTca                                          >  1:115044/1‑150 (MQ=255)
                                                                                                               gTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCTACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcagca                                       <  1:1103517/150‑1 (MQ=255)
                                                                                                               gTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcagca                                       <  1:1325035/150‑1 (MQ=255)
                                                                                                                 ttCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCACCCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCGTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTcatgagc                                     >  2:1118443/1‑145 (MQ=255)
                                                                                                                     gCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCa                                 >  1:1720272/1‑150 (MQ=255)
                                                                                                                     gCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCa                                 >  2:662830/1‑150 (MQ=255)
                                                                                                                        gCAGCTTGTCCACGCCTTGTTTCTTCAACGATGAAAACGTTTCAACCTGCACATCACCGTTTAAAGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCGAGTTTGTCCGCTTTGGTCAGCCGCACCAgca                              >  1:312023/1‑148 (MQ=255)
                                                                                                                            cTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGc                          <  2:532267/150‑1 (MQ=255)
                                                                                                                              tGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGa                        <  1:1304563/150‑1 (MQ=255)
                                                                                                                               gTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGat                       <  2:1447838/150‑1 (MQ=255)
                                                                                                                                tCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGata                      <  2:625078/150‑1 (MQ=255)
                                                                                                                                 ccACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCAGATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGctat                     >  2:19461/1‑150 (MQ=255)
                                                                                                                                    cGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACGTCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGc                  <  1:956938/150‑1 (MQ=255)
                                                                                                                                     aaacgttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACATACGCCGCTTGGCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATAGCt                 >  2:114960/7‑150 (MQ=255)
                                                                                                                                         gttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt             >  2:399555/3‑150 (MQ=255)
                                                                                                                                         gttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt             >  2:1298485/3‑150 (MQ=255)
                                                                                                                                             tcttcAACCAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAATGCGCATTACGTGCGCAGATTGCCAGTTTGTCCGCTATGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCg         >  1:960201/1‑150 (MQ=255)
                                                                                                                                                tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc      >  2:1583597/1‑150 (MQ=255)
                                                                                                                                                tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc      <  1:1182522/150‑1 (MQ=255)
                                                                                                                                                tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc      <  1:754048/150‑1 (MQ=255)
                                                                                                                                                 cAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCa     >  2:1296817/1‑150 (MQ=255)
                                                                                                                                                  aaCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAATGGGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccacc    >  2:851333/1‑147 (MQ=255)
                                                                                                                                                   aCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGGTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCACt   <  1:881350/150‑1 (MQ=255)
                                                                                                                                                    cGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCACTc  >  1:1358100/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCACTC  >  NC_000913/4050078‑4050375

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: