Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,654 A→T 6.3% E436V (GAG→GTG) ‡ 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,6540AT 6.3% 241.6 / 5.5 80E436V (GAG→GTG) ‡nrdEribonucleoside‑diphosphate reductase 2 subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (3/2);  total (44/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.18e-01

CTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCT  >  NC_000913/2802506‑2802799
                                                                                                                                                    |                                                                                                                                                 
ctctGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTagag                                                                                                                                                  <  1:782500/150‑1 (MQ=255)
     cTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGc                                                                                                                                             >  2:1356401/1‑150 (MQ=255)
       cAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCg                                                                                                                                           <  1:311215/150‑1 (MQ=255)
       cAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCg                                                                                                                                           >  1:1348853/1‑150 (MQ=255)
         gAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                         >  2:1155236/1‑150 (MQ=255)
         gAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                         >  1:329061/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                                                                    <  2:1348853/150‑1 (MQ=255)
                                                                              atTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgcc                                                                    >  1:153587/1‑150 (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)
                                                                                                                                                    |                                                                                                                                                 
CTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCT  >  NC_000913/2802506‑2802799

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: