Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,421,437 T→C 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.9% 318.4 / NA 88I38T (ATC→ACC) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (46/40);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAA  >  NC_000913/2421340‑2421533
                                                                                                 |                                                                                                 
gCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGAccc                                                                                                 <  2:1642251/100‑1 (MQ=255)
gCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGAccc                                                                                                 <  2:1023438/100‑1 (MQ=255)
   attattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACt                                                                                              <  1:1085347/100‑1 (MQ=255)
   attattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACt                                                                                              <  1:6358/100‑1 (MQ=255)
     tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAc                                                                                            <  2:553985/100‑1 (MQ=255)
     tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAc                                                                                            >  2:320309/1‑100 (MQ=255)
     tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAc                                                                                            <  1:1980387/100‑1 (MQ=255)
       ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTAccc                                                                                         <  2:825849/100‑1 (MQ=255)
          tcttaAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgc                                                                                      <  2:1236097/97‑1 (MQ=255)
           tttaAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  1:1195529/98‑1 (MQ=255)
           tctaAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  2:1130183/98‑1 (MQ=255)
           tctaAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  1:109231/98‑1 (MQ=255)
           tctaAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  1:1133584/98‑1 (MQ=255)
           tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgcc                                                                                     <  2:139554/100‑1 (MQ=255)
             taAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCAGa                                                                                   >  1:1236100/1‑100 (MQ=255)
              aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgccaga                                                                                  >  2:109232/1‑97 (MQ=255)
              aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgccaga                                                                                  >  2:1133582/1‑97 (MQ=255)
              aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCccgccaga                                                                                  >  1:1130185/1‑97 (MQ=255)
              aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGAt                                                                                  >  2:1940877/1‑100 (MQ=255)
               aaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATa                                                                                 >  2:1071830/1‑100 (MQ=255)
                  aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAg                                                                              <  2:258410/100‑1 (MQ=255)
                  aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAg                                                                              >  1:690460/1‑100 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcg                                                                            >  2:849026/1‑100 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgcg                                                                            >  1:1346386/1‑100 (MQ=255)
                    ccTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAgc                                                                             >  1:1345714/1‑99 (MQ=255)
                      ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTAGCAATGTGTCGATCCACTACCCCGCCGATAAAgcgc                                                                           <  2:1941940/99‑1 (MQ=255)
                          tttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCAc                                                                       <  2:2147151/99‑1 (MQ=255)
                            tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTa                                                                     >  1:789218/1‑99 (MQ=255)
                              cGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                  >  1:820292/1‑100 (MQ=255)
                                tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGaa                                                                >  2:1348987/1‑100 (MQ=255)
                                 aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                               >  2:1571860/1‑100 (MQ=255)
                                 aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                               <  1:800261/100‑1 (MQ=255)
                                  cGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                              >  2:1763927/1‑100 (MQ=255)
                                  cGTTTATCGGAATCAAGGAAGAAGAAATTAAGAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                              >  1:1549431/1‑100 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                            <  2:789216/100‑1 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                            <  1:771561/100‑1 (MQ=255)
                                    tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGata                                                            >  2:1876143/1‑100 (MQ=255)
                                      tATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                           >  2:50194/1‑99 (MQ=255)
                                      tATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                          >  1:848075/1‑100 (MQ=255)
                                      tATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                          >  1:906496/1‑100 (MQ=255)
                                       aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                         >  1:334283/1‑100 (MQ=255)
                                       aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                         >  1:731644/1‑100 (MQ=255)
                                       aTCGGAAGTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                          <  2:1935391/99‑1 (MQ=255)
                                         cGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                       >  2:1387589/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCACTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  2:1279077/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  2:1535585/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  1:111984/1‑100 (MQ=255)
                                          ggAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                      >  1:115931/1‑100 (MQ=255)
                                            aaTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                    >  2:2193393/1‑100 (MQ=255)
                                             aTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                   >  2:484818/1‑100 (MQ=255)
                                                     agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                           >  1:1507402/1‑100 (MQ=255)
                                                       aagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                         >  2:789382/1‑100 (MQ=255)
                                                       aagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                         <  1:1468592/100‑1 (MQ=255)
                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                        <  2:1456814/100‑1 (MQ=255)
                                                          aaaTTAACAACCGTCATGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      <  2:1346384/100‑1 (MQ=255)
                                                          aaaTTAACAACCGTCATGATATTGTTATAACTGTGACGACCCCAACACCCGCAGATCACGCGCGCACTAGAGCAGATATCAACGCTGCGCTGAATTATCg                                      <  2:1345712/100‑1 (MQ=255)
                                                          aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      <  1:320307/100‑1 (MQ=255)
                                                          aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                      <  2:2237943/100‑1 (MQ=255)
                                                           aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                     >  1:72079/1‑100 (MQ=255)
                                                             ttAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                   <  2:1070931/100‑1 (MQ=255)
                                                             ttAACAACCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                   <  1:1348989/100‑1 (MQ=255)
                                                                acaacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                <  1:648952/100‑1 (MQ=255)
                                                                  aacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                              >  2:1451437/1‑100 (MQ=255)
                                                                  aacCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                              <  2:334285/100‑1 (MQ=255)
                                                                   acCGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                             <  1:1071832/100‑1 (MQ=255)
                                                                     cGTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                           >  1:209495/1‑100 (MQ=255)
                                                                      gTCAGGATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa                          >  2:1138459/1‑100 (MQ=255)
                                                                          ggATATTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                      <  2:111985/100‑1 (MQ=255)
                                                                          atttatTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                      <  2:1787685/97‑1 (MQ=255)
                                                                              atTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                  <  2:1507398/100‑1 (MQ=255)
                                                                               ttGTTATCAATGTGACGACCCACTACCCCGCCGATAAGGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                 >  2:1234295/1‑100 (MQ=255)
                                                                               ttGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGCAACGAAAAACattat                 >  1:2274883/1‑100 (MQ=255)
                                                                                 gTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc               >  2:209276/1‑100 (MQ=255)
                                                                                   tATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg             <  2:906494/100‑1 (MQ=255)
                                                                                   tATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg             <  1:1743127/100‑1 (MQ=255)
                                                                                   tATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg             >  2:1749602/1‑100 (MQ=255)
                                                                                    aTCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc            <  2:2274872/100‑1 (MQ=255)
                                                                                    aTCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc            >  1:1354907/1‑100 (MQ=255)
                                                                                    aTCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc            <  1:938409/100‑1 (MQ=255)
                                                                                    aTCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc            >  2:1895770/1‑100 (MQ=255)
                                                                                      cAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt          >  1:1296285/1‑100 (MQ=255)
                                                                                       aaTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg         >  2:45487/1‑100 (MQ=255)
                                                                                       aaTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg         >  2:1068398/1‑100 (MQ=255)
                                                                                        atctgACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt        <  1:2071138/97‑1 (MQ=255)
                                                                                        aTGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt        <  1:1749607/100‑1 (MQ=255)
                                                                                         tgtgACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa       <  1:1138461/100‑1 (MQ=255)
                                                                                         tgtgACGACCCACTACCCCGCCGATAAAGCACGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt        <  2:690458/99‑1 (MQ=255)
                                                                                          gtgACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag      <  1:484816/100‑1 (MQ=255)
                                                                                            gacgacCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagat    >  2:2071130/1‑99 (MQ=255)
                                                                                            gacgacCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag    >  2:1600809/1‑100 (MQ=255)
                                                                                             aCGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga   <  1:586630/100‑1 (MQ=255)
                                                                                              cGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa  >  2:1321473/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
GCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAA  >  NC_000913/2421340‑2421533

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: