Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,930,097 (A)7→6 intergenic (+14/+42) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/2‑keto‑3‑deoxygluconate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,930,0910A.98.9% 417.5 / ‑2.3 95intergenic (+8/+48)purT/edaphosphoribosylglycinamide formyltransferase 2/2‑keto‑3‑deoxygluconate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (39/54);  minor base C (1/0);  total (40/55)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAG  >  NC_000913/1929999‑1930179
                                                                                            |                                                                                        
ggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCa                                                                                          >  1:62306/1‑93 (MQ=255)
ggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCa                                                                                          <  1:87399/93‑1 (MQ=255)
 gTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaa                                                                                         >  2:4090774/1‑93 (MQ=255)
 gTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaa                                                                                         <  2:2950145/93‑1 (MQ=255)
  tGGCACTGGCTACTGCAGAGAGCGTTGTTGACCCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAAAACAGGGTAAACCCGGGCaaa                                                                                        >  2:3606000/1‑93 (MQ=255)
   ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaa                                                                                       >  1:1152514/1‑93 (MQ=255)
   ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaa                                                                                       <  1:428108/93‑1 (MQ=255)
     cACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaa                                                                                     <  1:1079875/93‑1 (MQ=255)
     cACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaa                                                                                     <  1:4472388/93‑1 (MQ=255)
      aCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaac                                                                                    <  1:4272287/93‑2 (MQ=255)
       cTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaacg                                                                                   <  1:1193332/93‑3 (MQ=255)
        tGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaacgc                                                                                  <  2:471563/93‑4 (MQ=255)
        tGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaacgc                                                                                  <  2:3432543/93‑4 (MQ=255)
         ggCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaacgct                                                                                 >  2:4514943/1‑89 (MQ=255)
          gCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTa                                                                               >  2:1933779/1‑93 (MQ=255)
           cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTAc                                                                              <  2:1168063/93‑1 (MQ=255)
           cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTAc                                                                              <  2:1610553/93‑1 (MQ=255)
             aCTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaa                                                                            <  1:1191025/93‑1 (MQ=255)
              cTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaa                                                                           <  1:4625815/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                         <  1:2828604/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                         >  1:1046735/1‑93 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                         >  2:4331600/1‑93 (MQ=255)
                  agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATg                                                                       >  1:1088734/1‑93 (MQ=255)
                   gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGc                                                                      <  1:2705432/93‑1 (MQ=255)
                   gagagCTTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAGCCCGGCAAAAAACGCTACAAAAATGc                                                                      <  2:1419684/93‑1 (MQ=255)
                    agagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                     <  2:3680533/93‑1 (MQ=255)
                    agagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                     <  1:1357465/93‑1 (MQ=255)
                    agagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                     >  1:3529679/1‑93 (MQ=255)
                       gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGa                                                                  <  1:4008221/93‑1 (MQ=255)
                       gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGa                                                                  >  2:3112211/1‑93 (MQ=255)
                       gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGa                                                                  <  1:2769568/93‑1 (MQ=255)
                        tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACGGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGAt                                                                 >  1:308403/1‑93 (MQ=255)
                        tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGac                                                                 >  1:308402/1‑92 (MQ=255)
                        tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGAt                                                                 <  2:3890113/93‑1 (MQ=255)
                             ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCg                                                            <  1:2488420/93‑1 (MQ=255)
                             ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCg                                                            <  1:1322907/93‑1 (MQ=255)
                               gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGAt                                                          <  2:609241/93‑1 (MQ=255)
                                aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATc                                                         >  1:3093063/1‑93 (MQ=255)
                                 cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCg                                                        <  2:2064583/93‑1 (MQ=255)
                                 cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCg                                                        <  2:2341532/93‑1 (MQ=255)
                                   ccATTGAACGCGCGAAACACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                                                      >  2:4508118/1‑93 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                    <  1:3759032/93‑1 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                    >  2:147163/1‑93 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                    >  2:726149/1‑93 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                    <  2:3700203/93‑1 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                    >  1:4280196/1‑93 (MQ=255)
                                      ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                   <  1:1548301/93‑1 (MQ=255)
                                      ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                   >  2:2459855/1‑93 (MQ=255)
                                       tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                  <  2:1655301/93‑1 (MQ=255)
                                          aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttttt                                               >  2:2639672/1‑92 (MQ=255)
                                            gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                             <  2:4651323/93‑1 (MQ=255)
                                            gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                             <  1:3168315/93‑1 (MQ=255)
                                                gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAACCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                         <  2:1704805/93‑1 (MQ=255)
                                                gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                         <  2:3863915/93‑1 (MQ=255)
                                                gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                         >  2:612550/1‑93 (MQ=255)
                                                 aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAACGCTACAAAACTGCCCGATCCTCGATCCGGCATTTTGACTTTTa                                        >  1:126689/1‑93 (MQ=255)
                                                  aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                       <  1:4049477/93‑1 (MQ=255)
                                                   gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTTATTTTTACa                                      >  2:2717600/1‑93 (MQ=255)
                                                   gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                      >  2:2717599/1‑93 (MQ=255)
                                                      cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCt                                   <  1:3271415/93‑1 (MQ=255)
                                                      cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCt                                   <  1:2112657/93‑1 (MQ=255)
                                                        ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGCGCATTTTGACTTTTACAGCTTa                                 <  1:820521/93‑1 (MQ=255)
                                                           ccgGACAGGTAAAAGGACAGGGTAAAACCCGGCAAAAAACGCTACAAAAATCCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                              <  2:47403/93‑1 (MQ=255)
                                                             ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                            >  1:2093380/1‑93 (MQ=255)
                                                             ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                            >  1:371070/1‑93 (MQ=255)
                                                             ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                            >  1:2909931/1‑93 (MQ=255)
                                                              gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                           >  1:1913890/1‑93 (MQ=255)
                                                               aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTTACTTTTACAGCTTAGCGccct                          >  1:2601777/1‑91 (MQ=255)
                                                                cAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAACGCTACAAAAAGGCCCGACCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                         <  2:1809423/93‑1 (MQ=255)
                                                                cAGGTAAAAGTACAGGGTTAAAGCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                         >  1:842628/1‑93 (MQ=255)
                                                                 aGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                        >  1:2847276/1‑93 (MQ=255)
                                                                  ggTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                       >  1:3193354/1‑93 (MQ=255)
                                                                  ggTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                       >  2:1291695/1‑93 (MQ=255)
                                                                   gTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                      <  2:1454938/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGGCCCATCCTCGATCAGGCATTTTGACTTTTACATCTTATCGACTTCTACa                     >  2:1247126/1‑93 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                     <  1:663048/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                     <  1:4112731/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                     <  1:309562/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACGTTTACAGCTTAGCGCCTTCTACa                     <  2:4636805/93‑1 (MQ=255)
                                                                       aaGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                  >  1:4012447/1‑93 (MQ=255)
                                                                        aGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                 >  1:570165/1‑93 (MQ=255)
                                                                        aGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                 >  1:1062635/1‑93 (MQ=255)
                                                                         gTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                <  1:149765/93‑1 (MQ=255)
                                                                         gTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                <  1:2591474/93‑1 (MQ=255)
                                                                         gTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                >  1:97458/1‑93 (MQ=255)
                                                                          tACAGGGTTAACCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa               <  2:2893272/93‑1 (MQ=255)
                                                                          tACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa               <  2:2893269/93‑1 (MQ=255)
                                                                           aCAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc              <  2:728581/93‑1 (MQ=255)
                                                                           aCAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc              <  2:3829210/93‑1 (MQ=255)
                                                                            cAGGGTTAAACCCGGCAAAAAACGCTCAAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcg             <  1:3524915/93‑1 (MQ=255)
                                                                              gggTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTTACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg           >  2:2762039/1‑93 (MQ=255)
                                                                              gggTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGAGCCTTCTACAGCTTCActcg           >  2:2467025/1‑93 (MQ=255)
                                                                                  tAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCTGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg       <  1:1634413/93‑1 (MQ=255)
                                                                                  tAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg       <  1:1634552/93‑1 (MQ=255)
                                                                                   aaaCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc      <  1:365414/93‑1 (MQ=255)
                                                                                   aaaCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc      >  2:3099083/1‑93 (MQ=255)
                                                                                    aaCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt     >  1:4644138/1‑93 (MQ=255)
                                                                                      cccGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa   <  1:1027371/93‑1 (MQ=255)
                                                                                      cccGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa   >  1:3286763/1‑93 (MQ=255)
                                                                                       ccGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAg  >  1:3883227/1‑93 (MQ=255)
                                                                                            |                                                                                        
GGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAG  >  NC_000913/1929999‑1930179

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