Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,500,487 T→A 5.1% N3K (AAT→AAA U712_07720 → Inositol‑1‑monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,500,4870TA 5.1% 255.2 / 3.7 79N3K (AAT→AAAU712_07720Inositol‑1‑monophosphatase
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base A (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TGAAAACGAACCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCAATATTGACAAAGAAACC  >  CP006881/1500351‑1500625
                                                                                                                                        |                                                                                                                                          
tGAAAACGAACCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTg                                                                                                                                          >  3:47006/1‑139 (MQ=255)
tGAAAACGAACCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTg                                                                                                                                          >  1:263573/1‑139 (MQ=255)
  aaaaCGAACCGCATGAATGGTTCGTTTTTACGTGAGATGGGATATTTTACAATTTCCTTGGTTTTTCATTGCGGTTATCTGTATCCTATAAAATAAAAGGTCTACAGAGAAGAGAGGAACGAACGAATGACAAATTGGa                                                                                                                                        >  7:175572/1‑139 (MQ=255)
   aaaCGAACCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGa                                                                                                                                        <  5:27318/138‑1 (MQ=255)
         aCCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAAt                                                                                                                                 >  1:279334/1‑139 (MQ=255)
                   tGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGAtt                                                                                                                       >  8:33629/1‑139 (MQ=255)
                     gTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTg                                                                                                                      >  8:97833/1‑138 (MQ=255)
                     gTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTg                                                                                                                      >  2:102064/1‑138 (MQ=255)
                     gTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTg                                                                                                                      >  1:90602/1‑138 (MQ=255)
                      ttCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCa                                                                                                                    <  7:97833/139‑1 (MQ=255)
                      ttCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCa                                                                                                                    <  3:133349/139‑1 (MQ=255)
                      ttCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCa                                                                                                                    <  3:103210/139‑1 (MQ=255)
                      gtCGTTTTTACGTGAGTATGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCa                                                                                                                    <  8:240878/138‑1 (MQ=255)
                             tttcGTGAGAAGGTTTATTTTTCATTTTCCTTGCTTTTTCATTGGCGTTATCAGTATCCTATAAAATAATATGTCTACATAGAAGAGAGGTTCGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGaaa                                                                                                              <  4:219583/135‑1 (MQ=255)
                                cGTGAGAAGGGATATTTTTGATTTTCCTTGCTTTTTCATTGCCGTTATCATTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGa                                                                                                          <  4:2851/139‑1 (MQ=255)
                                  tGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTTGACGGAAATCTATGAGATTGCAAAGAAATGGATa                                                                                                        <  8:219024/139‑1 (MQ=255)
                                                 ttCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  5:285525/138‑1 (MQ=255)
                                                 ttCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  1:29258/138‑1 (MQ=255)
                                                 ttCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGGACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  1:1626/138‑1 (MQ=255)
                                                 ttCATTTTCCTTGCTTTTTCATTGCCGTTACCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  7:321596/138‑1 (MQ=255)
                                                 ttCATTTTCCTTGCTTTTTCATTGCCGATATCATTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  6:86301/138‑1 (MQ=255)
                                                   cATTTTCCTTGCTTTTTCATTGCCGTTATAAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATAGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCgag                                                                                       >  2:123429/1‑139 (MQ=255)
                                                     ttttCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     <  6:129236/139‑1 (MQ=255)
                                                     ttttCCTTGCTTTTTCATTGCCGTTATCAGTATACTATAAAATAAAATGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAAAGGATAAGAGAAGCAGGCGCGAGaa                                                                                     <  1:144745/139‑1 (MQ=255)
                                                      tttCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     <  2:90602/138‑1 (MQ=255)
                                                          cTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTaca                                                                                <  8:127994/139‑1 (MQ=255)
                                                             gCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACaa                                                                             >  8:337180/1‑139 (MQ=255)
                                                             gCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACaa                                                                             >  5:297440/1‑139 (MQ=255)
                                                              cTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAAt                                                                            <  2:33675/139‑1 (MQ=255)
                                                                         ccGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTaca                                                                                <  2:215054/124‑1 (MQ=255)
                                                                         ccGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTaca                                                                                >  1:215054/1‑124 (MQ=255)
                                                                           gTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAg                                                               <  4:47006/139‑1 (MQ=255)
                                                                           gTTATCACTATCCAATAAAATAAAATGTCTACATAGACGAGAGGTACGAACGATAGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATTCATGAAAg                                                               <  2:263573/139‑1 (MQ=255)
                                                                                   tATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  1:137979/104‑1 (MQ=255)
                                                                                   tATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          >  2:137979/1‑104 (MQ=255)
                                                                                   tATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     >  1:87981/1‑109 (MQ=255)
                                                                                   tATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     <  2:87981/109‑1 (MQ=255)
                                                                                          taaaataaaaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  8:312680/97‑1 (MQ=255)
                                                                                          taaaataaaaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          >  7:312680/1‑97 (MQ=255)
                                                                                          taaaataaaaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAAc                                                >  5:345270/1‑139 (MQ=255)
                                                                                          taaaataaaaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAAc                                                >  2:138392/1‑139 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATg                                                                                                            >  6:204975/1‑71 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATg                                                                                                            <  5:204975/71‑1 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          >  8:98627/1‑89 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGgcgc                                                                                          <  7:98627/89‑1 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     >  5:243518/1‑94 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGaa                                                                                     <  6:243518/94‑1 (MQ=255)
                                                                                                  aaGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGaaatca                                         >  6:272898/1‑138 (MQ=255)
                                                                                                          cATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGa                                  <  4:49874/137‑1 (MQ=255)
                                                                                                                agagagGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTa                          <  1:102064/139‑1 (MQ=255)
                                                                                                                  agagGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        <  1:138392/139‑1 (MQ=255)
                                                                                                                  agagGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        >  7:311616/1‑139 (MQ=255)
                                                                                                                  agagGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        <  6:345270/139‑1 (MQ=255)
                                                                                                                  agagGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        <  7:295789/139‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTaca                                                                                >  1:151889/1‑81 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTaca                                                                                <  2:151889/81‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGa                                  <  4:21781/127‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGa                                  >  3:21781/1‑127 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        >  7:193020/1‑137 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        <  8:193020/137‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        >  2:103411/1‑137 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGt                        <  1:103411/137‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCa                      >  4:122847/1‑139 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCa                      >  3:292179/1‑139 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCa                      >  6:185031/1‑139 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCa                      >  1:172571/1‑139 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCCGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCa                      >  2:40421/1‑139 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAAATGGACGGAAATCGATGAGATTGCAAAGAAATGGATa                                                                                                        >  1:299887/1‑57 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAAATGGACGGAAATCGATGAGATTGCAAAGAAATGGATa                                                                                                        <  2:304488/57‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAAATGGACGGAAATCGATGAGATTGCAAAGAAATGGATa                                                                                                        <  2:299887/57‑1 (MQ=255)
                                                                                                                    agGTACGAACGAATGACAAAATGGACGGAAATCGATGAGATTGCAAAGAAATGGATa                                                                                                        >  1:304488/1‑57 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATg                                                                                                            <  8:281371/48‑1 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATg                                                                                                            >  7:281371/1‑48 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGaaat                                           >  5:146183/1‑113 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGaaat                                           <  6:146183/113‑1 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCaa                  >  3:299407/1‑138 (MQ=255)
                                                                                                                         cgaacgaaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCaa                  >  3:46830/1‑138 (MQ=255)
                                                                                                                               aaTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCAATATTGa            >  8:139167/1‑138 (MQ=255)
                                                                                                                                aTGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCAATATTGAc           >  1:181064/1‑138 (MQ=255)
                                                                                                                                        ttGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCAATATTGACAAAGAAAcc  >  6:192709/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGAAAACGAACCGCATGAATGGTTCGTTTTTACGTGAGAAGGGATATTTTTCATTTTCCTTGCTTTTTCATTGCCGTTATCAGTATCCTATAAAATAAAAGGTCTACATAGAAGAGAGGTACGAACGAATGACAAATTGGACGGAAATCGATGAGATTGCAAAGAAATGGATAAGAGAAGCAGGCGCGAGAATTACACAATCCATGCATGAAAGCCTGACAATTGAAACGAAATCAAATCCGAATGATTTAGTCACCAATATTGACAAAGAAACC  >  CP006881/1500351‑1500625

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