Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,121,887 (GC)5→4 coding (400‑401/531 nt) hslV ← peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,8790G.96.9% 426.7 / 6.4 98coding (409/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (44/51);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0009134,121,8800C.99.0% 431.3 / ‑1.7 98coding (408/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (45/51);  minor base G (1/0);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121780‑4121972
                                                                                                   ||                                                                                            
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                               >  2:2652654/1‑100 (MQ=255)
 gATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgc                                                                                              >  2:1799820/1‑100 (MQ=255)
  aTGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcg                                                                                             <  2:2387421/100‑1 (MQ=255)
     gtgtGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgc                                                                                          >  1:3250540/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         >  1:2430731/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGACTGCAATATAAAACGCCTTTTCAGAAATTTTACGCGCGCTAAGTTCAGTGTTTGCTAACAg                                                                                               >  2:2936602/1‑94 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  2:3266513/1‑100 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  1:2211931/1‑98 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGACGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  1:1988218/1‑98 (MQ=255)
          gAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgc                                                                                     >  2:1261975/1‑97 (MQ=255)
          gAAATGGTTGGTATAGATGCAAATGTCGCCCGCGATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgc                                                                                     >  2:2143888/1‑97 (MQ=255)
            aaTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcag                                                                                   <  1:1860299/100‑6 (MQ=255)
             aTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagc                                                                                  >  2:1697986/1‑94 (MQ=255)
              tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 >  1:535774/1‑93 (MQ=255)
               ggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCt                                                                              >  1:1359351/1‑100 (MQ=255)
                  tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTggg                                                                           >  1:1289266/1‑100 (MQ=255)
                   ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          <  1:693946/100‑1 (MQ=255)
                   ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          >  1:2034411/1‑100 (MQ=255)
                   ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          >  1:1755085/1‑100 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:3157001/100‑1 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:84304/100‑1 (MQ=255)
                     tataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGt                                                                        >  1:1424169/1‑100 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       <  1:74235/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       <  1:74253/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       <  2:2069965/100‑1 (MQ=255)
                       taGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTaa                                                                      >  2:531985/1‑100 (MQ=255)
                        aGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAg                                                                     <  1:2896110/100‑1 (MQ=255)
                        aGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAg                                                                     <  1:2983718/100‑1 (MQ=255)
                          aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAggg                                                                   >  2:2799271/1‑100 (MQ=255)
                          aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAggg                                                                   <  1:1099915/100‑1 (MQ=255)
                          aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAggg                                                                   <  2:2029852/100‑1 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  <  2:2993946/100‑1 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  <  1:200688/100‑1 (MQ=255)
                              aaaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgc                                                               <  1:1256643/100‑1 (MQ=255)
                                  gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGa                                                           >  2:1481579/1‑100 (MQ=255)
                                   tCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAg                                                          >  2:1587357/1‑100 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGc                                                         >  2:683801/1‑100 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGc                                                         >  2:185308/1‑100 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  1:3005988/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCGGGGCGTAAGGGCCGCCGGAGcc                                                        >  1:1110405/1‑100 (MQ=255)
                                     gCCCGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        >  2:1146769/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  2:2884874/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  1:608742/100‑1 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  1:2725776/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  2:3192289/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    >  1:1544264/1‑100 (MQ=255)
                                          cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                   <  2:3244260/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 <  2:389073/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑ACGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 >  2:1844653/1‑100 (MQ=255)
                                             tatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCaa                                                <  2:2211905/100‑1 (MQ=255)
                                             tatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCaa                                                >  2:2485092/1‑100 (MQ=255)
                                               tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              <  2:2396728/100‑1 (MQ=255)
                                               tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              <  2:347125/100‑1 (MQ=255)
                                                 cAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGGGCGCGCCGCCGCCTGGGCGTAAGGGCCGCCGGACCCGATGGCAATAAg                                            >  2:1384381/1‑100 (MQ=255)
                                                   aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                          >  1:1871785/1‑100 (MQ=255)
                                                    cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATc                                         >  2:580116/1‑100 (MQ=255)
                                                     gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTCACA‑‑GCGCGCGCCGCAGCCTGGGCGTCAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                        >  2:1225985/1‑100 (MQ=255)
                                                      ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                       >  2:3204302/1‑100 (MQ=255)
                                                      ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                       >  2:1311185/1‑100 (MQ=255)
                                                      ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                       >  1:2384612/1‑100 (MQ=255)
                                                       cTTTTCAGCAATTTCACGTGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGACGCCGGAGACGATAGCAATAAGATCGtt                                      >  2:1619316/1‑100 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTCTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtt                                      <  1:2292487/100‑1 (MQ=255)
                                                          ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGGGCGCGGCCCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATCGCAATAAGATCGTTTTc                                   >  2:1155456/1‑100 (MQ=255)
                                                          ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTc                                   <  1:1280438/100‑1 (MQ=255)
                                                           tCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCt                                  <  2:1233521/100‑1 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 >  2:163745/1‑100 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 <  2:535766/100‑1 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 >  2:1436567/1‑100 (MQ=255)
                                                             aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTgg                                <  2:3228816/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               >  2:2604976/1‑100 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  1:2150855/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  2:2065665/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  1:124577/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               >  1:605346/1‑100 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGc                            <  2:603267/100‑1 (MQ=255)
                                                                  tttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGca                           <  1:2652687/100‑1 (MQ=255)
                                                                  tttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGca                           <  1:560186/100‑1 (MQ=255)
                                                                   ttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcac                          <  1:2434021/100‑1 (MQ=255)
                                                                    tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         <  2:3152952/100‑1 (MQ=255)
                                                                    tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         <  2:284457/100‑1 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                        >  1:3123874/1‑100 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  1:1704809/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  2:2349798/100‑1 (MQ=255)
                                                                        gggCGCTAAGTTCAGTGTTTTCTAGCA‑‑GCTCGCGCCTCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                        <  2:1254265/97‑1 (MQ=255)
                                                                        gggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGt                     <  2:1524492/100‑1 (MQ=255)
                                                                           cgcTAAGTTCAGTATTTTCTAACA‑‑GCGCGCGCCGCAGCTTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAc                  >  2:2017675/1‑100 (MQ=255)
                                                                             cTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCg                >  2:217900/1‑100 (MQ=255)
                                                                              tAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGt               >  2:2733990/1‑100 (MQ=255)
                                                                               aaGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGtt              >  2:1500451/1‑100 (MQ=255)
                                                                                  ttCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc           <  1:217905/100‑1 (MQ=255)
                                                                                     aGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt        <  1:2979877/100‑1 (MQ=255)
                                                                                       tgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtga      <  2:1424155/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     <  1:1227244/100‑1 (MQ=255)
                                                                                         ttttCTAACAGCGCGGGCCGCAGCCTGGGCGTAAGGTTCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgan     <  2:726755/99‑2 (MQ=255)
                                                                                         ttttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatg    <  2:2430703/100‑1 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAATGGCCGCCGGAGCCGATAGCAATACGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga   >  2:2286056/1‑100 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCGCCACGTCACCGGTAACGGtgacgt   <  1:2155230/100‑4 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga   <  1:59006/100‑1 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgt   <  2:2681953/100‑4 (MQ=255)
                                                                                           ttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtg  <  1:3192326/100‑5 (MQ=255)
                                                                                                   ||                                                                                            
CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121780‑4121972

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