Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 638,614 G→T T14K (ACA→AAA)  thiD ← phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730638,6140GT98.7% 254.1 / ‑3.1 76T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (35/40);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  CP000730/638477‑638751
                                                                                                                                         |                                                                                                                                          
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTCTCCATGGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGATTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                         >  2:521445/1‑141 (MQ=255)
 aTCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                                        <  2:420988/141‑1 (MQ=255)
     aTTGGTAATGGTGCAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGaa                                                                                                                                    <  1:532992/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                  <  1:84048/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAATTCACGGGACCGTGCCCGACGTATCTAATTCTTGGAACTTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                  <  1:302173/141‑1 (MQ=255)
         gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCgg                                                                                                                                <  2:22330/141‑1 (MQ=255)
          tAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGc                                                                                                                               >  2:160440/1‑141 (MQ=255)
            aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                             >  1:157100/1‑141 (MQ=255)
                tgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGt                                                                                                                         <  2:183542/141‑1 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                        >  1:373514/1‑141 (MQ=255)
                  tAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑a                                                                                                                      <  1:76780/141‑1 (MQ=255)
                   aaCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑aa                                                                                                                     <  2:800194/141‑1 (MQ=255)
                    aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑aaa                                                                                                                    <  2:240224/141‑1 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAAct                                                                                                                 <  1:388039/141‑1 (MQ=255)
                         gtgtGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActtt                                                                                                               >  2:844951/1‑141 (MQ=255)
                         gtgtGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActtt                                                                                                               >  2:260768/1‑141 (MQ=255)
                          tgtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttc                                                                                                              <  2:367557/141‑1 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttct                                                                                                             >  1:856999/1‑141 (MQ=255)
                            tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAActttctt                                                                                                            >  2:362606/1‑141 (MQ=255)
                               ctaTGTATCTTTATCCATAGTAACGATGGCCGTTAAAGCGACCATGCCATACGTATCTATTTCTTGGAACTTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTaa                                                                                                         <  2:402990/139‑1 (MQ=255)
                                   gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCACAACCGGCAATTGTT‑AAAACTTTCTTTAAAGcc                                                                                                     >  1:595238/1‑141 (MQ=255)
                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGCACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCATAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   >  2:711927/1‑141 (MQ=255)
                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   >  2:804317/1‑141 (MQ=255)
                                     aTCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCAt                                                                                                   <  2:263181/141‑1 (MQ=255)
                                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGa                                                                                                >  2:157794/1‑141 (MQ=255)
                                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGa                                                                                                <  1:863898/141‑1 (MQ=255)
                                              cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCa                                                                                          >  1:344488/1‑141 (MQ=255)
                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACa                                                                                  >  2:109061/1‑141 (MQ=255)
                                                        aTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACata                                                                                >  1:129248/1‑141 (MQ=255)
                                                         tGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGGTTTCAAATCTGCTTGCATACCTGAGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataa                                                                               >  2:296230/1‑141 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                            >  2:27324/1‑141 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                            >  1:832602/1‑141 (MQ=255)
                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                         <  2:597160/141‑1 (MQ=255)
                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                         >  2:261044/1‑141 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  2:583312/141‑1 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  2:583342/141‑1 (MQ=255)
                                                                   aGCGACCATGCCATACGTATCCAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                                                                     >  2:160881/1‑141 (MQ=255)
                                                                     cGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTAtt                                                                   <  2:703008/141‑1 (MQ=255)
                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtc                                                                  <  1:504672/141‑1 (MQ=255)
                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACGTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtc                                                                  >  2:733215/1‑141 (MQ=255)
                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAGCCGGCAATTGTG‑AAAACTTTCTGTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:736420/1‑141 (MQ=255)
                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTaaaaca                                                                >  1:91101/1‑136 (MQ=255)
                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:11065/1‑141 (MQ=255)
                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAAGTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                >  1:86333/1‑141 (MQ=255)
                                                                            gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                            <  2:4753/141‑1 (MQ=255)
                                                                              cATACGTATCTAATTCTTGGGACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTATAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  2:544948/141‑1 (MQ=255)
                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  1:109061/141‑1 (MQ=255)
                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  2:12021/141‑1 (MQ=255)
                                                                                     aTCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttttt                                                    <  2:75572/141‑1 (MQ=255)
                                                                                      tCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                  >  2:491782/1‑141 (MQ=255)
                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                 >  1:290595/1‑141 (MQ=255)
                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                 <  2:679943/141‑1 (MQ=255)
                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                 <  2:196704/141‑1 (MQ=255)
                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                         >  2:370836/1‑141 (MQ=255)
                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                         <  1:469116/141‑1 (MQ=255)
                                                                                               tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                         <  1:469149/141‑1 (MQ=255)
                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       <  1:14342/141‑1 (MQ=255)
                                                                                                 gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       <  2:212122/141‑1 (MQ=255)
                                                                                                   aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaaa                                     >  2:670800/1‑141 (MQ=255)
                                                                                                               cTGCTTGCATACCTGCTCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTc                         <  1:355985/141‑1 (MQ=255)
                                                                                                                 gCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCaa                       >  1:619242/1‑141 (MQ=255)
                                                                                                                 gCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCaa                       <  2:680942/141‑1 (MQ=255)
                                                                                                                    tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatat                    <  1:416308/141‑1 (MQ=255)
                                                                                                                    tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatat                    <  1:868243/141‑1 (MQ=255)
                                                                                                                          ccTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatg              <  2:659531/141‑1 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga             <  1:453676/141‑1 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga             >  2:625313/1‑141 (MQ=255)
                                                                                                                            tacgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat            >  1:196561/3‑141 (MQ=255)
                                                                                                                            tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat            >  2:448278/1‑141 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          >  2:249265/1‑141 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          >  1:289998/1‑141 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          <  2:518804/141‑1 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          <  1:358508/141‑1 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt         <  1:491686/141‑1 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa    >  2:290413/1‑141 (MQ=255)
                                                                                                                                    cACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa    >  1:744246/1‑141 (MQ=255)
                                                                                                                                      cTTTTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATg  <  1:739034/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTT‑AAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  CP000730/638477‑638751

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