Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,565,986 A→G F421L (TTC→CTC)  yjiN ← zinc‑type alcohol dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,565,9860AG99.0% 360.6 / NA 96F421L (TTC→CTC) yjiNzinc‑type alcohol dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (48/47);  total (49/47)
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.99e-01

TGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCCGCGACATATCCCGCG  >  NC_000913/4565850‑4566114
                                                                                                                                        |                                                                                                                                
tGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCa                                                                                                                                <  1:1508928/139‑1 (MQ=255)
  aTTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAg                                                                                                                              <  1:509465/139‑1 (MQ=255)
  aTTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAg                                                                                                                              >  1:488021/1‑139 (MQ=255)
  aTTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAg                                                                                                                              >  1:186182/1‑139 (MQ=255)
     ttccttccCCGTACCACGGGAGATAATGAAAAGAGGCCACATTAGTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGcc                                                                                                                           >  1:1261112/1‑139 (MQ=255)
      tccttccCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCg                                                                                                                          <  2:1246176/139‑1 (MQ=255)
            cccGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCt                                                                                                                    >  1:1532314/1‑139 (MQ=255)
                    aCGGGCGATAATGAAAAGAGGCCACATTATTACCATTCTGTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTGAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGAGAGc                                                                                                            <  2:1087077/139‑1 (MQ=255)
                    aCGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            <  2:1261112/139‑1 (MQ=255)
                     cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAg                                                                                                                      <  2:335555/128‑1 (MQ=255)
                     cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAg                                                                                                                      >  1:335555/1‑128 (MQ=255)
                     cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCa                                                                                                           <  2:186182/139‑1 (MQ=255)
                     cGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCa                                                                                                           <  2:1532314/139‑1 (MQ=255)
                     cGGGAGATAATGAAAAGAGGCCACATTAGTACCATTCTTTTAGTATGGTCTATTTGATCTCGCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCa                                                                                                           >  1:1156141/1‑139 (MQ=255)
                      gggAGATAATGAAAAGAGGACACATTATTACCATTCTTTTAGCATGGTTTATTTGATCTCTCTCGCGGTTAGCCACCTAGCTTTTCAAGTATTTCTCACAATTGTCGAGTGGGAACAAGACCGTGAGCAGCGAcaacaa                                                                                                          >  1:1758760/1‑139 (MQ=255)
                       ggAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCaaa                                                                                                         >  1:1258025/1‑139 (MQ=255)
                          gaTAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAatat                                                                                                      >  1:664533/1‑139 (MQ=255)
                          gaTAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAatat                                                                                                      >  1:285585/1‑139 (MQ=255)
                                        gCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCGACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAAt                                                                                        >  2:598481/1‑139 (MQ=255)
                                         ccACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATa                                                                                       >  2:541613/1‑139 (MQ=255)
                                               tattaCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCg                                                                                 <  2:113258/139‑1 (MQ=255)
                                                  taCCATTCTTTTAGTATGTTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                                 <  1:1139781/120‑1 (MQ=255)
                                                  taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAg                                                                                                                      >  1:1703086/1‑99 (MQ=255)
                                                  taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAg                                                                                                                      <  2:1703086/99‑1 (MQ=255)
                                                  taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            <  2:446090/109‑1 (MQ=255)
                                                  taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            >  1:446090/1‑109 (MQ=255)
                                                  taCCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAAt                                                                                                 >  2:1139781/1‑120 (MQ=255)
                                                    ccATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  2:285585/139‑1 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAata                                                                                                       >  1:1014483/1‑109 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAata                                                                                                       <  2:1014483/109‑1 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  1:1524589/136‑1 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  2:1524589/1‑136 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCa                                                                          <  1:867610/138‑1 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCa                                                                          >  2:867610/1‑138 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                                         >  2:1386748/1‑139 (MQ=255)
                                                       ttCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAg                                                                         >  2:1683047/1‑139 (MQ=255)
                                                                  tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                                              <  2:51706/139‑1 (MQ=255)
                                                                  tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                                              <  1:1683047/139‑1 (MQ=255)
                                                                  tGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                                              >  1:143481/1‑139 (MQ=255)
                                                                          ttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                                      <  2:143481/139‑1 (MQ=255)
                                                                          ttGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                                      <  2:664533/139‑1 (MQ=255)
                                                                              tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCa                                                                          >  2:1509161/1‑115 (MQ=255)
                                                                              tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCa                                                                          <  1:1509161/115‑1 (MQ=255)
                                                                              tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                                                  >  1:117888/1‑139 (MQ=255)
                                                                              tctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACt                                                  >  2:1490895/1‑139 (MQ=255)
                                                                               ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCa                                                                                                                                <  2:1050332/60‑1 (MQ=255)
                                                                               ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  2:878155/112‑1 (MQ=255)
                                                                               ctctctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  1:878155/1‑112 (MQ=255)
                                                                                  tctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  2:942805/109‑1 (MQ=255)
                                                                                  tctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  1:942805/1‑109 (MQ=255)
                                                                                  tctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAg                                              >  1:656571/1‑139 (MQ=255)
                                                                                   ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                                      >  2:479558/1‑98 (MQ=255)
                                                                                   ctcGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                                      <  1:479558/98‑1 (MQ=255)
                                                                                      gcgGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCt                                          <  2:860380/139‑1 (MQ=255)
                                                                                          ttAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGc                                      <  1:1490895/139‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            >  2:974569/1‑67 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            <  1:974569/67‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCa                                                                                                           >  2:1367329/1‑68 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCa                                                                                                           <  1:1367329/68‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGac                                                             >  2:1347208/1‑114 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGac                                                             <  1:1347208/114‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                                      <  1:194875/121‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATa                                                      >  2:194875/1‑121 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    <  2:1258025/139‑1 (MQ=255)
                                                                                            aGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    >  2:479368/1‑139 (MQ=255)
                                                                                              ccACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                                      <  1:1767063/87‑1 (MQ=255)
                                                                                              ccACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATAc                                                                                      >  2:1767063/1‑87 (MQ=255)
                                                                                                aCTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  2:951553/1‑95 (MQ=255)
                                                                                                aCTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  1:951553/95‑1 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            <  1:1103442/61‑1 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGc                                                                                                            >  2:1103442/1‑61 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATaaaa                                                                                                  >  1:1398647/1‑71 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATaaaa                                                                                                  <  2:1398647/71‑1 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAttt                              >  2:1537515/1‑139 (MQ=255)
                                                                                                  ttAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATAttt                              >  1:490850/1‑139 (MQ=255)
                                                                                                    aGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  1:855096/1‑91 (MQ=255)
                                                                                                    aGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  2:855096/91‑1 (MQ=255)
                                                                                                     gttgttCATGGCTTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATCAACTGCAGATCTTTGCCGATATTTAAc                           <  2:1156141/139‑1 (MQ=255)
                                                                                                           cATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    <  2:1705795/124‑1 (MQ=255)
                                                                                                           cATGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    >  1:1705795/1‑124 (MQ=255)
                                                                                                            aTGGATTTCTAAAAATTTCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGc                                                >  1:1022216/1‑111 (MQ=255)
                                                                                                            aTGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                                              >  1:990283/1‑97 (MQ=255)
                                                                                                            aTGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGa                                                              <  2:990283/97‑1 (MQ=255)
                                                                                                            aTGGATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGc                                                <  2:1022216/111‑1 (MQ=255)
                                                                                                               gATTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGcc                 <  1:598481/139‑1 (MQ=255)
                                                                                                                aTTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    <  1:953634/119‑1 (MQ=255)
                                                                                                                aTTTCTAAAAATTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    >  2:953634/1‑119 (MQ=255)
                                                                                                                          aTTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            <  2:760549/69‑1 (MQ=255)
                                                                                                                          aTTGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAAc                                                                            >  1:760549/1‑69 (MQ=255)
                                                                                                                            tGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAatat                                                                                                      >  2:119161/1‑41 (MQ=255)
                                                                                                                            tGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAatat                                                                                                      <  1:119161/41‑1 (MQ=255)
                                                                                                                            tGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGcc                 <  2:500236/126‑1 (MQ=255)
                                                                                                                            tGCCGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGcc                 >  1:500236/1‑126 (MQ=255)
                                                                                                                              ccGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCCGCGACATATCCcgcg  >  2:36924/1‑139 (MQ=255)
                                                                                                                               cGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    <  2:510114/104‑1 (MQ=255)
                                                                                                                               cGAGGGGGAGCAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCg                                    >  1:510114/1‑104 (MQ=255)
                                                                                                                                        |                                                                                                                                
TGATTTTCCTTCCCCGTACCACGGGAGATAATGAAAAGAGGCCACATTATTACCATTCTTTTAGTATGGTTTATTTGATCTCTCTCGCGGTTAGCCACTTAGTTTTTCATGGATTTCTAAAAATTGCCGAGGGGGAACAAGGCCGGGAGCTGCGACAGCAAATATAAAATTAGCCCAATACAACCGCCAACCAGCGTACCGTTGACACGAATAAACTGCAGATCTTTGCCGATATTTAACTCTATTTGCCGCGACATATCCCGCG  >  NC_000913/4565850‑4566114

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