Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,772,236 A→T 16.1% pseudogene (831/900 nt) yfjV ← pseudogene, CP4‑57 prophage;Phage or Prophage Related; putative pump protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,772,2360AT16.1% 282.9 / 43.1 113pseudogene (831/900 nt)yfjVpseudogene, CP4‑57 prophage;Phage or Prophage Related; putative pump protein
Reads supporting (aligned to +/- strand):  ref base A (46/48);  new base T (8/10);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 114 total aligned reads displayed.

GATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAT  >  NC_000913/2772137‑2772334
                                                                                                   |                                                                                                  
gATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACa                                                                                                    >  2:2113561/1‑100 (MQ=255)
  tGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAAt                                                                                                  >  2:2125187/1‑100 (MQ=255)
       aTGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCgg                                                                                             <  2:822912/100‑1 (MQ=255)
       aTGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCgg                                                                                             <  1:1267976/100‑1 (MQ=255)
         gTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTa                                                                                           >  2:2104398/1‑100 (MQ=255)
          ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGCGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTAc                                                                                          <  2:1460742/100‑1 (MQ=255)
          ttGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgc                                                                                          <  2:1476285/100‑3 (MQ=25)
          ctGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgc                                                                                          <  2:2049118/99‑3 (MQ=25)
           tGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgct                                                                                         <  2:489831/100‑4 (MQ=25)
           tGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgct                                                                                         <  1:1171512/100‑4 (MQ=25)
           tGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACg                                                                                         >  2:632938/1‑100 (MQ=255)
            gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                        <  1:1389681/100‑5 (MQ=25)
            gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                        <  1:657193/100‑5 (MQ=25)
            gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgcag                                                                                        >  1:2049132/1‑96 (MQ=25)
            gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAAAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACgg                                                                                        >  2:1850398/1‑100 (MQ=255)
               cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgttt                                                                                     <  1:694019/100‑8 (MQ=25)
               cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAg                                                                                     >  2:1908026/1‑100 (MQ=255)
               cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAg                                                                                     <  1:1089472/100‑1 (MQ=255)
               cATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAg                                                                                     >  1:1441867/1‑100 (MQ=255)
                atatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgtttg                                                                                    <  1:319670/100‑9 (MQ=25)
                 tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTa                                                                                   >  2:2199703/1‑100 (MQ=255)
                  atCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTaa                                                                                  <  1:1410087/100‑1 (MQ=255)
                  atCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTaa                                                                                  >  2:1887470/1‑100 (MQ=255)
                    cAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                <  1:1101772/100‑1 (MQ=255)
                    cAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                >  1:2317038/1‑100 (MQ=255)
                    cAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGACCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                <  2:337756/100‑1 (MQ=255)
                     aGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                <  2:49616/99‑1 (MQ=255)
                        aTCTACTTACAAAGTGACATACAGCCGCCCCGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc                                                                            >  2:39006/1‑100 (MQ=37)
                         tCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATcccc                                                                           <  1:745447/100‑1 (MQ=255)
                          cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCATCGGCACAGTCATGACAATGCCGGCACTGAAGTAATATCCCCa                                                                          <  2:1729460/100‑1 (MQ=255)
                          cTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                          <  1:1485030/100‑1 (MQ=255)
                           tAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCACCaa                                                                         >  2:1170426/1‑100 (MQ=255)
                             aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGt                                                                       >  2:562188/1‑100 (MQ=255)
                             aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGt                                                                       <  2:2040674/100‑1 (MQ=255)
                              ttACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTg                                                                      <  1:1780791/100‑1 (MQ=255)
                               tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                     >  2:2111228/1‑100 (MQ=255)
                               tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                     >  1:365781/1‑100 (MQ=255)
                               tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                     <  1:1304110/100‑1 (MQ=255)
                               tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                     >  1:1208405/1‑100 (MQ=255)
                                aCAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                                                    <  1:149796/100‑1 (MQ=255)
                                aCAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                                                    >  1:1871251/1‑100 (MQ=255)
                                 cAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAtt                                                                   >  2:519675/1‑100 (MQ=255)
                                     gTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                                               >  2:320192/1‑100 (MQ=255)
                                        aCAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCatat                                                            >  1:307396/1‑100 (MQ=255)
                                          agagagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                                          >  1:1627577/1‑100 (MQ=255)
                                           gagagCCGCCACGCGAGCGGGGCCAGAGTACTAAATAGCACCGTCACAGTAATGATAATGCCGGTACGGAAGTAATATCCCCGAGTGATTGTCATAtttt                                                         <  2:1521556/100‑1 (MQ=25)
                                           gagagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttt                                                         >  1:407079/1‑100 (MQ=255)
                                             gagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATTACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTc                                                       <  2:845935/100‑1 (MQ=255)
                                             gagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttta                                                       >  2:1046151/1‑99 (MQ=255)
                                             gagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTc                                                       >  2:242348/1‑100 (MQ=255)
                                             gagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTc                                                       <  1:1506338/100‑1 (MQ=255)
                                              agCCGCCACGCCAGCGCGGCCAGATTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTc                                                       <  2:884671/99‑1 (MQ=255)
                                              agCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCt                                                      <  2:86931/100‑1 (MQ=255)
                                               gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                     <  1:1743251/100‑1 (MQ=255)
                                                ccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTgg                                                    <  1:1511768/100‑1 (MQ=255)
                                                  gccacgccaGCGCGGCCAGAGTGACAAACAGCAGCGGCACAGTCATGACAATGCCGGTACGGAAGAAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                   <  2:1113092/99‑1 (MQ=255)
                                                  gccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGc                                                  <  1:679659/100‑1 (MQ=255)
                                                   ccacgccaGCGCGGCCAGAGTGACAAACAGCACCGTCACAGTCATTACAATGCCGGTACGTAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGc                                                  <  2:2326336/99‑1 (MQ=255)
                                                   ccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCa                                                 >  1:1887062/1‑100 (MQ=255)
                                                    cacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggtc                                                  >  2:301075/1‑96 (MQ=255)
                                                    cacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCaa                                                >  1:632537/1‑100 (MQ=255)
                                                      cgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGc                                              <  1:767430/100‑1 (MQ=255)
                                                      cgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGc                                              <  1:1785095/100‑1 (MQ=255)
                                                       gccaGCGCGGCCAGAGTGACAAATAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGca                                             >  1:1021103/1‑100 (MQ=255)
                                                       gccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGca                                             <  2:2317024/100‑1 (MQ=255)
                                                           gcgcGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                         <  2:307392/100‑1 (MQ=255)
                                                             gcgGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCa                                       <  2:257380/100‑1 (MQ=255)
                                                                gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGcc                                    >  1:2348395/1‑100 (MQ=255)
                                                                gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGcc                                    <  2:2149908/100‑1 (MQ=255)
                                                                 ccAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCca                                   >  1:127181/1‑100 (MQ=255)
                                                                        gACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagcag                             <  2:729188/99‑1 (MQ=255)
                                                                        gACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCaga                            >  1:729195/1‑100 (MQ=255)
                                                                        gACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCaga                            >  1:451230/1‑100 (MQ=255)
                                                                         aCAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCagag                           <  2:434398/100‑1 (MQ=255)
                                                                            ttccgtaCCGGCATTGTCATGACGGGGCCGGTACGGAAGTCATATCTACAAGTGATTGTCATATTTTTCTGGGCAAGCCCCAGCAGCCACCACAGAGTTg                        >  2:2331291/7‑100 (MQ=16)
                                                                              cAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                      <  1:202865/100‑1 (MQ=255)
                                                                               cgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                     >  1:191525/4‑100 (MQ=38)
                                                                               cgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTGTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                     >  1:1800888/4‑100 (MQ=37)
                                                                                gtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                    >  2:274545/3‑100 (MQ=38)
                                                                                gtaCCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                    >  1:1608214/3‑100 (MQ=255)
                                                                                  aCCGGCATTGTCATGCCTGTGCCGGTCCGGAAGAAATAACCCCAAATGATCGTCATATTTTTCTGTGCACGCCCATGCAGCCACAGCAGAGTTGcc                      >  2:2167373/1‑96 (MQ=255)
                                                                                  aCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGAc                  >  2:2203749/1‑100 (MQ=38)
                                                                                  aCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGAc                  >  1:2114409/1‑100 (MQ=38)
                                                                                   ccggcattgTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACt                 >  1:2057550/9‑100 (MQ=38)
                                                                                   ccGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACt                 <  2:2465525/100‑1 (MQ=255)
                                                                                    cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTg                <  1:2135036/100‑1 (MQ=255)
                                                                                    cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTg                <  2:787437/100‑1 (MQ=255)
                                                                                    cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTg                <  1:2135075/100‑1 (MQ=255)
                                                                                     ggCACAGTCATGACAATGCCTGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTTGGCAAGCACATGCAGCCACAGCAGGGTTGCCAGACTGc               <  2:1582430/100‑1 (MQ=255)
                                                                                     ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACttc               <  1:1730325/100‑3 (MQ=255)
                                                                                     ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGc               <  1:2151223/100‑1 (MQ=255)
                                                                                      gCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc              <  1:1957238/100‑1 (MQ=255)
                                                                                       cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCa             <  2:1761620/100‑1 (MQ=255)
                                                                                       cacaGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCa             >  1:394352/1‑100 (MQ=255)
                                                                                            tCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgg        >  1:253314/1‑100 (MQ=255)
                                                                                             cATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCggg       <  2:191523/100‑1 (MQ=255)
                                                                                             cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgg        >  2:1844359/1‑99 (MQ=255)
                                                                                              aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgggg      >  2:384197/1‑100 (MQ=255)
                                                                                                  cAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAt  >  1:687169/1‑100 (MQ=255)
                                                                                                  cAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAt  >  1:689269/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAT  >  NC_000913/2772137‑2772334

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