Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,957 T→A 7.0% D66E (GAT→GAA dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9570TA 7.0% 389.1 / 12.7 140D66E (GAT→GAAdinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base T (70/60);  new base A (2/8);  total (72/68)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 150 total aligned reads displayed.

AACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCA  >  NC_000913/3817808‑3818106
                                                                                                                                                     |                                                                                                                                                     
aaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                                       <  1:1584926/150‑1 (MQ=255)
aaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                                       >  1:1329915/1‑150 (MQ=255)
 aCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACATTTCGAGGAGATCCTCGatat                                                                                                                                                     <  2:196736/150‑1 (MQ=255)
  cGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGGCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                                                                                     <  1:421/150‑1 (MQ=255)
   ggAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCGCGATATg                                                                                                                                                    >  1:1257910/1‑150 (MQ=255)
    gAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATgg                                                                                                                                                   <  1:421826/150‑1 (MQ=255)
    gAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACATTGGGCACACAATTGTTGACCATATCGAGGATCTCCTCGAAATGGt                                                                                                                                                  <  2:49954/150‑1 (MQ=255)
       cAGAACAATGGGCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTTTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGTATCTCCTCGAAATGGTCa                                                                                                                                                <  1:1371661/150‑1 (MQ=255)
         caaCAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCatc                                                                                                                                              <  1:775407/149‑3 (MQ=255)
           aCAATGGTCAGCAAGACAACTTGGGAAACTACTGTGTTATTCAGAGTATCGTCACTATAAACCTGTATTACCGCGCGCCAAAGCAGCCTGTGAGAACACTGGTCACACAATTGATGACCATTTCTAGGAGATCCTCGATATGGTCAAAAt                                                                                                                                            >  1:266493/1‑150 (MQ=255)
             aaTGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                          <  1:667318/150‑1 (MQ=255)
              attgTCAGAAAGACAACTTTGGAAACTACTGGGTTATTCATAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                                         <  1:651250/147‑1 (MQ=255)
                ggTCAGCAAGACAAGTTGGGAAACTACGGGGTTATTCAGAGTATCGTCACTTTATACCTGTCTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                       <  2:1041660/150‑1 (MQ=255)
                ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                       >  1:1122669/1‑150 (MQ=255)
                ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                       >  2:172233/1‑150 (MQ=255)
                ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCCTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                       <  1:722699/150‑1 (MQ=255)
                 gTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTCTACCTGTATTAACGCGCGCCAAAGAAGCCGGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTc                                                                                                                                      <  1:1090050/150‑1 (MQ=255)
                  tCAGCAAGACAACTTGGTAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTCTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGAGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCa                                                                                                                                     >  1:288748/1‑150 (MQ=255)
                  tCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCa                                                                                                                                     <  1:99098/150‑1 (MQ=255)
                   cAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaa                                                                                                                                    >  1:1636690/1‑150 (MQ=255)
                    aGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGCAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaaa                                                                                                                                   >  2:100827/1‑150 (MQ=255)
                     gCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                                                  <  1:770233/150‑1 (MQ=255)
                     gCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                                                  >  2:726804/1‑150 (MQ=255)
                      cAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATAAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATg                                                                                                                                 <  2:1186000/150‑1 (MQ=255)
                             aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAag                                                                                                                          >  1:1484273/1‑150 (MQ=255)
                              cTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                                         <  1:314676/150‑1 (MQ=255)
                                tGGGAACCTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGc                                                                                                                       <  2:1492822/150‑1 (MQ=255)
                                tGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGc                                                                                                                       >  2:613450/1‑150 (MQ=255)
                                 gggAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCGAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                                      >  2:1702876/1‑150 (MQ=255)
                                 gggAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                                      <  2:635317/150‑1 (MQ=255)
                                       ctactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                                >  2:155817/1‑150 (MQ=255)
                                        tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAg                                                                                                               >  2:1556840/1‑150 (MQ=255)
                                        tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAg                                                                                                               <  1:739870/150‑1 (MQ=255)
                                         actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                                    <  1:408633/144‑1 (MQ=255)
                                         actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                                    >  2:408633/1‑144 (MQ=255)
                                         actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAAAGAGCGCCAAAGCAGCCTGTGAGAACAGTGGACACACAATTAATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATGGGCTCAAATGCCAAAAGAACATTAAAAgc                                                                                                              >  2:283866/1‑149 (MQ=255)
                                            gggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAGATGCCAAAAGAGCATTAAAAGACAt                                                                                                           <  1:1078416/150‑1 (MQ=255)
                                               ttATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGt                                                                                                        <  2:1442097/150‑1 (MQ=255)
                                                  ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAct                                                                                                     >  1:1455583/1‑150 (MQ=255)
                                                    cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAatct                                                                                                   >  1:381453/1‑150 (MQ=255)
                                                      gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTcc                                                                                                 >  1:1007769/1‑150 (MQ=255)
                                                       agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCAAAAATAGCATTAAAAGACATCGTACTCTccc                                                                                                >  2:871366/1‑150 (MQ=255)
                                                         tATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGc                                                                                              <  1:1159448/150‑1 (MQ=255)
                                                            cGTCAATTTATACCTGTATTAACGCGCGCCAAAGAAGCATGTGAGAACAGTGGACACACAATTGATGACCATTTCGAGTAGATCATCGATATCGAAAAAATTGGCACAAATCCCAAAAGAGAATAAAAACACAAATTACTCTccacctaa                                                                                           >  2:462349/1‑144 (MQ=255)
                                                               cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGTGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                        <  2:974108/150‑1 (MQ=255)
                                                               cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                        <  2:1586847/150‑1 (MQ=255)
                                                                       aCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCACCCGCTATGCCTGTTACCt                                                                                >  1:1257483/1‑150 (MQ=255)
                                                                         cTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                              >  2:832429/1‑150 (MQ=255)
                                                                         cTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCAATGCCTGTTACCTgg                                                                              <  2:399628/150‑1 (MQ=255)
                                                                             aTTAACGCGCGACAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagt                                                                          <  2:1198719/150‑1 (MQ=255)
                                                                              ttAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta                                                                         <  2:36487/150‑1 (MQ=255)
                                                                               tAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTAc                                                                        <  1:1058630/150‑1 (MQ=255)
                                                                                aaCGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACa                                                                       <  1:413791/150‑1 (MQ=255)
                                                                                aaCGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACa                                                                       <  2:1036773/150‑1 (MQ=255)
                                                                                  cgcgcgCCAAAGAAGCCTGCGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                                     >  1:472100/1‑150 (MQ=255)
                                                                                    cgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAAc                                                                   <  1:1033026/150‑1 (MQ=255)
                                                                                    cgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAAc                                                                   <  2:1288473/150‑1 (MQ=255)
                                                                                      cgcCAAAGAAGCCTGAGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGCTCCGCGATATGGTCAAAATTGGCTCAACTGCCAAAAGAGCATTAAAAGACATCGTACTCACCCGCTATGCCTGTTACCTGGTAGTACAAAACgg                                                                 >  2:265924/1‑150 (MQ=255)
                                                                                       gcaaaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATACTCGATATGGTCAAAATAGGCTCAAATGCCAAAAGAGCAATAAAAGACATCGTACTCTCCCCATAAGCCTGATACCAGGAAGTACAAAAcccc                                                                >  2:63507/4‑147 (MQ=255)
                                                                                       gcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                                                >  1:1003762/1‑150 (MQ=255)
                                                                                       gcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                                                <  1:614188/150‑1 (MQ=255)
                                                                                         cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCgc                                                              >  2:1138278/1‑149 (MQ=255)
                                                                                          aaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAc                                                             >  1:758937/1‑150 (MQ=255)
                                                                                           aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAcc                                                            <  2:242752/150‑1 (MQ=255)
                                                                                                  cTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCg                                                                                               >  2:317817/1‑108 (MQ=255)
                                                                                                   tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGTAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                                    >  2:1595426/1‑150 (MQ=255)
                                                                                                   tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                                    >  2:1690471/1‑150 (MQ=255)
                                                                                                   tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                                    >  2:1005007/1‑150 (MQ=255)
                                                                                                   tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                                    >  2:1011448/1‑150 (MQ=255)
                                                                                                      gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                                                 >  2:1627336/1‑150 (MQ=255)
                                                                                                        gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGt                                               <  1:1406267/150‑1 (MQ=255)
                                                                                                          aCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCa                                             >  2:505885/1‑150 (MQ=255)
                                                                                                           cAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCAt                                            >  2:738044/1‑150 (MQ=255)
                                                                                                            aGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCAtt                                           <  2:1492155/150‑1 (MQ=255)
                                                                                                              tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGc                                         <  2:1329915/150‑1 (MQ=255)
                                                                                                                 tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggc                                      >  1:874600/1‑150 (MQ=255)
                                                                                                                   acacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcag                                    <  1:721343/150‑1 (MQ=255)
                                                                                                                    cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagg                                   >  2:664114/1‑150 (MQ=255)
                                                                                                                      caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggc                                 >  2:64425/1‑150 (MQ=255)
                                                                                                                       aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggca                                <  1:155817/150‑1 (MQ=255)
                                                                                                                             tgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTa                          >  1:551521/1‑150 (MQ=255)
                                                                                                                              gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAt                         >  1:528003/1‑150 (MQ=255)
                                                                                                                              gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAt                         >  2:5662/1‑150 (MQ=255)
                                                                                                                               aCCATATCGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt                        >  2:942935/1‑150 (MQ=255)
                                                                                                                                ccATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAttt                       >  1:891896/1‑150 (MQ=255)
                                                                                                                                 cATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtttt                      >  1:1316906/1‑150 (MQ=255)
                                                                                                                                  aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg                     >  1:559291/1‑150 (MQ=255)
                                                                                                                                    atCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCt                   <  1:1706968/149‑1 (MQ=255)
                                                                                                                                      cGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTAt                 >  1:174009/1‑150 (MQ=255)
                                                                                                                                      cGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGACAAAAGAGCATAAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCTGCAGGGCAGACTTATTTTGCTAt                 <  1:838527/150‑1 (MQ=255)
                                                                                                                                       gaggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAACCCGGTCATTGCGGCAGGGAAGACTTATTTTGCTATc                >  1:834078/1‑150 (MQ=255)
                                                                                                                                        aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATcc               <  2:1636690/150‑1 (MQ=255)
                                                                                                                                        aggagATCCTCGATATGGTCAAAATTGGCTCAAATGACAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATcc               >  2:1127236/1‑150 (MQ=255)
                                                                                                                                         ggagATCCTCGATATTGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCACGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCa              <  2:288748/150‑1 (MQ=255)
                                                                                                                                             aTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAATAGCATTAATAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAcc          <  1:63507/150‑1 (MQ=255)
                                                                                                                                             aTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAcc          <  2:528003/150‑1 (MQ=255)
                                                                                                                                             aTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAcc          >  2:530139/1‑150 (MQ=255)
                                                                                                                                                cTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCAATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGa       >  2:1439241/1‑150 (MQ=255)
                                                                                                                                                    atatGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGc   >  2:333348/1‑150 (MQ=255)
                                                                                                                                                     tatGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCa  <  1:355892/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCA  >  NC_000913/3817808‑3818106

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