Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,050,224 A→G 8.6% L181S (TTG→TCG) ‡ 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,2240AG 8.6% 323.1 / 16.7 116L181S (TTG→TCG) ‡yihAcell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (53/53);  new base G (4/6);  total (57/59)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 118 total aligned reads displayed.

TCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050075‑4050373
                                                                                                                                                     |                                                                                                                                                     
tCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTCTTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCa                                                                                                                                                       <  1:207999/150‑1 (MQ=255)
   aTTGAATTTACGGGCCGGATACGCCACATCCGGCACCAGCATTAAGGCAAGAAAATAATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCAAGTTTCTGCCGCAGCTTGTCCACGACTTGTTTCTTCACCg                                                                                                                                                    >  2:1570693/1‑150 (MQ=255)
   aTTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTAATCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCTGCTTGTCCACGCCTTGTTTCTTCAACg                                                                                                                                                    <  1:463737/150‑1 (MQ=255)
    ttGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACga                                                                                                                                                   >  2:530985/1‑150 (MQ=255)
     tGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCAACGCCTTGTTTCTTCAACgaa                                                                                                                                                  >  1:1202611/1‑150 (MQ=255)
      gAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 <  2:1536879/150‑1 (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)
              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)
                 ccGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa                                                                                                                                      <  1:158642/150‑1 (MQ=255)
                    gATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAAcc                                                                                                                                   <  2:1019867/150‑1 (MQ=255)
                       aCGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                <  1:469906/150‑1 (MQ=255)
                         gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcac                                                                                                                              >  2:1325035/1‑150 (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)
                                      aaaGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAAc                                                                                                                 <  2:100171/149‑1 (MQ=255)
                                          cATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCATTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCa                                                                                                             >  1:151869/1‑150 (MQ=255)
                                          aaTTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCa                                                                                                             <  2:1269382/149‑1 (MQ=255)
                                                ggCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAg                                                                                                       <  1:748401/150‑1 (MQ=255)
                                                ggCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAg                                                                                                       >  1:1255421/1‑150 (MQ=255)
                                                ggCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAg                                                                                                       >  1:361924/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                                                                                               <  2:1303541/150‑1 (MQ=255)
                                                        aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                               <  2:361924/150‑1 (MQ=255)
                                                        aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                               <  2:1255421/150‑1 (MQ=255)
                                                        aaTTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGc                                                                                               <  1:920690/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:1584866/150‑1 (MQ=255)
                                                              tCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCata                                                                                         <  1:1632423/150‑1 (MQ=255)
                                                               cGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACTTCACCGTTAAACGCCATTACAGCTTCACGCACCatat                                                                                        >  2:877582/1‑150 (MQ=255)
                                                                gCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTTCCGCTTCACGCCCCATAtt                                                                                       >  1:834732/1‑150 (MQ=255)
                                                                  cGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCa                                                                                     <  1:1307821/150‑1 (MQ=255)
                                                                    tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                   <  2:1589493/150‑1 (MQ=255)
                                                                    tCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCAGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAt                                                                                   <  1:130452/150‑1 (MQ=255)
                                                                      cTGCGTTTCTTCTACAGGCTGCATCTCGCTAATCCAGGTATCCAGTTTCTTCCGCAGCTTGCCCACGCCTTGTTTCTTCAACGAAGATAACGTTTCAACCTGCACATCACCCTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTg                                                                                 <  1:1179499/150‑1 (MQ=255)
                                                                        gCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTgcg                                                                               >  2:1182522/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)
                                                                                 ctACAGGCTGCATCTCGCTAAACCAGGTATCCAGTAACTCCCGCAGCTTGTCCACGCCTTGTTTCATAAACGAAGAAAACGTATCAAACTGCACATCACCGTTAAACGCCAGTACACCAACCCGAACCATAATCAATTGCGCTTTACGTg                                                                      >  2:548729/1‑150 (MQ=255)
                                                                                 ctACAGGCTGCATCTCGCTAAACCAGGTAACCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTg                                                                      >  1:1034968/1‑150 (MQ=255)
                                                                                 caaCAGGCTGCATCTTGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTg                                                                      <  2:254914/148‑1 (MQ=255)
                                                                                    cAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGTGCTTTACGTgcgc                                                                   <  1:223989/150‑1 (MQ=255)
                                                                                            aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGcc                                                           <  2:618017/150‑1 (MQ=255)
                                                                                            aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGcc                                                           >  1:390986/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                                                       <  2:1341740/150‑1 (MQ=255)
                                                                                                cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                       >  2:1350622/1‑150 (MQ=255)
                                                                                                cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                       <  1:367034/150‑1 (MQ=255)
                                                                                                cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGtt                                                       <  1:502005/150‑1 (MQ=255)
                                                                                                  cTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTg                                                     <  1:695788/150‑1 (MQ=255)
                                                                                                   tAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGt                                                    <  1:1016933/150‑1 (MQ=255)
                                                                                                          ggcaTCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                             <  2:1641811/147‑1 (MQ=255)
                                                                                                          ggTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                             >  1:523652/1‑150 (MQ=255)
                                                                                                          ggTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCttt                                             <  2:1550766/150‑1 (MQ=255)
                                                                                                           gTATCCAGTTTCTGCCGCAGCTTGTCCACTCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCAAGCACCATATTCAATTGCGCTTTAAGTGCGCCGATTGCCAGTTTGTCCGCTTTg                                            >  2:1104091/1‑150 (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)
                                                                                                                                 tGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGa                      <  1:1304563/150‑1 (MQ=255)
                                                                                                                                  gTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGat                     <  2:1447838/150‑1 (MQ=255)
                                                                                                                                   tCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGata                    <  2:625078/150‑1 (MQ=255)
                                                                                                                                    ccACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGatat                   <  1:887800/150‑1 (MQ=255)
                                                                                                                                    ccACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCAGATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGctat                   >  2:19461/1‑150 (MQ=255)
                                                                                                                                       cGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACGTCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGc                <  1:956938/150‑1 (MQ=255)
                                                                                                                                         ccTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCTCATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCATGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTg              >  1:231298/1‑150 (MQ=255)
                                                                                                                                            tGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAATGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           >  1:1465656/1‑150 (MQ=255)
                                                                                                                                            gttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCt           >  2:399555/3‑150 (MQ=255)
                                                                                                                                                tcttcAACCAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAATGCGCATTACGTGCGCAGATTGCCAGTTTGTCCGCTATGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCg       >  1:960201/1‑150 (MQ=255)
                                                                                                                                                   tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc    <  1:754048/150‑1 (MQ=255)
                                                                                                                                                   tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc    >  2:1583597/1‑150 (MQ=255)
                                                                                                                                                   tcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccc    <  1:1182522/150‑1 (MQ=255)
                                                                                                                                                    cAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCa   >  2:1296817/1‑150 (MQ=255)
                                                                                                                                                     aaCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAAATGGGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGccacc  >  2:851333/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050075‑4050373

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: