Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,421,437 T→C 100% I38T (ATC→ACC)  folX → D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,4370TC98.3% 408.5 / ‑2.9 116I38T (ATC→ACC) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (68/46);  total (68/48)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 120 total aligned reads displayed.

GCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAA  >  NC_000913/2421340‑2421533
                                                                                                 |                                                                                                
  tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCAc                                                                                              >  1:1419640/1‑100 (MQ=255)
  tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCAc                                                                                              <  2:221245/100‑1 (MQ=255)
     tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAc                                                                                           >  1:71742/1‑100 (MQ=255)
       ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAccc                                                                                         >  2:897312/1‑100 (MQ=255)
       ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAccc                                                                                         <  2:3412158/100‑1 (MQ=255)
        tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAcccc                                                                                        >  2:1438571/1‑100 (MQ=255)
         cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccg                                                                                       <  1:365760/100‑1 (MQ=255)
         cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccg                                                                                       <  2:3144237/100‑1 (MQ=255)
         cGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccg                                                                                       >  2:2674822/1‑100 (MQ=255)
          gTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgc                                                                                      >  1:3295514/1‑100 (MQ=255)
           tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  1:3468033/100‑1 (MQ=255)
           tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  1:591905/100‑1 (MQ=255)
           tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     >  1:965899/1‑100 (MQ=255)
            ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgccg                                                                                    <  2:1741908/100‑1 (MQ=255)
            ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgccg                                                                                    <  2:3544337/100‑1 (MQ=255)
             taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGa                                                                                   >  1:1991490/1‑100 (MQ=255)
             taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGa                                                                                   >  1:3196516/1‑100 (MQ=255)
             taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGa                                                                                   >  1:209626/1‑100 (MQ=255)
             taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGa                                                                                   >  2:2481412/1‑100 (MQ=255)
               aaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATa                                                                                 >  2:3291454/1‑100 (MQ=255)
               aaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATa                                                                                 >  2:805774/1‑100 (MQ=255)
                aGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATaa                                                                                <  1:532410/100‑1 (MQ=255)
                 gAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATaaa                                                                               <  1:2183046/100‑1 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcg                                                                            >  1:3197857/1‑100 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcg                                                                            <  2:209626/100‑1 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcg                                                                            >  1:774347/1‑100 (MQ=255)
                       tCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcgcgc                                                                         <  1:2008922/100‑1 (MQ=255)
                       tCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcgcgc                                                                         <  1:1111762/100‑1 (MQ=255)
                        cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCa                                                                        >  1:1298611/1‑100 (MQ=255)
                         gTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCAc                                                                       >  2:1453381/1‑100 (MQ=255)
                            tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCGCTACCCCGCCGATAAAGCGCGCACTAg                                                                    >  1:3083127/1‑100 (MQ=255)
                             gCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                   <  1:2839678/100‑1 (MQ=255)
                              cGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACACCCGTCATGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATACAGCGCGCACTAGCg                                                                  >  2:951325/1‑100 (MQ=255)
                              cGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                  >  2:2715460/1‑100 (MQ=255)
                              cGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                  >  1:2511185/1‑100 (MQ=255)
                               gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                 >  1:171119/1‑100 (MQ=255)
                                tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGaa                                                                <  1:2914107/100‑1 (MQ=255)
                                 aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                               <  1:2668655/100‑1 (MQ=255)
                                 aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                               <  1:2481428/100‑1 (MQ=255)
                                  cGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                               >  1:147049/1‑99 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGGAATTAACAATCGTCAGGATATAATTATCAATGTGACAACAAACTACCCCGCCGGTAAAGCGCGCACTAGCGAAGaga                                                            >  2:3645602/1‑98 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                            >  2:1085223/1‑100 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                            <  1:448204/100‑1 (MQ=255)
                                     ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                           >  2:1177634/1‑100 (MQ=255)
                                       aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                         >  1:705220/1‑100 (MQ=255)
                                         cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                       >  1:10656/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  2:3637387/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  1:1060027/1‑100 (MQ=255)
                                           gAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                     >  2:2690163/1‑100 (MQ=255)
                                            aaTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                    >  2:2390669/1‑100 (MQ=255)
                                             aTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                   <  1:160239/100‑1 (MQ=255)
                                               tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                 >  1:56837/1‑100 (MQ=255)
                                               tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                 >  1:2024869/1‑100 (MQ=255)
                                                aaGGAAGAAGAAATTAACCACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcgc                                                <  2:1458822/100‑1 (MQ=255)
                                                   gaagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                             >  2:64652/1‑100 (MQ=255)
                                                    aagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                            >  1:3645763/1‑100 (MQ=255)
                                                    aagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                            >  1:2927267/1‑100 (MQ=255)
                                                    aagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                            <  1:2973712/100‑1 (MQ=255)
                                                     agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                           <  1:1068031/100‑1 (MQ=255)
                                                     agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                           >  1:3044256/1‑100 (MQ=255)
                                                     agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                           >  1:3465028/1‑100 (MQ=255)
                                                     agaagaaAGTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCCACGATGCGCTGAAt                                           >  2:3004179/1‑100 (MQ=255)
                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                        <  2:1651145/100‑1 (MQ=255)
                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                        >  2:652968/1‑100 (MQ=255)
                                                          aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      >  2:1925684/1‑100 (MQ=255)
                                                          aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      <  2:1998039/100‑1 (MQ=255)
                                                           aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                     >  1:240571/1‑100 (MQ=255)
                                                           aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                     <  2:774342/100‑1 (MQ=255)
                                                           aaTGAACAACCGTCAGGATATTGTTATCAAGGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      <  2:2518097/99‑1 (MQ=255)
                                                             ttAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                   <  1:1453389/100‑1 (MQ=255)
                                                               aacaacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCATCGATGCGCTGAATTATCGCACCg                                 <  2:1018654/100‑1 (MQ=255)
                                                                 caacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                               >  1:1105554/1‑100 (MQ=255)
                                                                 caacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                               >  1:74483/1‑100 (MQ=255)
                                                                 caacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                               >  2:1849739/1‑100 (MQ=255)
                                                                    ccGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGACGATAAAGAGCGCACTAGCGAAGATATAAACGATGCGCTGAATTATCGCACCGTAAc                             <  2:1304343/99‑1 (MQ=255)
                                                                       tCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                         <  2:705217/100‑1 (MQ=255)
                                                                         aGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                       >  2:2933715/1‑100 (MQ=255)
                                                                          ggATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                      >  1:201086/1‑100 (MQ=255)
                                                                          ggATATTGTTATCAATGTAACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                      <  2:2768963/100‑1 (MQ=255)
                                                                            atatTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                     <  1:1701774/99‑1 (MQ=255)
                                                                             tatTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                   >  2:1862554/1‑100 (MQ=255)
                                                                               ttGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                 >  1:2052233/1‑100 (MQ=255)
                                                                                tGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                >  1:2708413/1‑100 (MQ=255)
                                                                                tGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                >  1:2251046/1‑100 (MQ=255)
                                                                                tGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                 >  2:2035578/1‑99 (MQ=255)
                                                                                 gTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc               <  2:2052221/100‑1 (MQ=255)
                                                                                   tATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg             >  2:754888/1‑100 (MQ=255)
                                                                                    aTCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc            <  2:1171950/100‑1 (MQ=255)
                                                                                     tCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa           <  1:2654077/100‑1 (MQ=255)
                                                                                      cAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt          >  1:3120943/1‑100 (MQ=255)
                                                                                      cAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt          <  1:833005/100‑1 (MQ=255)
                                                                                      cAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt          <  1:833359/100‑1 (MQ=255)
                                                                                       aaTGTGACGACCCACTACCCCGCCGATAGAGCGCGCACTAGCGAAGATATCAACGATACGATGAATTCTCGCACCGTAACGAAAAACATTATTCAGCATg         >  2:2556518/1‑100 (MQ=255)
                                                                                        aTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt        <  1:3627862/100‑1 (MQ=255)
                                                                                        aTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt        >  1:1928671/1‑100 (MQ=255)
                                                                                         tgtgACGACCCAGTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa       >  1:2402564/1‑100 (MQ=255)
                                                                                           ngacgacCCACTACCCCGCCGATAAAGCGCGCACTAGCAAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga     >  2:2238878/2‑100 (MQ=255)
                                                                                            gacgacCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag    >  1:1366752/1‑100 (MQ=255)
                                                                                            gacgacCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag    <  2:2312187/100‑1 (MQ=255)
                                                                                              cGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa  >  1:2257552/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAA  >  NC_000913/2421340‑2421533

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