Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,401,844 A→G 100% I471T (ATC→ACC)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,8440AG98.6% 288.6 / ‑5.8 90I471T (ATC→ACC) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (42/44);  minor base T (2/1);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTG  >  NC_000913/1401748‑1401940
                                                                                                |                                                                                                
gacgacGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAg                                                                                               <  2:1896257/100‑1 (MQ=255)
   gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                            <  2:910502/100‑1 (MQ=255)
   gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                            <  1:101259/100‑1 (MQ=255)
   gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                            <  1:2742772/100‑1 (MQ=255)
   gacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGc                                                                                            <  1:1673978/100‑1 (MQ=255)
     cgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAg                                                                                          <  1:52774/100‑1 (MQ=255)
      ggTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGt                                                                                         >  1:16747/1‑100 (MQ=255)
       gTATGGATGACATACTCATGGGCATTTCTTGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAg                                                                                          <  1:475035/98‑1 (MQ=255)
       gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGt                                                                                         <  1:2078178/99‑1 (MQ=255)
       gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg                                                                                        >  1:1474665/1‑100 (MQ=255)
        tATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGa                                                                                       <  2:1474650/100‑1 (MQ=255)
          tGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATGTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATa                                                                                     >  2:1939353/1‑100 (MQ=255)
           ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                    >  2:1096423/1‑100 (MQ=255)
           ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc                                                                                    >  2:1175084/1‑100 (MQ=255)
             aTGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGAGTCACGTTTTTCGGGGTAGGGCAGTGATACgg                                                                                  >  1:842211/1‑100 (MQ=255)
               gACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGtt                                                                                >  2:1600908/1‑100 (MQ=255)
                aCATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTg                                                                               >  2:2153084/1‑100 (MQ=255)
                 cATACTCATGGGCATTCCTCGGTTGTTGTTTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGc                                                                              <  1:2318695/100‑1 (MQ=255)
                  aTACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTg                                                                               <  2:2279370/98‑1 (MQ=255)
                   tACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgt                                                                            <  2:196599/100‑1 (MQ=255)
                   tACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgt                                                                            >  2:44755/1‑100 (MQ=255)
                     cTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgtgt                                                                          >  1:2179732/1‑100 (MQ=255)
                     cTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgtgt                                                                          >  2:1931848/1‑100 (MQ=255)
                      tCATGGGTATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGTGCAGTGATACGGTTGCGTGTc                                                                         <  2:479268/100‑1 (MQ=255)
                       cATGGGCATTCGTCGGTTGTTGTGTTTGCGTTGGTGTTGTTATTTTAAAGGTGACGGTGTTACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgtgtgc                                                                        >  2:1349367/1‑98 (MQ=255)
                       cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                        >  2:1673943/1‑100 (MQ=255)
                       cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc                                                                        >  1:323985/1‑100 (MQ=255)
                            gCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCATGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTAc                                                                   <  2:48705/100‑1 (MQ=255)
                              aTTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACtt                                                                 <  1:1516155/100‑1 (MQ=255)
                               ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg                                                                >  1:2317695/1‑100 (MQ=255)
                                 ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt                                                              >  2:144284/1‑100 (MQ=255)
                                  cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg                                                             <  2:1475146/100‑1 (MQ=255)
                                  cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg                                                             >  1:1044512/1‑100 (MQ=255)
                                   tCGGTCTTTGTGTTTGCGTTGTTGTTGTTATTTTTATGGTGACTGTGTCACGTTTTTCGGTGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTgn                                                            <  2:1373974/100‑2 (MQ=255)
                                    cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTTGTTTTAGTGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt                                                           >  2:1391328/1‑100 (MQ=255)
                                           tgtgTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttg                                                    >  1:686001/1‑100 (MQ=255)
                                             tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgtt                                                  >  1:2163189/1‑100 (MQ=255)
                                             tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgtt                                                  >  1:2075669/1‑100 (MQ=255)
                                              gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCTGGGTAGGGCAGTGATACGGTTGAGTGTTCGTTACTTGCTGATGCTgttgtt                                                  >  2:3032227/1‑99 (MQ=255)
                                               ttttCGTTGTTTTTGTTATTTTAAAGGTTTCGGTGTCACTTTTTTCTGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTgg                                                <  1:124270/100‑1 (MQ=255)
                                               tttGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGTTTGCGTGTCTGTTACTTGCTGATGCTCttgttgn                                                <  2:616931/100‑2 (MQ=255)
                                                ttGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGc                                               <  2:2986348/100‑1 (MQ=255)
                                                 tGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCa                                              >  2:1663889/1‑100 (MQ=255)
                                                  gCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCAt                                             <  2:2559168/100‑1 (MQ=255)
                                                     ttgttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTTTTTTTTTTTCAGTGATACGGTTGCGTGTCCGTTACTTGCTGGTGCTCTTGTTGGCATTCt                                          >  2:1241467/1‑100 (MQ=255)
                                                     ttgttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCt                                          <  1:2841459/100‑1 (MQ=255)
                                                        ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTTGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                       <  1:704702/100‑1 (MQ=255)
                                                        ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                       <  2:771319/100‑1 (MQ=255)
                                                        ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc                                       <  2:2362879/100‑1 (MQ=255)
                                                         ttttgttAGTTTAAAGGTGACGGTGTCACGTTTTTCATGGTAGGTCAGTGATACTTTTGCGTGTCCGTTACTTGCTTATGCTCTTGTTGGCATTCTTGc                                       <  2:1727319/97‑1 (MQ=255)
                                                         tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa                                      <  1:833068/100‑1 (MQ=255)
                                                         tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa                                      <  2:174307/100‑1 (MQ=255)
                                                         tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa                                      >  2:2460261/1‑100 (MQ=255)
                                                           ttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAAt                                    >  1:2203750/1‑100 (MQ=255)
                                                            tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATa                                   >  1:496315/1‑100 (MQ=255)
                                                            tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATa                                   >  2:2242037/1‑100 (MQ=255)
                                                             gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg                                  >  1:1512418/1‑100 (MQ=255)
                                                             gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg                                  >  1:1247974/1‑100 (MQ=255)
                                                             gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGGTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg                                  >  1:2902247/1‑100 (MQ=255)
                                                              ttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGc                                 >  1:2593611/1‑100 (MQ=255)
                                                                atgttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGccc                               <  2:2528332/97‑1 (MQ=255)
                                                                  tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  2:1605869/100‑1 (MQ=255)
                                                                  tttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  1:1119316/100‑1 (MQ=255)
                                                                   ttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGTCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg                             <  2:472501/99‑1 (MQ=255)
                                                                    tAAAGGTGACGGTGTCACGTTTTTCTGGGTAGTGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGaa                           <  1:2803001/100‑1 (MQ=255)
                                                                    tAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGaa                           >  2:1696677/1‑100 (MQ=255)
                                                                     aaaGGTGACGGTGTCACGATTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCACTGAAt                          >  2:1481976/1‑100 (MQ=255)
                                                                      aaGGTGACTGTGTCACGTTTTTCGTGGTCGGGCAGTCATACTGTTGCGTGTCCGTTACTTGCTGATTCTCTTGTTGGCATTCTTGAAATAGCCCTGaa                           <  2:2622711/98‑1 (MQ=255)
                                                                      aaGGTGACGGTGTCACGTTTTTCGGGTTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATc                         <  1:829766/100‑1 (MQ=255)
                                                                      aaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATc                         <  1:1608440/100‑1 (MQ=255)
                                                                      aaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATc                         <  1:1935206/100‑1 (MQ=255)
                                                                       aGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCa                        <  1:2314472/100‑1 (MQ=255)
                                                                        ggTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCaa                       >  2:2599509/1‑100 (MQ=255)
                                                                        ggTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCaa                       >  2:611924/1‑100 (MQ=255)
                                                                        ggTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCaa                       >  2:2598777/1‑100 (MQ=255)
                                                                         gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg                      <  1:1805296/100‑1 (MQ=255)
                                                                         gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCATGAATCAAg                      >  1:2162582/1‑100 (MQ=255)
                                                                          tGACGGTGTTACGTTTTTCGGGGTATGGCAGTGATACGGTTGCGTGTCCGTTAGTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                     <  1:197857/100‑1 (MQ=255)
                                                                          tGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                     <  1:2768291/100‑1 (MQ=255)
                                                                          tGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                     <  1:735376/100‑1 (MQ=255)
                                                                          tGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAgg                     <  2:1317902/100‑1 (MQ=255)
                                                                             cGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAag                  >  2:931915/1‑100 (MQ=255)
                                                                                tgtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAgaga               >  1:2007626/1‑100 (MQ=255)
                                                                                   cACGTTTTTCGGGGTATGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTg            >  2:138199/1‑100 (MQ=255)
                                                                                   cACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTg            <  2:170065/100‑1 (MQ=255)
                                                                                   cACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTg            >  1:520628/1‑100 (MQ=255)
                                                                                     cGTTTTTCGGGGTAGGGCAGTGATATGGTTGCGTGTTCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGAc          <  2:1702591/100‑1 (MQ=255)
                                                                                       tttttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTACAATAGCCCTGAATCAAGGAAGAAATTGAc          >  1:1242447/1‑98 (MQ=255)
                                                                                         tttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGACATTGACtgtg      >  2:1990172/1‑100 (MQ=255)
                                                                                            cGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGAATGTGGtt   >  1:2149151/1‑100 (MQ=255)
                                                                                             ggggTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtg  >  1:2779711/1‑100 (MQ=255)
                                                                                             ggggTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCCGATGCTCTTGTTGGCATTCTTGAAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTtg  <  2:2738514/100‑1 (MQ=255)
                                                                                                |                                                                                                
GACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTG  >  NC_000913/1401748‑1401940

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