Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,143,778 A→C E52A (GAA→GCA)  frwB → putative enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,143,7780AC93.2% 452.6 / 13.8 148E52A (GAA→GCA) frwBputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base A (6/4);  new base C (79/59);  total (85/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CAGCAGAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCAACAT  >  NC_000913/4143682‑4143874
                                                                                                |                                                                                                
cagcagAGGCGCTGGAAAGTGCCGCGAACGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGa                                                                                               >  1:1171353/1‑100 (MQ=255)
cagcagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGa                                                                                               >  1:14002/1‑100 (MQ=255)
cagcagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGa                                                                                               <  1:3132033/100‑1 (MQ=255)
cagcagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCCGTGAAGCGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGa                                                                                               >  2:2121695/1‑100 (MQ=255)
 agcagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGAt                                                                                              >  1:1329764/1‑100 (MQ=255)
   cagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGAtgt                                                                                            <  2:3748083/100‑1 (MQ=255)
   cagAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGAtgt                                                                                            >  2:3333416/1‑100 (MQ=255)
    agagGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGAtgtg                                                                                           <  2:2800671/100‑1 (MQ=255)
     gagGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTgg                                                                                          <  2:2992018/100‑1 (MQ=255)
       ggCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGACCTGACTGCGGCAGATGTGGCg                                                                                        >  2:3958954/1‑100 (MQ=255)
       ggCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCg                                                                                        <  2:1017878/100‑1 (MQ=255)
         cgcTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAg                                                                                      <  2:13811/100‑1 (MQ=255)
          gcTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAGACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAgc                                                                                     <  1:2264116/100‑1 (MQ=255)
          gcTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAgc                                                                                     >  2:3358772/1‑100 (MQ=255)
          gcTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAgc                                                                                     <  2:3955794/100‑1 (MQ=255)
           cTGGAACGTGCCGCGAAAACGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGGCTGGAAAATGACCTGACTGCCGCAGATGCGGCCAgcg                                                                                    >  1:3405593/1‑100 (MQ=255)
           cTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAgcg                                                                                    >  2:1269184/1‑100 (MQ=255)
             ggAAAGTGCCGCGAAAACGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATCGGTCTGGAAAATGACCTGACTCCGGCAGCTGTGTCGaaccct                                                                                  >  2:2259856/1‑95 (MQ=255)
                 aGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGata                                                                              <  1:311514/100‑1 (MQ=255)
                   tGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATg                                                                            >  2:3118961/1‑100 (MQ=255)
                    gCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATgg                                                                           <  2:2911364/100‑1 (MQ=255)
                     ccGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGt                                                                          <  1:1772441/100‑1 (MQ=255)
                      cgcgAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTg                                                                         >  1:2834698/1‑100 (MQ=255)
                      cgcgAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTg                                                                         >  1:1266955/1‑100 (MQ=255)
                       gcgaaagcgaaaGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGa                                                                        <  1:3932127/100‑1 (MQ=255)
                        cgaaagcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGAt                                                                       <  2:824604/100‑1 (MQ=255)
                         gaaagcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGAtt                                                                      >  1:2887669/1‑100 (MQ=255)
                          aaagcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTc                                                                     <  1:716483/100‑1 (MQ=255)
                          aaagcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTc                                                                     <  2:3526804/100‑1 (MQ=255)
                          aaagcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTc                                                                     <  2:1334669/100‑1 (MQ=255)
                            agcgaaagGCTGGGAAGTGAAAGTAGAAACCCGGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTg                                                                   <  2:2544250/100‑1 (MQ=255)
                            agcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTg                                                                   >  1:179378/1‑100 (MQ=255)
                            agcgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGCTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTg                                                                   <  1:3849076/100‑1 (MQ=255)
                              cgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGAc                                                                 <  1:2993344/100‑1 (MQ=255)
                              cgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGAc                                                                 <  1:3943310/100‑1 (MQ=255)
                              cgaaagGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGAc                                                                 <  2:534235/100‑1 (MQ=255)
                                   ggCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAg                                                            >  2:3813291/1‑100 (MQ=255)
                                   ggCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAg                                                            >  2:1612867/1‑100 (MQ=255)
                                   ggCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAg                                                            <  2:811642/100‑1 (MQ=255)
                                   ggCTGGGAAGAGAAAGAAGAAACCCAGGGGTCAATTGATCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAg                                                            <  2:440710/100‑1 (MQ=255)
                                      tGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGata                                                         >  2:2466293/1‑100 (MQ=255)
                                       gggAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGatat                                                        >  1:1388237/1‑100 (MQ=255)
                                        ggAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATc                                                       >  2:2713097/1‑100 (MQ=255)
                                        ggAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATc                                                       <  2:3671742/100‑1 (MQ=255)
                                        ggAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATc                                                       >  1:740559/1‑100 (MQ=255)
                                         gAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCg                                                      <  2:3457397/100‑1 (MQ=255)
                                          aaGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCgg                                                     <  2:276905/100‑1 (MQ=255)
                                           aGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGc                                                    <  1:1490365/100‑1 (MQ=255)
                                            gTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCa                                                   >  1:116872/1‑100 (MQ=255)
                                            gTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGGGATTCTGACCAAAGATATCGGCa                                                   >  1:3342843/1‑100 (MQ=255)
                                            gTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCAAGCGCTGATATAGTGATTCTGACCAAAGATATCGGCa                                                   >  1:3799772/1‑100 (MQ=255)
                                             tGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGGGATTCTGACCAAAGATATCGGCAt                                                  >  1:2015659/1‑100 (MQ=255)
                                             tGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCAt                                                  >  2:1217562/1‑100 (MQ=255)
                                              gAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATc                                                 >  1:88826/1‑100 (MQ=255)
                                              gAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATc                                                 >  1:3841956/1‑100 (MQ=255)
                                              gAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATc                                                 >  2:558140/1‑100 (MQ=255)
                                              gAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATc                                                 >  2:1856034/1‑100 (MQ=255)
                                                aaGTAGAAACCCAGGGGTCAATGGTTCGGGAAAATGAACTGACTTCGACAGATTTGGCGATCGCTAATAAGAGGAATCTGACCAAAGAGATCGTCATCaa                                               >  2:2195417/1‑100 (MQ=255)
                                                  gTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGt                                             <  2:1422457/100‑1 (MQ=255)
                                                   tAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGtt                                            >  1:2515702/1‑100 (MQ=255)
                                                   tAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGtt                                            >  2:1981057/1‑100 (MQ=255)
                                                    aGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGttt                                           >  1:2423318/1‑100 (MQ=255)
                                                     gAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTg                                          >  2:2112147/1‑100 (MQ=255)
                                                      aaaCCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGa                                         >  1:3419097/1‑100 (MQ=255)
                                                        aCCCCGGGGGCACTTGGTCTGGAACATGCACTAACTGCGGAAGATGTGGCGAGCGGTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAag                                       >  2:894622/1‑100 (MQ=255)
                                                          ccAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAagag                                     >  1:272813/1‑100 (MQ=255)
                                                           cAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGc                                    <  2:1572086/100‑1 (MQ=255)
                                                            aGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCg                                   <  1:1300392/100‑1 (MQ=255)
                                                            aGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCg                                   >  1:3933043/1‑100 (MQ=255)
                                                               ggTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAagagatcgg                                >  2:1719973/1‑95 (MQ=255)
                                                               ggTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGttt                                >  2:2399032/1‑100 (MQ=255)
                                                                gTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGtttt                               <  1:2844013/100‑1 (MQ=255)
                                                                gTCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGATGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGtttt                               >  2:2693124/1‑100 (MQ=255)
                                                                 tCAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTg                              <  1:990100/100‑1 (MQ=255)
                                                                 tCAATTGGGCTGGAAAATGAACTGACTGCGGCAGATGTGGCGGGCGATGATATGGGGATTCTGTCCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTg                              <  1:990086/100‑1 (MQ=255)
                                                                  cAATTGGTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGc                             <  1:2397622/100‑1 (MQ=255)
                                                                  cAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGc                             >  2:465687/1‑100 (MQ=255)
                                                                      tGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGc                         <  1:1240829/100‑1 (MQ=255)
                                                                       ggTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCa                        >  1:1037352/1‑100 (MQ=255)
                                                                       ggTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCa                        >  1:1620864/1‑100 (MQ=255)
                                                                       ggTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCa                        <  2:1015753/100‑1 (MQ=255)
                                                                        gTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCaa                       >  2:3448884/1‑100 (MQ=255)
                                                                        gTCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCaa                       >  2:1423600/1‑100 (MQ=255)
                                                                         tgtGGAAAATGAACTGGCTGCGGCAGAGGTGGCGAGCGCTGATAGGGTGATTCTGTCAGAGGATATCGGCATGAAGTTTGGAGGGCGTTTGGCGGGCAAg                      <  2:3091493/98‑1 (MQ=255)
                                                                         tCTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGTTGATTCTGACCAAGGATATCGGCCTCAAGTTTGAAGAGCGTTTTGCTGGCAAg                      <  2:19631/100‑1 (MQ=255)
                                                                          cTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGa                     >  1:9142/1‑100 (MQ=255)
                                                                          cTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGa                     >  1:1933766/1‑100 (MQ=255)
                                                                          cTGGAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCAGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGa                     >  1:3971834/1‑100 (MQ=255)
                                                                            ggAAAATGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGAcc                   <  1:1640323/100‑1 (MQ=255)
                                                                                 aTGAACTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGt              >  1:4073686/1‑100 (MQ=255)
                                                                                    aaCTGACTGCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTgcg           >  2:2729323/1‑100 (MQ=255)
                                                                                       tgactgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGt        >  1:321142/1‑100 (MQ=255)
                                                                                        gactgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTc       >  1:1483675/1‑100 (MQ=255)
                                                                                        gactgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTc       >  1:3761692/1‑100 (MQ=255)
                                                                                        gactgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTc       <  1:1443322/100‑1 (MQ=255)
                                                                                        gactgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTc       >  2:683964/1‑100 (MQ=255)
                                                                                         actgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCa      <  1:1477119/100‑1 (MQ=255)
                                                                                          ctgCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCaa     <  2:1078371/100‑1 (MQ=255)
                                                                                            gCGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCAACa   <  2:2386126/100‑1 (MQ=255)
                                                                                             cGGCAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCAACAt  >  1:4003398/1‑100 (MQ=255)
                                                                                                |                                                                                                
CAGCAGAGGCGCTGGAAAGTGCCGCGAAAGCGAAAGGCTGGGAAGTGAAAGTAGAAACCCAGGGGTCAATTGGTCTGGAAAATGAACTGACTGCGGAAGATGTGGCGAGCGCTGATATGGTGATTCTGACCAAAGATATCGGCATCAAGTTTGAAGAGCGTTTTGCTGGCAAGACCATCGTGCGCGTCAACAT  >  NC_000913/4143682‑4143874

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