Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,275 T→A 6.6% intergenic (+24/‑172) atoC → / → atoD DNA‑binding transcriptional activator/ornithine decarboxylase inhibitor AtoC/acetyl‑CoA:acetoacetyl‑CoA transferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,323,2750TA 6.6% 281.7 / 5.8 89intergenic (+24/‑172)atoC/atoDDNA‑binding transcriptional activator/ornithine decarboxylase inhibitor AtoC/acetyl‑CoA:acetoacetyl‑CoA transferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (44/39);  new base A (4/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATCATTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAA  >  NC_000913/2323126‑2323422
                                                                                                                                                     |                                                                                                                                                   
atcatTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAt                                                                                                                                                     <  2:1245821/150‑1 (MQ=255)
   atTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTg                                                                                                                                                  <  1:133573/150‑1 (MQ=255)
    ttATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGc                                                                                                                                                 <  2:34831/150‑1 (MQ=255)
     tATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGca                                                                                                                                                <  2:604252/150‑1 (MQ=255)
      aTGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcac                                                                                                                                               <  2:219636/150‑1 (MQ=255)
           aGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                                                                          >  2:1000882/1‑150 (MQ=255)
            gTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcg                                                                                                                                         >  2:866621/1‑150 (MQ=255)
             tGCTGGGACAACAAGAAGGAAATCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                        <  1:670633/150‑1 (MQ=255)
             tGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                        >  1:1362086/1‑150 (MQ=255)
              gCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCACTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCa                                                                                                                                       >  1:1300719/1‑150 (MQ=255)
                 ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                                                    >  1:232152/1‑150 (MQ=255)
                    acaacaAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                                                                 <  1:827562/150‑1 (MQ=255)
                       acaaGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                                                              <  2:198629/150‑1 (MQ=39)
                        caaGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGAGCAATTTTCTGCa                                                                                                                             >  2:742793/1‑150 (MQ=255)
                           gaagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAg                                                                                                                          >  1:589371/1‑150 (MQ=255)
                            aagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                                                                         <  1:406804/150‑1 (MQ=255)
                              ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                                                                       >  1:785370/1‑150 (MQ=255)
                                 aaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTc                                                                                                                    <  2:256822/150‑1 (MQ=255)
                                  aCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCa                                                                                                                   >  1:1115907/1‑150 (MQ=255)
                                     gAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                                                                <  1:434131/150‑1 (MQ=255)
                                       aCCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtctc                                                                                                              <  2:267145/150‑1 (MQ=255)
                                         ccGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                                                            <  2:534712/150‑1 (MQ=255)
                                         ccGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                                                            <  2:845232/150‑1 (MQ=255)
                                          cGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGcaagtt                                                                                                                    <  1:371077/141‑6 (MQ=38)
                                          cGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGcaagtt                                                                                                                    >  2:371077/1‑136 (MQ=38)
                                             aCTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                                                                        <  2:785370/150‑1 (MQ=255)
                                              cTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaa                                                                                                       >  1:480065/1‑150 (MQ=255)
                                                    tAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcccc                                                                                                 >  2:1161882/1‑150 (MQ=255)
                                                       tGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                                                              <  2:1362086/150‑1 (MQ=255)
                                                             gCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTTTCCTTATAAATCCCCATCCAAttt                                                                                        >  2:1316630/1‑150 (MQ=255)
                                                             gCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                                                        >  1:358729/1‑150 (MQ=255)
                                                             gCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                                                        <  1:1369704/150‑1 (MQ=255)
                                                                  aGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                                                                   <  1:865247/150‑1 (MQ=255)
                                                                  aGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                                                                   <  1:816678/150‑1 (MQ=255)
                                                                      gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                               <  1:986265/150‑1 (MQ=255)
                                                                      gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                               <  2:884377/150‑1 (MQ=255)
                                                                         gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                                                            >  2:1332240/1‑150 (MQ=255)
                                                                             aTTGATGTATAAACTCCAGTAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                                                        <  1:427190/150‑1 (MQ=255)
                                                                                gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTa                                                                     >  1:260702/1‑150 (MQ=255)
                                                                                 aTGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCACTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                                                                    <  1:1000882/150‑1 (MQ=255)
                                                                                  tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCGGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                   >  2:704601/1‑150 (MQ=255)
                                                                                  tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                   >  1:32557/1‑150 (MQ=255)
                                                                                    tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                                                                 >  1:408847/1‑150 (MQ=255)
                                                                                    tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                                                                 >  1:92111/1‑150 (MQ=255)
                                                                                     ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaa                                                                >  2:876682/1‑150 (MQ=255)
                                                                                      taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                                                               >  2:1162164/1‑150 (MQ=255)
                                                                                            ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                                                                  <  2:14066/141‑1 (MQ=255)
                                                                                            ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                                                                  >  1:14066/1‑141 (MQ=255)
                                                                                               ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                                                      <  2:379744/150‑1 (MQ=255)
                                                                                               ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                                                      >  2:382401/1‑150 (MQ=255)
                                                                                                gAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                     >  1:929914/1‑150 (MQ=255)
                                                                                                 aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc                                                    <  2:545471/150‑1 (MQ=255)
                                                                                                  aTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa                                                   >  1:335632/1‑150 (MQ=255)
                                                                                                   tACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa                                                  <  2:415012/150‑1 (MQ=255)
                                                                                                     cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata                                                >  2:899543/1‑150 (MQ=255)
                                                                                                     cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata                                                >  1:550502/1‑150 (MQ=255)
                                                                                                      ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCATTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCGCTTCCAATTTATCCCTTGATATTGAATTAGTTAAATAACTAAATCCaataa                                               <  2:1300719/150‑1 (MQ=255)
                                                                                                        tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                             <  1:38852/150‑1 (MQ=255)
                                                                                                           cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCa                                          <  2:232152/150‑1 (MQ=255)
                                                                                                           cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCa                                          <  2:32936/150‑1 (MQ=255)
                                                                                                            gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAt                                         >  1:1063359/1‑150 (MQ=255)
                                                                                                             aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtt                                        >  1:970721/1‑150 (MQ=255)
                                                                                                             aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtt                                        >  2:1377250/1‑150 (MQ=255)
                                                                                                                cggcggATGTATAACACCAAAACTTGCTATTCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTg                                     <  2:849749/150‑1 (MQ=255)
                                                                                                                cggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTg                                     >  1:1361004/1‑150 (MQ=255)
                                                                                                                  gcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGc                                   <  1:1162164/150‑1 (MQ=255)
                                                                                                                    ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTTTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCAc                                 <  1:1316630/150‑1 (MQ=255)
                                                                                                                    ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCAc                                 >  2:1137710/1‑150 (MQ=255)
                                                                                                                       tGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcc                              <  2:480065/150‑1 (MQ=255)
                                                                                                                         tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                            >  1:1290442/1‑150 (MQ=255)
                                                                                                                          ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCt                           <  1:357318/150‑1 (MQ=255)
                                                                                                                          ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCt                           <  1:1325115/150‑1 (MQ=255)
                                                                                                                           taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCtt                          >  1:537092/1‑150 (MQ=255)
                                                                                                                            aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTg                         >  2:502800/1‑150 (MQ=255)
                                                                                                                              cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCt                       >  2:1086200/1‑150 (MQ=255)
                                                                                                                               acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                             <  1:45652/127‑1 (MQ=255)
                                                                                                                               acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                             >  2:45652/1‑127 (MQ=255)
                                                                                                                                ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTAt                     >  1:493484/1‑150 (MQ=255)
                                                                                                                                ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTAt                     >  1:514444/1‑150 (MQ=255)
                                                                                                                                ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTAt                     >  2:1388777/1‑150 (MQ=255)
                                                                                                                                ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATACCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGAACTACCCTTGCTAt                     >  1:1273774/1‑150 (MQ=255)
                                                                                                                                   aaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGc                  <  1:1089235/150‑1 (MQ=255)
                                                                                                                                    aaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCGCTCCCCTTGCTATTGcc                 <  2:260702/150‑1 (MQ=255)
                                                                                                                                     agcggCTATGCAGAAAATTGGGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCt                >  2:1329630/5‑150 (MQ=38)
                                                                                                                                       cggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGa              >  1:479438/3‑150 (MQ=38)
                                                                                                                                       cggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGa              >  1:1076811/3‑150 (MQ=38)
                                                                                                                                          cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTg           <  1:382401/150‑1 (MQ=255)
                                                                                                                                             tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTAc        >  1:362629/1‑150 (MQ=255)
                                                                                                                                              gCAGAAATTTGCACAGTGAGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTAcc       >  1:69189/1‑150 (MQ=255)
                                                                                                                                                   aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACaa  >  1:470203/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCATTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAA  >  NC_000913/2323126‑2323422

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