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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730638,6140GT99.1% 219.3 / ‑6.2 69T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (33/34);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  CP000730/638477‑638751
                                                                                                                                         |                                                                                                                                         
cATCCATTGGTAATTGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                        <  1:461409/141‑1 (MQ=255)
 aTCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGGAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                                       <  2:263181/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                 <  1:551351/141‑1 (MQ=255)
       tGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                 <  1:274851/141‑1 (MQ=255)
         gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCGCTTTTGTCGGAGCCgg                                                                                                                               >  1:8952/1‑141 (MQ=255)
              ggTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATGCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAAAACCGGCAAtg                                                                                                                          >  2:424891/1‑140 (MQ=255)
               gtgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTg                                                                                                                         >  1:326961/1‑141 (MQ=255)
                tgtAACATCGTGTGACCATGTATCTTTATCCATATCCACGACGGCAGTTAAAGCGACCATGCCATACGTATCCAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGt                                                                                                                        <  1:514101/141‑1 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCTCTTTTGTCAGAACCGGCAATTGtt                                                                                                                       <  1:487567/141‑1 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTCTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                       <  1:248110/141‑1 (MQ=255)
                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAc                                                                                                                  >  1:408163/1‑141 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGGAATTGTTAAAAct                                                                                                                 >  1:319657/1‑141 (MQ=255)
                        cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGTAATTGTTAAAActt                                                                                                                <  2:453679/141‑1 (MQ=255)
                        cGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                                >  1:166781/1‑141 (MQ=255)
                        cGTGTGACCATGTATCTTTATCCAGAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                                <  2:25965/141‑1 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttct                                                                                                             <  2:180392/141‑1 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttct                                                                                                             >  1:525982/1‑141 (MQ=255)
                           gtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttct                                                                                                             >  1:525862/1‑141 (MQ=255)
                            tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttctt                                                                                                            <  2:519896/141‑1 (MQ=255)
                              aCCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                                                                          >  2:105413/1‑141 (MQ=255)
                                 aTGTATCTTTATCCATAGTAGCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAg                                                                                                       <  2:654075/141‑1 (MQ=255)
                                   gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTGTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                                                                     <  1:244059/141‑1 (MQ=255)
                                      tCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTACTTGGGTACTGGCTCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAtt                                                                                                  >  1:655125/1‑141 (MQ=255)
                                        tttATCCATAGTAGCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGa                                                                                                <  2:202502/141‑1 (MQ=255)
                                         ttATCCATAGTACCGATGGCAGTTCAAGCGCCCATGCCATACGTATATAATTCTTGGAACGTTTCCCAATCTGCTGGTATACTGGCGCCAGCAGTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                                                               <  2:387870/141‑1 (MQ=255)
                                         ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                                                               <  1:214283/141‑1 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCGTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                              <  1:340487/141‑1 (MQ=255)
                                            tCCATAGTAACGATGGCAGTTGACACGACCATGCCTTACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCtt                                                                                            <  2:498690/141‑1 (MQ=255)
                                              cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCa                                                                                          >  2:658257/1‑141 (MQ=255)
                                                gaGTAACGATGGCAGTTAAAGCGACCATGCCATACGTAGCTAATTCTTGGAACGTTTTCAAATCTGATTGCTTACTTGCGCCAGCTCTTTTGTCACAACCGGCAATAGTTACAACTTTCTTTAAAGCCATTCAGCTTCACt                                                                                        >  2:444513/2‑141 (MQ=255)
                                                   tAACGATGGCAGTTAACGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                     <  2:591306/141‑1 (MQ=255)
                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACa                                                                                  >  1:466275/1‑141 (MQ=255)
                                                       gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACat                                                                                 <  2:428640/141‑1 (MQ=255)
                                                           gCAGTTAAAGCGACCATGCCATACGTATCTAATTATTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTCCctaata                                                                             >  2:483527/1‑141 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                            >  2:292216/1‑141 (MQ=255)
                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATCatat                                                                         >  2:66172/1‑141 (MQ=255)
                                                               ttAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                         <  2:209548/141‑1 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                        <  2:363604/141‑1 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACAGAATAATAtt                                                                        >  1:415802/1‑141 (MQ=255)
                                                                        ccATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcat                                                                <  2:202966/141‑1 (MQ=255)
                                                                          aTGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATACTATTGTATtcatca                                                              >  1:396032/1‑141 (MQ=255)
                                                                           tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                             >  1:158819/1‑141 (MQ=255)
                                                                            gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                            >  2:356359/1‑141 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTATATAATAATATTGTATTCatcatat                                                           <  1:427742/141‑1 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                           >  1:564108/1‑141 (MQ=255)
                                                                              cATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCATCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatc                                                          <  1:196115/141‑1 (MQ=255)
                                                                                      tCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                  <  1:450711/141‑1 (MQ=255)
                                                                                                ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAGTTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCACATCATTTTCAACCTAATTGaa                                        >  2:608072/1‑141 (MQ=255)
                                                                                                ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaa                                        >  2:206440/1‑141 (MQ=255)
                                                                                                 gAACGTTTTCAAAGCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAAAAATATTGTATTCATCATATCATTTTTAACCTAATTGaaa                                       <  1:157980/141‑1 (MQ=255)
                                                                                                  aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaa                                      >  1:510925/1‑141 (MQ=255)
                                                                                                    cGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTAAACCCTGACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAat                                    >  1:455411/1‑141 (MQ=255)
                                                                                                       tttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                 >  2:158277/1‑141 (MQ=255)
                                                                                                        ttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTa                                <  1:358047/141‑1 (MQ=255)
                                                                                                          cAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAg                              >  1:327859/1‑141 (MQ=255)
                                                                                                               cTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCCTTTTTAACCTAATTGAAAAATATTAAGCATTc                         >  1:586662/1‑141 (MQ=255)
                                                                                                                tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCa                        >  1:605271/1‑141 (MQ=255)
                                                                                                                  cTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAat                      <  2:284432/141‑1 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga             <  1:292216/141‑1 (MQ=255)
                                                                                                                             gcgcAAGCACTTTTGTCAGAACCTGCAATTGTTAAATCTTTCTTTAAAGCCTTTGAGCTTCACTCCTACATAATAATCTTGTATTCATCTTATCATTTTTAACCTAATGGAAAAATATTAAGCATTCAATATTTGATGAtt           <  2:437999/141‑1 (MQ=255)
                                                                                                                             gagcgatcACTTTTGTCAGAACGAGAAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAtt           <  2:109193/134‑1 (MQ=255)
                                                                                                                              cgcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg          >  2:360596/1‑141 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt         <  1:147952/141‑1 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt         <  2:657449/141‑1 (MQ=255)
                                                                                                                                 cAGCACTGGTGTCAGAACCGGCAATTGTAAAAACTTTCTTTAAAGCCAGTGAGCTTCACTCCGACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgttg       <  2:615751/141‑1 (MQ=255)
                                                                                                                                   gCACTTTTGTCAGAACGGGCAATTGTTAAAACTTTCTTAAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGCAAAATATTAAGCATTCAATATTTGATGATTGTTGaa     >  1:371914/1‑141 (MQ=255)
                                                                                                                                   gCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTGGaa     >  2:520236/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa    >  1:566141/1‑141 (MQ=255)
                                                                                                                                     ccTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAt   <  2:408163/140‑1 (MQ=255)
                                                                                                                                      cTTGTGTCATAACCGGCAATTGTTAAACCTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATg  <  2:224992/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATG  >  CP000730/638477‑638751

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