Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,659 G→A 6.0% A438T (GCC→ACC)  nrdE → ribonucleoside‑diphosphate reductase 2 subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,802,6590GA 6.0% 260.7 / 4.5 83A438T (GCC→ACC) nrdEribonucleoside‑diphosphate reductase 2 subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base A (3/2);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGGTT  >  NC_000913/2802511‑2802806
                                                                                                                                                    |                                                                                                                                                   
cTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGc                                                                                                                                                    >  2:1356401/1‑150 (MQ=255)
  cAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCg                                                                                                                                                  <  1:311215/150‑1 (MQ=255)
  cAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCg                                                                                                                                                  >  1:1348853/1‑150 (MQ=255)
    gAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                                >  1:329061/1‑150 (MQ=255)
    gAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                                >  2:1155236/1‑150 (MQ=255)
      aaTTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTgcg                                                                                                                                              >  2:1023213/1‑150 (MQ=255)
       aTTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                             <  1:1458017/150‑2 (MQ=255)
       aTTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGATCGTTGAATATTGCGAACACAATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTgcgc                                                                                                                                             >  2:437083/1‑150 (MQ=255)
        ttttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcggg                                                                                                                                            <  1:363848/150‑3 (MQ=255)
         tttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCgg                                                                                                                                           <  1:375760/150‑1 (MQ=255)
          ttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGt                                                                                                                                          <  2:451664/150‑1 (MQ=255)
          ttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCATTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGt                                                                                                                                          <  1:736892/150‑1 (MQ=255)
           tGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCAGGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCCCCCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGtt                                                                                                                                         >  1:1330369/1‑150 (MQ=255)
                gTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACg                                                                                                                                    <  2:631493/150‑1 (MQ=255)
                   aaCAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCa                                                                                                                                 <  1:3298/150‑1 (MQ=255)
                    aCAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAg                                                                                                                                <  1:1495510/150‑1 (MQ=255)
                     cAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                                >  2:123093/1‑133 (MQ=255)
                     cAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                                <  1:123093/133‑1 (MQ=255)
                      aGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTa                                                                                                                              >  1:410429/1‑150 (MQ=255)
                          ccTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAg                                                                                                                          >  1:1435334/1‑150 (MQ=255)
                          ccTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAg                                                                                                                          >  2:300878/1‑150 (MQ=255)
                             cAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGata                                                                                                                       >  1:283831/1‑150 (MQ=255)
                             cAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGata                                                                                                                       >  2:1453519/1‑150 (MQ=255)
                                 ttATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAg                                                                                                                   <  1:304725/149‑1 (MQ=255)
                                  tATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGt                                                                                                                  >  1:580473/1‑150 (MQ=255)
                                   aTGCCGTGAATCTCGTCTATACCCGCACAGGCCATCATATTTCCTCCAATTTAGGTTCGTTGAATATTGCGCACTCCATGGATTCCCCCGATTTTGCGCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTc                                                                                                                 <  2:387758/150‑1 (MQ=255)
                                   aTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGCTATGAGTc                                                                                                                 >  2:624429/1‑150 (MQ=255)
                                   aTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTc                                                                                                                 <  2:1435334/150‑1 (MQ=255)
                                    tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                                <  1:1352899/150‑1 (MQ=255)
                                     gACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCat                                                                                                               >  2:1372609/1‑150 (MQ=255)
                                      aCGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCata                                                                                                              >  1:1232548/1‑150 (MQ=255)
                                        gagaATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATc                                                                                                            <  2:851436/150‑1 (MQ=255)
                                         agaATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCATGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATcc                                                                                                           >  1:921889/1‑150 (MQ=255)
                                            aTCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCa                                                                                                        <  1:29213/150‑1 (MQ=255)
                                                  aCTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGc                                                                                                  <  2:756751/150‑1 (MQ=255)
                                                  aCTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCCGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGc                                                                                                  >  2:620298/1‑150 (MQ=255)
                                                   cTATACCCGCACAGGCCATGCTTTTTCCTGCAATTTAGGTTCGTTGAATATTGCGCGCCCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGAGATGAGTCATATCCGCCGCGTGcc                                                                                                 >  2:1307119/1‑150 (MQ=255)
                                                     ataCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGt                                                                                               >  1:801821/1‑150 (MQ=255)
                                                      taCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTc                                                                                              >  2:708/1‑150 (MQ=255)
                                                          cGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATc                                                                                          <  1:1000152/150‑1 (MQ=255)
                                                              caGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAg                                                                                      <  2:1400058/150‑1 (MQ=255)
                                                                   cATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGa                                                                                 <  1:1372609/150‑1 (MQ=255)
                                                                      gATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAAt                                                                              >  2:1439410/1‑150 (MQ=255)
                                                                      gATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAAt                                                                              >  1:1099576/1‑150 (MQ=255)
                                                                        tatTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgc                                                                            >  1:1514250/1‑150 (MQ=255)
                                                                         atTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgcc                                                                           >  1:153587/1‑150 (MQ=255)
                                                                         atTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgcc                                                                           <  2:1348853/150‑1 (MQ=255)
                                                                           ttCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgccgc                                                                         <  1:300878/150‑1 (MQ=255)
                                                                            tCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgccgcc                                                                        <  1:1023213/150‑1 (MQ=255)
                                                                             ccTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCt                                                                       <  1:831881/150‑1 (MQ=255)
                                                                                         gTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCg                                                           >  1:508119/1‑150 (MQ=255)
                                                                                            cGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGAc                                                        >  2:26283/1‑150 (MQ=255)
                                                                                             gTTGAATATTGCGCACACCATGGATTCCCCCGAATTTGCCCGCACGGTAGAGACTGCCCTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACt                                                       >  1:8984/1‑150 (MQ=255)
                                                                                               tGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgg                                                     >  2:1191165/1‑150 (MQ=255)
                                                                                               tGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgg                                                     >  1:107879/1‑150 (MQ=255)
                                                                                               tGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGAACGGTAGAGACTGGCGTGCGAGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACCCCATCGGACTgg                                                     >  1:1096526/1‑150 (MQ=255)
                                                                                                gAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCATGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGGAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTggg                                                    <  2:921889/150‑1 (MQ=255)
                                                                                                 aaTATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgggg                                                   >  2:878618/1‑150 (MQ=255)
                                                                                                  atatTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGc                                                  >  1:521760/1‑150 (MQ=255)
                                                                                                        cgcACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGa                                            >  1:317114/1‑150 (MQ=255)
                                                                                                          cacacaATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAt                                          >  2:676382/1‑150 (MQ=255)
                                                                                                          cacaCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAt                                          >  1:86209/1‑150 (MQ=255)
                                                                                                           acgccATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAtt                                         <  2:793159/147‑1 (MQ=255)
                                                                                                              ccATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTac                                      <  1:1423483/150‑1 (MQ=255)
                                                                                                                    tttcggccGATTTTGCCCGCACGGTAGAGAGTGGCGTGCGCTGTTTAACTGCAGTATCAGATATGAGTCATATACGCAGCTTGCCGTCCATCGAAGCCGGAAATGCCGGCTCGCACTCCATCGGACTGGGGCAGATGAATTTACACGGCt                                <  2:1096526/144‑1 (MQ=255)
                                                                                                                          cctaTTTTGACATCACGTTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGAAGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTgg                          >  2:48815/4‑150 (MQ=255)
                                                                                                                              ttttGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGcgcg                      <  2:329061/150‑1 (MQ=255)
                                                                                                                                 tGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                   >  1:601839/1‑150 (MQ=255)
                                                                                                                                 accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                   >  1:1358338/5‑150 (MQ=255)
                                                                                                                                  gCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGaa                  >  2:488060/1‑150 (MQ=255)
                                                                                                                                  ccgcGCACGGCAGTCTCTACCGTGAGCGGTTTAACGGCAGTATAAGATATCACTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCAGCATCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGaa                  >  2:1101712/4‑150 (MQ=255)
                                                                                                                                   cgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg                 >  1:1278034/3‑150 (MQ=255)
                                                                                                                                    ccGCACGGTAGATACTGCCGTGCGCGGTTTAACGGCATTATCATATATTAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAgg                <  2:521760/150‑1 (MQ=255)
                                                                                                                                    ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAgg                <  1:1439410/150‑1 (MQ=255)
                                                                                                                                    ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAgg                <  2:283831/150‑1 (MQ=255)
                                                                                                                                    ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAgg                <  1:710284/150‑1 (MQ=255)
                                                                                                                                     cGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGc               <  1:1453519/150‑1 (MQ=255)
                                                                                                                                     cGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGc               >  1:224160/1‑150 (MQ=255)
                                                                                                                                     cGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGc               <  2:1514250/150‑1 (MQ=255)
                                                                                                                                           gTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGGCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCt         <  1:1155236/150‑1 (MQ=255)
                                                                                                                                                gaCTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAGTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATgg    >  2:57802/1‑150 (MQ=255)
                                                                                                                                                  cTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGGtt  >  2:595540/1‑150 (MQ=255)
                                                                                                                                                  cTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGGtt  >  2:472218/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGGTT  >  NC_000913/2802511‑2802806

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