Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,451,538 C→T 6.5% intergenic (‑189/+59) tynA ← / ← paaZ tyramine oxidase, copper‑requiring/oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,451,5380CT 6.5% 356.0 / 9.5 107intergenic (‑189/+59)tynA/paaZtyramine oxidase, copper‑requiring/oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (50/50);  new base T (3/4);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGAATACCTCT  >  NC_000913/1451392‑1451684
                                                                                                                                                  |                                                                                                                                                  
tctGTGGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGt                                                                                                                                                  >  2:1061688/1‑149 (MQ=255)
  tgtgGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTca                                                                                                                                                >  3:1232790/1‑149 (MQ=255)
   gtgGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTcac                                                                                                                                               >  4:1123020/1‑149 (MQ=255)
      ggCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaa                                                                                                                                            >  4:587659/1‑149 (MQ=255)
       gCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacatta                                                                                                                                           <  1:1231095/149‑4 (MQ=255)
       gCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacatta                                                                                                                                           <  1:1367691/149‑4 (MQ=255)
        cGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAAc                                                                                                                                          <  3:587659/149‑1 (MQ=255)
           tttAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  3:1123020/149‑1 (MQ=255)
           tttAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  4:200589/149‑1 (MQ=255)
             tAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGt                                                                                                                                     <  2:1046664/149‑1 (MQ=255)
                tCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACt                                                                                                                                         <  3:620747/142‑1 (MQ=255)
                tCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACt                                                                                                                                         >  4:620747/1‑142 (MQ=255)
                      tCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaaaa                                                                                                                            <  1:1414574/149‑1 (MQ=255)
                      tCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaaaa                                                                                                                            <  4:1232790/149‑1 (MQ=255)
                          aTTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAgg                                                                                                                          >  4:355844/1‑147 (MQ=255)
                          aTTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAgg                                                                                                                          <  3:355844/147‑1 (MQ=255)
                           ttATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCa                                                                                                                       <  1:1061688/149‑1 (MQ=255)
                            tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaaa                                                                                                                                           >  2:1302266/1‑128 (MQ=255)
                            tATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaaaa                                                                                                                                           <  1:1302266/128‑1 (MQ=255)
                             aTCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       >  3:427907/1‑131 (MQ=255)
                             aTCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  4:427907/131‑1 (MQ=255)
                                 gTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacattatgc                                                                                                                                        >  4:1101129/1‑120 (MQ=255)
                                 gTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacattatgc                                                                                                                                        <  3:1101129/126‑7 (MQ=255)
                                 gTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaa                                                                                                                              >  4:1019117/1‑136 (MQ=255)
                                 gTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaa                                                                                                                              <  3:1019117/136‑1 (MQ=255)
                                     tGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                                      >  2:1314740/1‑124 (MQ=255)
                                     tGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAg                                                                                                                                      <  1:1314740/124‑1 (MQ=255)
                                            aaGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGt                                                                                                                                                  <  1:834781/105‑1 (MQ=255)
                                            aaGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGt                                                                                                                                                  >  2:834781/1‑105 (MQ=255)
                                                ggATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGcc                                                                                                             >  4:1020248/1‑138 (MQ=255)
                                                ggATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGcc                                                                                                             <  3:1020248/138‑1 (MQ=255)
                                                      ccTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCAcc                                                                                            >  1:103722/1‑149 (MQ=255)
                                                       cTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCGCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCtt                                                                                                    >  2:937566/1‑140 (MQ=255)
                                                       cTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCGCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCtt                                                                                                    <  1:937566/140‑1 (MQ=255)
                                                       cTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  3:1193628/105‑1 (MQ=255)
                                                       cTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       >  4:1193628/1‑105 (MQ=255)
                                                       cTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCg                                                                                           >  4:1170665/1‑149 (MQ=255)
                                                           aCCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       >  3:418613/1‑101 (MQ=255)
                                                           aCCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  4:418613/101‑1 (MQ=255)
                                                           aCCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGcc                                                                                                       >  2:351613/1‑133 (MQ=255)
                                                           aCCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGcc                                                                                                       <  1:351613/133‑1 (MQ=255)
                                                              tGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacattat                                                                                                                                          >  3:1454610/1‑91 (MQ=255)
                                                               gCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCg                                                                                   >  4:731694/1‑149 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacatta                                                                                                                                           >  2:23848/1‑87 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACGTTGTcacatta                                                                                                                                           <  1:23848/90‑4 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       <  2:269532/94‑1 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTa                                                                                                                                       >  1:269532/1‑94 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGccc                                                                                                      <  2:241893/127‑1 (MQ=255)
                                                                  aGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGccc                                                                                                      >  1:241893/1‑127 (MQ=255)
                                                                        tttGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaaaa                                                                                                                            <  1:1107606/99‑1 (MQ=255)
                                                                        tttGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATaaaa                                                                                                                            >  2:1107606/1‑99 (MQ=255)
                                                                              aGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGtgc                                                                    >  1:1267497/1‑149 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACaa                                                                                                                                             <  3:149/68‑1 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGc                                                                                                                        >  2:599311/1‑89 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGc                                                                                                                        <  1:599311/89‑1 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTc                                                                                                                <  2:1205447/97‑1 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTc                                                                                                                >  1:1205447/1‑97 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGtgc                                                                    <  2:352651/141‑1 (MQ=255)
                                                                                      cAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGtgc                                                                    >  1:352651/1‑141 (MQ=255)
                                                                                       aGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCAcc                                                                                            <  3:857128/116‑1 (MQ=255)
                                                                                         cTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCttca                                                                                                  >  4:503836/1‑108 (MQ=255)
                                                                                         cTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCttca                                                                                                  <  3:503836/108‑1 (MQ=255)
                                                                                          ttttCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGcc                                                           <  3:1055085/146‑1 (MQ=255)
                                                                                          ttttCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGcc                                                           >  4:1055085/1‑146 (MQ=255)
                                                                                               aTACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGt                                                   >  4:723564/1‑149 (MQ=255)
                                                                                                tACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTc                                                  >  2:331745/1‑149 (MQ=255)
                                                                                                 aCAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATAACCGTGCTGCCTGGCCACCAGCGTCa                                                 >  4:145421/1‑149 (MQ=255)
                                                                                                     tGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAAt                                             <  3:731694/149‑1 (MQ=255)
                                                                                                      ggATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAAtt                                            <  2:1267497/149‑1 (MQ=255)
                                                                                                         tGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATAACCGTGCTGCCTGGCCACCAGCGTCAGAATTGaa                                         <  3:145421/149‑1 (MQ=255)
                                                                                                            ccTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATc                                                                                    <  2:1446979/103‑1 (MQ=255)
                                                                                                            ccTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATc                                                                                    >  1:1446979/1‑103 (MQ=255)
                                                                                                             cTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGtgctgc                                                                 >  2:1002961/1‑121 (MQ=255)
                                                                                                             cTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGtgctgc                                                                 <  1:1002961/121‑1 (MQ=255)
                                                                                                             cTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGa                                          <  4:710981/144‑1 (MQ=255)
                                                                                                             cTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACa                                     >  4:1196351/1‑149 (MQ=255)
                                                                                                                   aTAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                          >  2:1322556/1‑106 (MQ=255)
                                                                                                                   aTAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCa                                                                          <  1:1322556/106‑1 (MQ=255)
                                                                                                                         gTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCgggg                         <  3:1196351/149‑1 (MQ=255)
                                                                                                                          tACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGt                        >  2:980254/1‑149 (MQ=255)
                                                                                                                            cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAgc                                   <  1:506695/136‑1 (MQ=255)
                                                                                                                            cacaTCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAgc                                   >  2:506695/1‑136 (MQ=255)
                                                                                                                               atcatGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGAt                   >  1:244715/1‑149 (MQ=255)
                                                                                                                               atcatGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGAt                   >  4:361517/1‑149 (MQ=255)
                                                                                                                               atcatGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGAt                   <  1:980254/149‑1 (MQ=255)
                                                                                                                                tcatCCATACTGTGACAACGTCACAAAACACAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCGTCATCCACTGAATAATCGACAAAATCACCGTGCTGCCTCGCCACGAGCGTCAGAATTGAATACAGCGCCGCCGGGGTTTGATg                  >  4:689770/1‑149 (MQ=255)
                                                                                                                                     cATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGAt             <  3:723564/149‑1 (MQ=255)
                                                                                                                                      aTAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGAtt            <  1:331745/149‑1 (MQ=255)
                                                                                                                                       tAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTg           >  3:517922/1‑149 (MQ=255)
                                                                                                                                         aTGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGaa         <  4:517922/149‑1 (MQ=255)
                                                                                                                                          tgtgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGccac                                                         <  1:955708/100‑1 (MQ=255)
                                                                                                                                          tgtgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGccac                                                         >  2:955708/1‑100 (MQ=255)
                                                                                                                                            tgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAAtt                                            >  1:134142/1‑111 (MQ=255)
                                                                                                                                            tgACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAAtt                                            <  2:134142/111‑1 (MQ=255)
                                                                                                                                             gACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTg                     <  2:201702/133‑1 (MQ=255)
                                                                                                                                               cAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATc                                                                           <  2:1462456/77‑1 (MQ=255)
                                                                                                                                               cAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATc                                                                           >  1:1462456/1‑77 (MQ=255)
                                                                                                                                               cAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGAATAcc     >  4:1282197/1‑147 (MQ=255)
                                                                                                                                               cAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGAATAcc     <  3:1282197/147‑1 (MQ=255)
                                                                                                                                                aaCGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGAATACctct  >  4:768166/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTGAAATAAAAGGGCAACTATTCGCCGTTGCCCTTCATTCACCGATTAATCGACAAAATCACCGTGCTGCCTGGCCACCAGCGTCAGAATTGAATACAGCGCCACCGGGGTTTGATGCTGATTGAATACCTCT  >  NC_000913/1451392‑1451684

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