Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,810,321 T→C 6.2% intergenic (+17/+22) coaD → / ← mutM pantetheine‑phosphate adenylyltransferase/formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,810,3210TC 6.2% 935.8 / 31.9 273intergenic (+17/+22)coaD/mutMpantetheine‑phosphate adenylyltransferase/formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (125/131);  new base C (9/8);  total (134/139)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 291 total aligned reads displayed.

TTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAG  >  NC_000913/3810173‑3810469
                                                                                                                                                    |                                                                                                                                                    
tttCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGt                                                                                                                                                      <  2:679042/149‑1 (MQ=255)
     gATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGcc                                                                                                                                                 <  3:206901/149‑1 (MQ=255)
     gATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGcc                                                                                                                                                 >  2:897017/1‑149 (MQ=255)
           gTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCgg                                                                                                                                           >  4:393399/1‑149 (MQ=255)
               aaaGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTa                                                                                                                                                     >  2:728040/1‑135 (MQ=255)
               aaaGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTACGCCGGATGGCATACCATCCGgcat                                                                                                                                        >  1:771894/1‑146 (MQ=255)
               aaaGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTACGCCGGATGGCATACCATCCGgcat                                                                                                                                        <  2:771894/148‑3 (MQ=255)
                    gTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGTCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATc                                                                                                                                              >  3:377360/1‑137 (MQ=255)
                     tGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                                                 <  3:393399/149‑1 (MQ=255)
                       gTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTAt                                                                                                                                                    >  2:1012967/1‑128 (MQ=255)
                          gTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACt                                                                                                                            <  3:1130533/149‑1 (MQ=255)
                          gTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCAt                                                                                                                                                    <  4:1283350/125‑1 (MQ=255)
                           tttATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTAt                                                                                                                                                    >  2:512528/1‑124 (MQ=255)
                           tttATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTAt                                                                                                                                                    <  1:512528/124‑1 (MQ=255)
                                ctctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                                                                      >  2:1448343/1‑149 (MQ=255)
                                      aTCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGt                                                                                                                                                      <  2:1309884/111‑1 (MQ=255)
                                      aTCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTa                                                                                                                              <  2:384661/135‑1 (MQ=255)
                                      aTCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTa                                                                                                                              >  1:384661/1‑135 (MQ=255)
                                        cGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGa                                                                                                              <  1:690009/149‑1 (MQ=255)
                                          ttGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                                                                          >  2:1005450/1‑135 (MQ=255)
                                          ttGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                                                                          <  1:1005450/135‑1 (MQ=255)
                                          ttGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACa                                                                                                            >  1:614557/1‑149 (MQ=255)
                                           tGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACaa                                                                                                           <  2:614557/149‑1 (MQ=255)
                                               gAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaa                                                                                                       >  1:856267/1‑149 (MQ=255)
                                                  agagGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACg                                                                                                    >  4:1424395/1‑149 (MQ=255)
                                                      gTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGcc                                                                                                                                                 <  4:369196/100‑1 (MQ=255)
                                                      gTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGc                                                                                                >  4:1050479/1‑149 (MQ=255)
                                                      gTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGc                                                                                                <  1:1448343/149‑1 (MQ=255)
                                                          cgcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGc                                                                                            >  2:729996/1‑149 (MQ=255)
                                                           gcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCt                                                                                           >  1:736649/1‑149 (MQ=255)
                                                                 tCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCAt                                                                                     >  3:29105/1‑149 (MQ=255)
                                                                  cAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATg                                                                                    <  1:315528/149‑1 (MQ=255)
                                                                   aGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCg                                                                                                                                            >  4:685430/1‑92 (MQ=255)
                                                                   aGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCg                                                                                                                                            <  3:685430/92‑1 (MQ=255)
                                                                      gCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTa                                                                                <  4:237658/149‑1 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTa                                                                                                                                                     <  3:642818/72‑1 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTa                                                                                                                                                     >  4:642818/1‑72 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTAt                                                                                                                                                    >  1:35097/1‑73 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGa                                                                                                              <  2:1187257/111‑1 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgc                                                                                         >  3:321691/1‑132 (MQ=255)
                                                                              aCCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                                                                        >  4:658966/1‑149 (MQ=255)
                                                                                ccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgc                                                                                                                                          <  3:487596/81‑1 (MQ=255)
                                                                                ccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgc                                                                                                                                          >  4:487596/1‑81 (MQ=255)
                                                                                ccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATg                                                                                                                                   <  3:347444/88‑1 (MQ=255)
                                                                                 cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAAt                                                                     >  3:899299/1‑149 (MQ=255)
                                                                                     tCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGc                                                                                                                                                  <  3:1300011/68‑1 (MQ=255)
                                                                                        tGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACtt                                                                                                                           <  2:1048803/88‑1 (MQ=255)
                                                                                         gCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc                                                             <  1:1300407/149‑1 (MQ=255)
                                                                                         gCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc                                                             <  3:1279586/149‑1 (MQ=255)
                                                                                            ggAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                                                                  <  1:1093932/125‑1 (MQ=255)
                                                                                             gagaATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca                                                         >  1:1183064/1‑149 (MQ=255)
                                                                                               gaATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGcaca                                                       <  1:728948/149‑1 (MQ=255)
                                                                                                   gTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                                                                          <  2:994687/78‑1 (MQ=255)
                                                                                                   gTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg                                                   >  1:79667/1‑149 (MQ=255)
                                                                                                    tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGGACTGCCGACaa                                                                                                           >  2:1417735/1‑92 (MQ=255)
                                                                                                    tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGc                                                                 >  4:72243/1‑134 (MQ=255)
                                                                                                    tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGc                                                                 <  3:72243/134‑1 (MQ=255)
                                                                                                    tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAccc                                                    <  3:604857/147‑1 (MQ=255)
                                                                                                     ccATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCAt                                                                                     >  4:957340/1‑113 (MQ=255)
                                                                                                      cATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGCATTACTTCTGGCACTGCCGACaa                                                                                                           <  2:583301/90‑1 (MQ=255)
                                                                                                      cATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACaa                                                                                                           >  1:583301/1‑90 (MQ=255)
                                                                                                      cATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGaca                                                >  4:377625/1‑149 (MQ=255)
                                                                                                       aTCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGtt                                                                                  <  3:1422800/114‑1 (MQ=255)
                                                                                                       aTCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg                                                           <  2:349491/137‑1 (MQ=255)
                                                                                                          aGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACaa                                                                                                           <  2:1204732/86‑1 (MQ=255)
                                                                                                              gcTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCat                                        <  3:658966/149‑2 (MQ=255)
                                                                                                               cTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAcc                                                     <  2:1331599/135‑1 (MQ=255)
                                                                                                               cTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCCGCGTACCGCACACCCGACACGGCTCAcc                                       <  4:899299/149‑1 (MQ=255)
                                                                                                                   tGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTc                                                                                                                          >  1:236714/1‑62 (MQ=255)
                                                                                                                    ggCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGt                                                               <  4:940588/120‑1 (MQ=255)
                                                                                                                    ggCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGt                                                               >  3:940588/1‑120 (MQ=255)
                                                                                                                     gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACgg                                             >  4:1166949/1‑137 (MQ=255)
                                                                                                                       gAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgc                                                                                                                                          <  2:286809/42‑1 (MQ=255)
                                                                                                                       gAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGt                                                                              <  2:783967/102‑1 (MQ=255)
                                                                                                                          gTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATc                                                                    <  2:1118439/109‑1 (MQ=255)
                                                                                                                           ttAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCtt                                     >  3:967795/1‑139 (MQ=255)
                                                                                                                            tagcgtagcgTTTATGCCGGATGGTATGCCATCCGGCGCGCATGTATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc                                                            <  1:917990/115‑1 (MQ=255)
                                                                                                                                gtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAAcc                       <  2:478922/148‑1 (MQ=255)
                                                                                                                                gtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAAcc                       >  1:478922/1‑148 (MQ=255)
                                                                                                                                  agcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGc                    >  2:411303/1‑149 (MQ=255)
                                                                                                                                     gtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAAt                 <  1:1009899/149‑1 (MQ=255)
                                                                                                                                      tttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTc                                          >  1:1200662/1‑123 (MQ=255)
                                                                                                                                         aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCt             >  3:1364346/1‑149 (MQ=255)
                                                                                                                                         aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCt             >  4:309403/1‑149 (MQ=255)
                                                                                                                                           gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGa           <  4:1199682/149‑1 (MQ=255)
                                                                                                                                           gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGa           <  3:1170843/149‑1 (MQ=255)
                                                                                                                                           gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCt                      <  4:795698/138‑1 (MQ=255)
                                                                                                                                           gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCt                      >  3:795698/1‑138 (MQ=255)
                                                                                                                                            ccGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGCCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAg          >  2:844520/1‑149 (MQ=255)
                                                                                                                                             cGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTg                                                                                                 >  4:1300947/1‑61 (MQ=255)
                                                                                                                                             cGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTc               >  1:1399871/1‑143 (MQ=255)
                                                                                                                                             cGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGa           <  3:865472/147‑1 (MQ=255)
                                                                                                                                             cGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGc         <  4:329046/149‑1 (MQ=255)
                                                                                                                                             cGGATGGTATGCCATCCGGCGCGCATGAATTACCTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGa           >  4:865472/1‑147 (MQ=255)
                                                                                                                                             cGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGc                                            >  1:1311972/1‑114 (MQ=255)
                                                                                                                                              ggATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCg        <  4:1364346/149‑1 (MQ=255)
                                                                                                                                               gATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGc                    <  1:1045888/136‑1 (MQ=255)
                                                                                                                                                  ggTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGaaa     >  4:265624/1‑148 (MQ=255)
                                                                                                                                                    tATGTCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAg  >  3:17884/1‑149 (MQ=255)
                                                                                                                                                    tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAg  >  3:257068/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAG  >  NC_000913/3810173‑3810469

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: