Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 50,588 G→A 6.3% A212V (GCG→GTG)  apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091350,5880GA 6.3% 197.9 / 3.0 64A212V (GCG→GTG) apaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/2);  total (31/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.48e-01

TTTATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAGCCCCCGCAATTCCGGTGA  >  NC_000913/50440‑50733
                                                                                                                                                    |                                                                                                                                                 
tttATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCg                                                                                                                                                   <  5:145260/149‑1 (MQ=255)
  tATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGc                                                                                                                                                  >  1:213098/1‑148 (MQ=255)
    tCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCaaa                                                                                                                                               >  5:170341/1‑149 (MQ=255)
       tCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCa                                                                                                                                            >  7:97980/1‑149 (MQ=255)
        cccAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAAc                                                                                                                                              <  1:218194/146‑1 (MQ=255)
               caggcaggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCa                                                                                                                                    <  7:268606/149‑1 (MQ=255)
                    aggTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTggg                                                                                                                               <  1:67680/149‑1 (MQ=255)
                      gTTTATGTCCCACCCCAGCCGCAGCCGGTATCCCGCGCGTATATACCGTCCGGCGTCCCTTTGCCCTCCCGCGAGGCCCCGTGACCAAAGGCGATGATGTATTCTTCAGCGACAGGGCCAGGATTCGCAAACCCCGGTTTCAGTgggg                                                                                                                              >  3:30408/1‑148 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTgggg                                                                                                                              >  8:75139/1‑148 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTgggg                                                                                                                              >  7:211140/1‑148 (MQ=255)
                      gTTAATGTACCACCCCAGCAGCAGCCGGTATACAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTgggg                                                                                                                              >  7:132196/1‑148 (MQ=255)
                          aTGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAggg                                                                                                                         >  2:135577/1‑149 (MQ=255)
                          aTGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAggg                                                                                                                         <  7:3355/149‑1 (MQ=255)
                            gTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc                                                                                                                       >  8:288559/1‑149 (MQ=255)
                             tACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc                                                                                                                       <  6:110149/148‑1 (MQ=255)
                                    ccAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCg                                                                                                               >  6:26185/1‑149 (MQ=255)
                                    ccAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCg                                                                                                               <  5:105957/149‑1 (MQ=255)
                                               cGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGGGTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTg                                                                                                    <  7:247408/149‑1 (MQ=255)
                                                 gTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGAGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCACAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCt                                                                                                  >  8:73894/1‑149 (MQ=255)
                                                      ccgggcgTTTATTGCTTCCGGCGTACCTTTGCCCTGCTGCGATGGCCATTTAGCATTTGAGTTTCTGTATTCTTCAGCGTCTTGGCCAGGAATCGCAAACCACGGTTTCTGTGGGGCCGGTGTCTCTTCCGGCGATTCTTTGCTGTACa                                                                                             <  2:17653/145‑1 (MQ=255)
                                                       aGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACa                                                                                             >  7:81508/1‑148 (MQ=255)
                                                         cgcgTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACatat                                                                                          >  4:311363/1‑149 (MQ=255)
                                                            gTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATc                                                                                         >  1:249889/1‑147 (MQ=255)
                                                            gTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATc                                                                                         >  8:328547/1‑147 (MQ=255)
                                                                   ccTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGAcc                                                                                >  5:137824/1‑149 (MQ=255)
                                                                       ccGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGtt                                                                             <  7:288559/148‑1 (MQ=255)
                                                                         ggCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCgg                                                                          <  6:170341/149‑1 (MQ=255)
                                                                            gTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCCGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGaa                                                                       <  5:189804/149‑1 (MQ=255)
                                                                            gTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCggg                                                                         >  6:25250/1‑147 (MQ=255)
                                                                            gTACCTTTGCCCTCCAGCGATGCCCAGAGACCAAACGCGATGCTGTATTCTTCACCGACAGGGCCAGGAATCGCAAACCACGGTGTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGaa                                                                       <  2:81458/149‑1 (MQ=255)
                                                                                cTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCa                                                                    >  1:102589/1‑148 (MQ=255)
                                                                                  ttGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCACAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTg                                                                                                    >  1:41035/1‑114 (MQ=255)
                                                                                  ttGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCACAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTg                                                                                                    <  2:41035/114‑1 (MQ=255)
                                                                                     cccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAAc                                                               <  1:7813/148‑1 (MQ=255)
                                                                                     cccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAAc                                                               <  7:328547/148‑1 (MQ=255)
                                                                                      ccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc                                                                                                                       >  1:243565/1‑91 (MQ=255)
                                                                                      ccTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc                                                                                                                       <  2:243565/91‑1 (MQ=255)
                                                                                          cAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCAt                                                           <  3:123549/147‑1 (MQ=255)
                                                                                           aGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACg                                                        >  1:333534/1‑149 (MQ=255)
                                                                                            gCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACgg                                                       <  8:132196/149‑1 (MQ=255)
                                                                                                     aGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCa                                                                                                                                    <  1:19658/63‑1 (MQ=255)
                                                                                                     aGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCa                                                                                                                                    >  2:19658/1‑63 (MQ=255)
                                                                                                      gTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGt                                             >  2:271441/1‑149 (MQ=255)
                                                                                                      gTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGt                                             <  5:25250/149‑1 (MQ=255)
                                                                                                       tGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGAAAAAACGCATACGGGTAAAAGCGtt                                            <  5:26185/149‑1 (MQ=255)
                                                                                                         aCCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTgg                                          <  8:81508/149‑1 (MQ=255)
                                                                                                         aCCAAAGGCGATGCTGTATTCTGCAGCGACAGGGCCAGGAACCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGAGTTg                                           <  2:213098/148‑1 (MQ=255)
                                                                                                           cAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGAAGGGGCCTCTTCCGGCGATTCTTTGCTGTACCTATCCAGTTGACCGTTCTGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGt                                         >  5:72101/1‑148 (MQ=255)
                                                                                                            aaaGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCCCGGTTTCAGTGGGGCAGGGGGCTCTTCCGGGGATGCTTTGGTGTAGATATGCAGTTGAGCGGTCGGGAAGCAACAACGCATACGGGTAATCGCGgtggtgc                                       >  4:260164/1‑148 (MQ=255)
                                                                                                            aaaGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGa                                       >  8:328613/1‑149 (MQ=255)
                                                                                                              aGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGAt                                      >  5:224460/1‑148 (MQ=255)
                                                                                                                 cGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCTCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATaaaa                                  <  5:269374/149‑1 (MQ=255)
                                                                                                                 cGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATaaaa                                  >  7:297015/1‑149 (MQ=255)
                                                                                                                     gCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACg                                >  4:14490/1‑147 (MQ=255)
                                                                                                                      cTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCa                              <  2:102589/148‑1 (MQ=255)
                                                                                                                        gTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAATCGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTc                           >  3:17212/1‑149 (MQ=255)
                                                                                                                         tATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCt                          <  8:97980/149‑1 (MQ=255)
                                                                                                                             cttcAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTcc                       <  8:211140/148‑1 (MQ=255)
                                                                                                                             cttcAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTcc                       <  2:333534/148‑1 (MQ=255)
                                                                                                                             cttcAGCGACAGGGCCAGGAATCACAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTcc                       <  7:73894/148‑1 (MQ=255)
                                                                                                                               tcAGCGACAGGGCCAGGGATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAg                    >  6:126347/1‑149 (MQ=255)
                                                                                                                                    gACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGa                                                                        <  2:196237/92‑1 (MQ=255)
                                                                                                                                          gCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAGCCCCCGCAAtt         <  8:297015/149‑1 (MQ=255)
                                                                                                                                             aGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAGCCCCCGCAATTCCg      <  1:135577/149‑1 (MQ=255)
                                                                                                                                                 aTCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAGCCCCCGCAATTCCGGTGa  <  6:224460/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTTATCTTCCCAGCGCAGGCAGGTTAATGTACCACCCCAGCAGCAGCCGGTATCCAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCGGCGATTCTTTGCTGTACATATCCAGTTGACCGTTCGGGAAGCAAAAACGCATACGGGTAAAAGCGTTGGTGATAAAACGCAGTCTTCCCAGCCCCCGCAATTCCGGTGA  >  NC_000913/50440‑50733

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: