Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,959,047 C→T 6.0% A74T (GCC→ACC)  cutC ← copper homeostasis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,959,0470CT 6.0% 207.4 / 2.8 66A74T (GCC→ACC) cutCcopper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCAT  >  NC_000913/1958899‑1959186
                                                                                                                                                    |                                                                                                                                           
gcgcGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGCCCGTGCGCACATCCTCAAGAATggc                                                                                                                                             <  6:411614/149‑1 (MQ=255)
gcgcGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggc                                                                                                                                             <  7:404088/149‑1 (MQ=255)
gcgcGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggc                                                                                                                                             <  3:267898/149‑1 (MQ=255)
 cgcgATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggcg                                                                                                                                            <  3:271121/149‑1 (MQ=255)
 cgcgATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggca                                                                                                                                            <  8:104407/149‑2 (MQ=255)
  gcgATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggcgg                                                                                                                                           <  5:202256/149‑1 (MQ=255)
   cgATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATggcggc                                                                                                                                          <  4:267413/149‑1 (MQ=255)
     aTGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCaa                                                                                                                                        <  7:225292/149‑1 (MQ=255)
          aaGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCAATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCa                                                                                                                                   >  6:375888/1‑149 (MQ=255)
                   ccAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGc                                                                                                                            <  4:99579/147‑1 (MQ=255)
                   ccAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTg                                                                                                                          <  8:24875/149‑1 (MQ=255)
                          aCCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCaa                                                                                                                    >  4:32258/1‑148 (MQ=255)
                            cGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCTACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCaaaaa                                                                                                                 <  3:4951/149‑1 (MQ=255)
                            cGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCaaaaa                                                                                                                 <  7:450960/149‑1 (MQ=255)
                                      cATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACcgc                                                                                                       >  6:134588/1‑149 (MQ=255)
                                      cATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACcgc                                                                                                       >  7:28456/1‑149 (MQ=255)
                                      cATTATTTTTTCCATTCGTAGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACcgc                                                                                                       >  8:419160/1‑149 (MQ=255)
                                       attattTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGTGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACcgcg                                                                                                      >  7:232219/1‑149 (MQ=255)
                                        ttattTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGt                                                                                                     >  3:382789/1‑149 (MQ=255)
                                         tattTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTg                                                                                                    >  2:66237/1‑149 (MQ=255)
                                          attTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTg                                                                                                    <  8:230352/148‑1 (MQ=255)
                                             ttttCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGc                                                                                                 <  5:134588/148‑1 (MQ=255)
                                             ttttCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGc                                                                                                 <  3:32258/148‑1 (MQ=255)
                                                   aTTCGTGGCATATCGACATTCCCGTCATCATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATcac                                                                                                             <  5:368572/130‑1 (MQ=255)
                                                   aTTCGTGGCATATCGACATTCCCGTCAGCATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATcac                                                                                                             >  6:368572/1‑130 (MQ=255)
                                                   aTTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTAt                                                                                          <  1:416702/149‑1 (MQ=255)
                                                   aTTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTAt                                                                                          >  5:168074/1‑149 (MQ=255)
                                                   aTTCGTGGCATATCCACATTCCCGTCAACATCGAGAACTCCCGTCACCAGTCCAGGAAAACCTAATTCGCTTACCGTGCGAACATCCTCAAGAATGTCGGCAAACTCACAGTCGCTGTAGAAAAAATcacc                                                                                                            >  7:151061/1‑131 (MQ=255)
                                                   aTTCGTGGAATATCGACATTCCCGTAAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGTACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCAGCGTCGCTGTAGGATAAATcacc                                                                                                            <  8:151061/131‑1 (MQ=255)
                                                    ttCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATc                                                                                         <  2:418375/149‑1 (MQ=255)
                                                       gTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGa                                                                                      <  4:240599/149‑1 (MQ=255)
                                                        tgtcATATAGACATTACCGTCATCATTGAGAAGTTCCGTGACCAGTCCATGATTACCTAATTCTCGGACCGTGCGCACATACTCAAGAATGGTGTCAAACTCTCCTTCTCTGTAGCAAAAATCACCACCGCTTTGGCGAATTATCGGAt                                                                                     <  8:232219/146‑1 (MQ=255)
                                                        tGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGAt                                                                                     <  3:73274/149‑1 (MQ=255)
                                                            atatCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGca                                                                                  >  8:293647/1‑148 (MQ=255)
                                                               tCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCCAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATcac                                                                                                             >  6:177905/1‑118 (MQ=255)
                                                                cGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAgg                                                                              >  1:339999/1‑148 (MQ=255)
                                                                 gACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                                             <  8:28456/148‑1 (MQ=255)
                                                                  aCATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTTGGGGAATTATCGGCTGCACAgggct                                                                           >  6:309344/1‑147 (MQ=255)
                                                                  aCATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGAt                                                                           <  7:141975/149‑1 (MQ=255)
                                                                   cATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATc                                                                          <  7:419160/149‑1 (MQ=255)
                                                                          gTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCCGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCAccc                                                                   <  5:375888/149‑1 (MQ=255)
                                                                           tCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                                                  >  4:409019/1‑149 (MQ=255)
                                                                              aCATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTg                                                               <  2:339999/149‑1 (MQ=255)
                                                                                         gCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGTCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGAtt                                                    >  8:42327/1‑149 (MQ=255)
                                                                                           ccGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGAttt                                                   >  1:175652/1‑148 (MQ=255)
                                                                                           ccGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGAttt                                                   >  2:321501/1‑148 (MQ=255)
                                                                                            cGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCa                                                 <  7:293647/149‑1 (MQ=255)
                                                                                            cGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGAGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCa                                                 <  6:80416/149‑1 (MQ=255)
                                                                                             gTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                                                <  6:168074/149‑1 (MQ=255)
                                                                                                aCCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGTCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTa                                              <  7:42327/148‑1 (MQ=255)
                                                                                                aCCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTa                                              <  4:382789/148‑1 (MQ=255)
                                                                                                aCCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGAACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGTCGAATTATCGGATGCACAGGGATCGTCACCCGCTCGCGCACGGATTTCAGTa                                              >  2:442008/1‑148 (MQ=255)
                                                                                                     tCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGCCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCa                                        >  6:273599/1‑149 (MQ=255)
                                                                                                       cAGGAAAACCTACTTCGCGGACCGTGCGCACATCCTCCCGAATGGCGGCACACTCACCGTCGCGGTAGCAAACATCCCCCCCGCGTGGGCGACTTATCGGCTCCACAGGGAGCGTCACCCGCTTGCGCACGGATTTCAGTACACCCAac                                      <  3:375567/149‑1 (MQ=255)
                                                                                                        aGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGTGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAacg                                     >  2:275213/1‑149 (MQ=255)
                                                                                                        aGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTa                                              <  7:135791/140‑1 (MQ=255)
                                                                                                        aGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTa                                              >  8:135791/1‑140 (MQ=255)
                                                                                                            aaaCCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACgg                                 >  7:381649/1‑149 (MQ=255)
                                                                                                                cTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGt                              >  8:90203/1‑148 (MQ=255)
                                                                                                                cTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCGGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCACCGACGGCGt                              >  2:97734/1‑148 (MQ=255)
                                                                                                                   aTTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTaa                           <  3:409019/148‑1 (MQ=255)
                                                                                                                    ttCGCGGACCGTGCGCACATCCTCAAGAATGGTGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGc                         <  1:275213/149‑1 (MQ=255)
                                                                                                                       gcgGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGcccc                      >  7:360924/1‑149 (MQ=255)
                                                                                                                               tGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCtttt              <  7:90203/149‑1 (MQ=255)
                                                                                                                                    acaTCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGc         >  4:220059/1‑149 (MQ=255)
                                                                                                                                    acaTACTCAATAATGTCCCCAAACCCACCGTCCCTGTAGCAAAAATCACAACCGCGTGGTCGAATTATCGGATGCACAGCGATCGTCACCCGATGGCGAACGGAATTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCAATTGGGGc         >  8:180315/1‑149 (MQ=255)
                                                                                                                                      aTCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTg       >  6:200022/1‑149 (MQ=255)
                                                                                                                                       tCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCGGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTgc      <  1:66237/149‑1 (MQ=255)
                                                                                                                                          tCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGCTTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCa   >  5:414858/1‑149 (MQ=255)
                                                                                                                                          tCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCa   >  5:249517/1‑149 (MQ=255)
                                                                                                                                          tCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCa   >  4:329988/1‑149 (MQ=255)
                                                                                                                                            aaGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCAt  <  3:323437/148‑1 (MQ=255)
                                                                                                                                               aaTGCCGGCAAATTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGGTAATAATGGGATGCACAGGGATCGGCACCGGCTGGCGCACGGATTTCATTACACCAAACGACGGCGTTAAGCCCCCCTCTTTTggg           <  5:54498/136‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCCCGTCACCAGTCCAGGAAAACCTAATTCGCGGACCGTGCGCACATCCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAAGCCCCCCTCTTTTGGGGCTGCGCAT  >  NC_000913/1958899‑1959186

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