Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 509,371 A→C 14.8% V670G (GTG→GGG)  copA ← copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913509,3710AC14.8% 183.5 / 3.1 72V670G (GTG→GGG) copAcopper transporter
Reads supporting (aligned to +/- strand):  ref base A (35/24);  new base C (4/9);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 7.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGTTTCAATGGCAACATCACTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCG  >  NC_000913/509136‑509604
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cgcgGTTGCGATGGCAACATCCCTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCa                                                                                                                                                                                                                                        <  1:43053/239‑1 (MQ=255)
         aaTGGCAACCTCCGCGCCGCCAACCATCGCAATGCCGACAGCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGGCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGatacc                                                                                                                                                                                                                               <  2:1080570/239‑3 (MQ=255)
                 cATCACTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCCGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGCGGGTGGGTTATCCCCGGTGAACATTCCCAGACGTAATCCCGCcttt                                                                                                                                                                                                                       >  2:707440/1‑238 (MQ=255)
                    cccTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGc                                                                                                                                                                                                                    <  2:522181/237‑1 (MQ=255)
                    cACTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGc                                                                                                                                                                                                                    <  1:716142/239‑1 (MQ=255)
                     aCTGCCGCCCCCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCa                                                                                                                                                                                                                   <  1:394412/239‑1 (MQ=255)
                            ccACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGGTGGGTTATCCCCGGTCAACATCCCCCGACGAAATACCGCCTCTTTCAGGCGttt                                                                                                                                                                                                            >  2:563696/1‑238 (MQ=255)
                               cccAGCGCACGGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCCCTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTCCCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAg                                                                                                                                                                                                         <  2:766917/239‑1 (MQ=255)
                                ccATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAgc                                                                                                                                                                                                        <  2:420472/239‑1 (MQ=255)
                                     gcatggCCGGCATCCGCCGGAGCCAGCGCTGGGTGATCGTTAGTCCCCTCGCACACCATTGTCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTCCCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGcgccgc                                                                                                                                                                                                   <  2:728833/234‑1 (MQ=255)
                                     gCACTGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGcgccgc                                                                                                                                                                                                   <  2:61986/239‑1 (MQ=255)
                                     gCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGcgccgc                                                                                                                                                                                                   <  1:565643/239‑1 (MQ=255)
                                         tGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGGTATCCCCGGTCAACATCACCCGACGATATCCCGCTTTATGCAGGCGTTGCCGCGCCGCCAc                                                                                                                                                                                                >  1:296507/1‑238 (MQ=255)
                                            cGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTa                                                                                                                                                                                            <  1:146169/239‑1 (MQ=255)
                                                   cGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCa                                                                                                                                                                                         <  1:298191/235‑1 (MQ=255)
                                                         aGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCg                                                                                                                                                                                                                            <  2:1131994/194‑1 (MQ=255)
                                                         aGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCg                                                                                                                                                                                                                            >  1:1131994/1‑194 (MQ=255)
                                                               cgcTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAcc                                                                                                                                                                                                                                         <  1:1133802/175‑1 (MQ=255)
                                                               cgcTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAcc                                                                                                                                                                                                                                         >  2:1133802/1‑175 (MQ=255)
                                                                gcTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATcgc                                                                                                                                                                         <  1:563696/238‑1 (MQ=255)
                                                                          cgtTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCCCCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCcagc                                                                                                                                                              <  1:297170/239‑1 (MQ=255)
                                                                                     tCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCCGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTAACACCACGCAACGGAACGCGTACTTGCAGCCCGGGTAcc                                                                                                                                                     >  2:417679/1‑237 (MQ=255)
                                                                                       gCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCt                                                                                                                                                  <  1:1129463/238‑1 (MQ=255)
                                                                                           accaTTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAcc                                                                                                                                                                                                                                         <  1:407070/147‑1 (MQ=255)
                                                                                             caTTGCCACCTGACGTCCTTCCCTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTa                                                                                                                                                                      <  1:26272/212‑1 (MQ=255)
                                                                                                tGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAg                                                                                                                                                                                                         <  2:564481/174‑1 (MQ=255)
                                                                                                tGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAg                                                                                                                                                                                                         >  1:564481/1‑174 (MQ=255)
                                                                                                     ccTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGCGGTTGGGTTTTCCCCTGTCAACATCCCCAGACGATATc                                                                                                                                                                                                                               >  2:407070/1‑147 (MQ=255)
                                                                                                       tGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGGGGTGCAGCGCCGCCACGCCATCCCCACTCCACGGGTCGCCTACTGCCAGCCAGGGTACCCCCTTTCCGTTAACCCCca                                                                                                                                  >  2:782557/1‑238 (MQ=255)
                                                                                                       tGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGGTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACCCTATCACTACGCAACGGAATGCGGACTGCCACCAGGGCGGCCCCTTTtccacaccagcgct                                                                                                                                  >  2:26272/1‑227 (MQ=255)
                                                                                                         aCGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAcc                                                                                                                                                                                                                                         >  1:727712/1‑133 (MQ=255)
                                                                                                         aCGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAcc                                                                                                                                                                                                                                         <  2:727712/133‑1 (MQ=255)
                                                                                                         aCGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCCGACGATATCCCGCTTTATTCAGGCGTTGCAGCGCCGCCACGCCATCACTACGCATCGGATCGCGGACTGCCAGCAGGGGTACCGCTTTTCCGTCCAccgcccgc                                                                                                                                >  2:829031/1‑238 (MQ=255)
                                                                                                                 cACTTTGCAGATGTGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCCGTCAACATCACCAGACGATATCCCGCTTTATTCAGGCGTTGCAGCGCCGCCACGCTATCACTACCCAACCGAACGCGGACTGCCCACAGGGCCACCGCCTTTCCGTCAACCGCCCGCAGccca                                                                                                                          >  2:1153829/1‑238 (MQ=255)
                                                                                                                        cAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCCGCACCGGCGTTGcc                                                                                                                 >  1:369407/1‑238 (MQ=255)
                                                                                                                        cAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCCCGGACTGCCAGCAGGGCTACCGCCTTTCCGTCAACCGCCAGCAGCACAGGCGTTGcc                                                                                                                 >  1:1096216/1‑238 (MQ=255)
                                                                                                                              tGTTAGACCGCTTCGCTTTTCCCGTCCGGCAACACCCGGGCGATCCCCTCATCAACCCCTGTTTCTTGGGCTTCCGCATTGCCGGTGTTGGGGTTTCCCCCGGCCACCACCCCCAGCCGATACCCCGCTTTATGCGGGCGTGCCAGCCCCGCCACGCTATCCCTCCGCAACGGACCGCGTACTGCCAGCAGGGCTCCCGCTTCCCCGTAAACCGCCAGCAGCACAGGCGTTGccc                                                                                                                <  2:584265/235‑1 (MQ=255)
                                                                                                                                tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCCGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCCTTTCCGGCAACCGCCCGCAGAACAgggcttggg                                                                                                                 >  2:560524/1‑225 (MQ=255)
                                                                                                                                tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCCGCAGCACAGGCGTTGCCCCTTGCg                                                                                                          >  1:654411/1‑239 (MQ=255)
                                                                                                                                    gcctcTTCGGCTTTACCGTCCGGCACCCCCCCGGCGATCACCTCACCAACCCCTCCTCTTTTGGCGCCCGCATGGGCGGTGGTTGGTTTCCCCCCGGCAAACCCCCCCAGGCGATACCCCGCTTTATGCGGGCGTTGCAGCGCCGCCACGCTATCACTCCCCACCGGATCGCGTCCTGCCAGCAGGGCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGc                                                                                                      <  2:221381/235‑1 (MQ=255)
                                                                                                                                         ttCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCCACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGa                                                                                                  >  1:779727/1‑238 (MQ=255)
                                                                                                                                          gggttttccccGCCCGGCAGCCCCCGGGCAACCCCCTCACCAACCCCTGCTTCTTTGCCGTCCGCTTGGGCGGTGGTGGGTTTACCCCGGGCAACCACCCCCAGACGATCCCCCGCTTTATGCAGGCGTTGCAGCGCCGCCCCGCTACCCCTCCGCACCGGATCGCGTACGGCCAGCAGGGCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGccc                                                                                                                <  2:536591/214‑1 (MQ=255)
                                                                                                                                                      tCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTcc                                                                                      >  1:445921/1‑237 (MQ=255)
                                                                                                                                                      tCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGGCAACCGCCAGCAGCCCAGGCGTTGCCCCTTGCGATGc                                                                                                      >  1:536591/1‑221 (MQ=255)
                                                                                                                                                       ccGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCt                                                                                   >  1:813193/1‑239 (MQ=255)
                                                                                                                                                                ccccGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGcgccgc                                                                                                                                                                                                   <  1:357180/118‑1 (MQ=255)
                                                                                                                                                                ccccGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGcgccgc                                                                                                                                                                                                   >  2:357180/1‑118 (MQ=255)
                                                                                                                                                                ccccGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGGAATCTCCGCTTCGATAGCt                                                                          >  1:205052/1‑239 (MQ=255)
                                                                                                                                                                         tccccTCATCATCCCTTGCTTTTTGGGCATCCGATTTGGCGGTGTTTGGTTATCCCCGGTCCACCACCCCCAGCCGATCCCCCGTTTTAGGCGGGCTTTGCAGCGCCGCCACGTTTCCCCTCCGCAACGGATCGCGTCCTGCCAGCAGGCCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCa                                                                 <  2:205052/236‑1 (MQ=255)
                                                                                                                                                                           aCCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCCCGCTATCACTACGCAACGGATCCCGTACTGCCAGCCGGGCTACCGCTTTCCCGTCCACCGCCAGCAGCACAGGCGTTGCCCCTTTCCATGCCTGGGCCGTAATTCCCGCTTCCATAACTTTTGTACCaa                                                                >  2:423249/1‑238 (MQ=255)
                                                                                                                                                                           aCCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCaa                                                                >  1:245767/1‑238 (MQ=255)
                                                                                                                                                                           aCCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCaa                                                                >  1:478381/1‑238 (MQ=255)
                                                                                                                                                                           aCCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCaa                                                                >  1:769001/1‑238 (MQ=255)
                                                                                                                                                                             cTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCCCGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTTCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGGAccca                                                                >  1:430698/1‑234 (MQ=255)
                                                                                                                                                                               catcaATCCTTGTTTCTTGGGCGTTCGCTTTGCCGGTGTTTGGTTTCCCCCGGGCAAACACCCCCAGCCGATACCCCGCTTTATGCGGGCGTTGCAGCGCCGCCACGCTATCACTCCGCACCGGATCGCGTACTGCCAGCAGGGCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCtt                                                                         <  2:150395/225‑1 (MQ=255)
                                                                                                                                                                                          tcttgttGGGCGTTCTCTTGGGGGGTGTTGGGTTTACCCCCGGTAAACACCCCAGGACGATACCCGGCTTAATGCGGGCGTTGCGGCGCCGCCACGCTATCCCTCCGCACCGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTaa                                                <  2:245767/238‑1 (MQ=255)
                                                                                                                                                                                          gCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCaaa                                                               >  1:150395/1‑223 (MQ=255)
                                                                                                                                                                                          gCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCCACCGCTTTCCCGTCAACCGCCCGCAGCACAGGCGTTGCCCCTTTCGATGCCTGAGCAGTAATTCTCGCTTCGCTAGCTTTGGGACCAACCTGTGGCTCAGTTa                                                 >  2:33425/1‑238 (MQ=255)
                                                                                                                                                                                              ttttGGCTATGGCATTGGCGGTGGTTGGTTTACCCCGGGTCAACCTCCCCAGACGATCCCCCGTTTTATGCGGGCGTTGCAGCCCCGCCACGCTTCCACCCCGCAACGGATCGCGTCCTGCCAGCGGGGCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTa                                                                    <  2:106678/214‑1 (MQ=255)
                                                                                                                                                                                                      aTCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAgcg                                                                                                                                                                                                       <  2:825072/76‑1 (MQ=255)
                                                                                                                                                                                                      aTCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAgcg                                                                                                                                                                                                       >  1:825072/1‑76 (MQ=255)
                                                                                                                                                                                                        cGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTTCCCCTTGCGATGCCTGAGCACTAATCTCCGCGTCGATAGCCTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTg                                   >  2:1231733/1‑238 (MQ=255)
                                                                                                                                                                                                        cGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGtt                                                          >  1:106678/1‑215 (MQ=255)
                                                                                                                                                                                                                  gtggtTGGGGTATCCCCGGTCAACATCACCCGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCCCCAGCAGCACAGGCGTTGCCCCTTTCGATGCCTGAGCAGTAATCTCCGCTTCCATAGCTTTGGTACCACCCCGTTGGTCATTTAAAAGCGCCTGATTGCCCCGCCATa                         >  2:630474/1‑238 (MQ=255)
                                                                                                                                                                                                                     ttaggtttACCCCGGTTCAACACCACCAGACGATACCCCGTTTTACGCCGGCGTTGCAGCGCCGCCACGTTACCCCTACCCAACGGATCGCGTACTGCCAGCGGGGCTCCCGTTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGcc                                 <  2:484750/221‑1 (MQ=255)
                                                                                                                                                                                                                     gtTGGGGTATCCCCGGTCAACATCCCCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCCCGCTATCCCTACCCAACGGATCGCGTCCTGCCAGCCCGGCTACCCCTTTTCCGTCAACCGCCAGCAGCCCAGGCGTTTCCCCTTGCGCTGCCTGGGCGGTACTCTCCACCTCGCTCACTTTGGGACCCACCCGTTTGTCACTTAACACCGCCCGCTTGcccccagacaaagc                     >  2:194714/1‑228 (MQ=255)
                                                                                                                                                                                                                        gggTTACCCCGGGTCAACATCACCAGACGATCTCCCGTTTTATGCGGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTAAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCAt                   <  1:423249/238‑1 (MQ=255)
                                                                                                                                                                                                                          tttACCCCCGGACAACATCCCCAGACGATACCCCGTTTTAGGCAGGCTTGGCAGCCCCCCCACGCTTTCCCTCCGCAACGGTTCGCGTACTGCCAGCGGGGCTCCCGCTTCCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGAt                                                                              <  2:1183511/176‑1 (MQ=255)
                                                                                                                                                                                                                               ccccGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAAc                       >  1:484750/1‑227 (MQ=255)
                                                                                                                                                                                                                                    gTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTa                                                                    >  1:1183511/1‑177 (MQ=255)
                                                                                                                                                                                                                                    gTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTc      >  1:804306/1‑239 (MQ=255)
                                                                                                                                                                                                                                           cATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCg  >  1:1073002/1‑238 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
CGCGGTTTCAATGGCAACATCACTGCCGCCACCCATCGCAATGCCGACATCCGCCTGAGCCAGCGCTGGCGCGTCGTTAATGCCGTCGCCCACCATTGCCACCTGACGTCCTTCACTTTGCAGA‑‑TGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCG  >  NC_000913/509136‑509604

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