Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,125,066 G→C 100% E729Q (GAG→CAG)  pycA → pyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,125,0660GC97.2% 214.3 / ‑2.8 68E729Q (GAG→CAG) pycApyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (37/28);  total (39/29)
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.85e-01
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.

TGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAA  >  CP000730/1124930‑1125192
                                                                                                                                        |                                                                                                                              
tGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGtt                                                                                                                            >  1:563875/1‑141 (MQ=255)
 gAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTa                                                                                                                           >  1:246898/1‑141 (MQ=255)
  aGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTaa                                                                                                                          >  2:310299/1‑141 (MQ=255)
          aaCATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAACTCAGcg                                                                                                                  >  2:427694/1‑140 (MQ=255)
          aaCATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTGATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGCTAAAATCAGCt                                                                                                                  >  2:597615/1‑141 (MQ=255)
              tttATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTGAAATCAGCTGTAg                                                                                                              >  1:199445/1‑141 (MQ=255)
                 ttaCTTTAGAGTATTATCTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCCTATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGAtt                                                                                                           <  1:379015/140‑1 (MQ=255)
                 ataCTTTAGATTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGAtt                                                                                                           <  1:631472/141‑1 (MQ=255)
                 ataCTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGAtt                                                                                                           <  2:550660/141‑1 (MQ=255)
                 ataCTTTAGAGTATTATGGCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGAtt                                                                                                           <  1:631303/141‑1 (MQ=255)
                  taCTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGAttt                                                                                                          <  2:233189/141‑1 (MQ=255)
                  taCTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACAAATTGATTGGTAAGTTAAAATCAACTGTAGAttt                                                                                                          >  2:347137/1‑141 (MQ=255)
                      ttAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCa                                                                                                      >  1:134141/1‑141 (MQ=255)
                          agTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAAtcc                                                                                                  >  1:301620/1‑139 (MQ=255)
                          agTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTc                                                                                                  >  1:429067/1‑141 (MQ=255)
                              ttatGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAAATCATCt                                                                                              >  2:87904/1‑141 (MQ=255)
                              ttatGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAAtttctct                                                                                              >  2:666174/1‑136 (MQ=255)
                              ttatGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCt                                                                                              >  2:340470/1‑141 (MQ=255)
                              ttatGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCt                                                                                              >  1:437061/1‑141 (MQ=255)
                              ttatGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATGAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTACAAACAGCTGTAGGTTTACCAATTCATCt                                                                                              >  2:433519/1‑141 (MQ=255)
                              ttatGTCAAAATAGCTAAAGAGTTAGAACGTGCAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGAGTGGTCAGTTAAAAACAGCTGTAGATTTACCAATTCAGCt                                                                                              >  2:165673/1‑141 (MQ=255)
                                     aaaCTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACt                                                                                       <  2:232694/141‑1 (MQ=255)
                                            cTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATa                                                                                <  1:409321/141‑1 (MQ=255)
                                             tAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATAc                                                                               <  1:541402/141‑1 (MQ=255)
                                                  agTTAGAACGTGAAGGTTGCCATATTTTAGCGCTTAAAGATATGGCAGGCTTATCAAAACCTAAAGCCGCTTACGAATTGATTGCTCAGTTAAAATCAGCAGTAGATCTACCAATTCATCTTCACCCTCATGATACaaatc                                                                          >  2:116851/1‑138 (MQ=255)
                                                     tAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTa                                                                       <  2:478115/141‑1 (MQ=255)
                                                       gAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAGTTCATCTTCACACTCAGCATACAAGTGGTAAt                                                                     >  2:594359/1‑141 (MQ=255)
                                                       gAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAAt                                                                     >  1:45207/1‑141 (MQ=255)
                                                         aCGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGCTTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATgg                                                                   <  2:44824/141‑1 (MQ=255)
                                                         aCGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCCCACTCATGATACAAGCGGCAATgg                                                                   >  1:246734/1‑141 (MQ=255)
                                                         aCGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATgg                                                                   <  1:314019/141‑1 (MQ=255)
                                                         aCGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATgg                                                                   <  2:193474/141‑1 (MQ=255)
                                                                 gTTTCCATATTTTGGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTaa                                                           <  2:15497/141‑1 (MQ=255)
                                                                 gTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTaa                                                           >  2:593705/1‑141 (MQ=255)
                                                                   ttCCCTATTTGAGCGATTAAAGGTTTGGCAGTCTTATTAAACCCTAAAGCCGCTTACGACTTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTACCa                                                         <  2:169193/141‑1 (MQ=255)
                                                                      cATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATac                                                      >  1:662594/1‑141 (MQ=255)
                                                                      cATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATac                                                      >  1:639736/1‑141 (MQ=255)
                                                                      cATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATac                                                      >  1:249073/1‑141 (MQ=255)
                                                                      cATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATAGATTGGTCAGTTAAAGTCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATac                                                      >  2:296738/1‑141 (MQ=255)
                                                                        tatTTTAGTGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATacaa                                                    >  1:207497/1‑141 (MQ=255)
                                                                        tatTTGAGTGATTAATGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATGGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATacaa                                                    <  2:509530/141‑1 (MQ=255)
                                                                         atTTTAGCGATTAAAGATATGGCATGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATacaaa                                                   <  1:582986/141‑1 (MQ=255)
                                                                         atTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATacaaa                                                   <  1:35933/141‑1 (MQ=255)
                                                                            ttAGCGATTAAAGATATGGCAGGGTTATTAAAACCGAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATacaaacaa                                                <  2:249325/141‑1 (MQ=255)
                                                                                cGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCaa                                            <  1:560959/141‑1 (MQ=255)
                                                                                       aGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGc                                     >  1:293123/1‑141 (MQ=255)
                                                                                        gATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCt                                    >  1:642009/1‑141 (MQ=255)
                                                                                         atatGGCAGGCTTATTAAAACCTACAGCAGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGAGGTAATGGATTATTAACATACAGACAAACAATAGatgatg                                   >  2:396653/1‑141 (MQ=255)
                                                                                            tGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGACTCGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCACGAGACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGgtg                                <  1:57570/141‑1 (MQ=255)
                                                                                             ggCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGATGgtgt                               >  2:658722/1‑141 (MQ=255)
                                                                                               cAGGCTTATTAAAACCTAAAGCCGGTTACTAATTGATTTGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAGATGGTAATGGTTTATTATCATACAAACAAGCAATAGATGCTGGTGTCg                             >  1:160876/1‑141 (MQ=255)
                                                                                               cAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCAAGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCg                             >  2:659932/1‑141 (MQ=255)
                                                                                                  gcctattaAAACCGAACGCCGCCTACGAATTGATTGGTCAGTTAAAAACAGCTGTACATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGata                          <  2:148631/138‑1 (MQ=255)
                                                                                                  gCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGata                          <  2:76992/141‑1 (MQ=255)
                                                                                                           aaCCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACACACAAGCAATAGCTGCTGGTGTCGATATCATTGATa                 >  2:163334/1‑141 (MQ=255)
                                                                                                           aaCCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTAATGATACAAGTGGTAATGGATTATTAACATACAAACAAGCAATCGATTCTGGTGTCGAAATCATTGATa                 >  2:268082/1‑141 (MQ=255)
                                                                                                            aCCTAATGCCGCTGACGAAGGTATTTGTCAGTTAAAATCGGGTGGAGGTTTACCAATTCATCTTCACACTCATTAGACTAGTGGTAATGTTTTATTAACATACAAACAAGCAATAGATGCTGGGGTCGATATCATTGATAc                <  2:199445/141‑1 (MQ=255)
                                                                                                               tAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATActgc             >  1:406900/1‑141 (MQ=255)
                                                                                                               tAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATACATCTTCACACTCATGATACAAGTGGTAATGGTTTATTACCATACAACCAAGCAATACAAGCTGGTATCGATATCATTGTTAcggc             >  1:77691/1‑141 (MQ=255)
                                                                                                                aaaGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTGCCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATActgct            >  1:434102/1‑141 (MQ=255)
                                                                                                                aaaGCCCCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACATGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATActgct            <  2:542957/141‑1 (MQ=255)
                                                                                                                       cTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGCTGCtt     >  2:644695/1‑141 (MQ=255)
                                                                                                                       cTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGATTATTAGCATACAAACAAGCAATAGATGCTAATGTAGATATCCTTGATACTGACGTTTCtt     >  2:457707/1‑141 (MQ=255)
                                                                                                                        ttACGAATTTATTGGTCAGATAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTc    <  2:357522/141‑1 (MQ=255)
                                                                                                                        ttACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTc    <  1:320056/141‑1 (MQ=255)
                                                                                                                          aCGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCaa  <  1:631414/141‑1 (MQ=255)
                                                                                                                          aCGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCaa  <  2:526493/141‑1 (MQ=255)
                                                                                                                          aCGAATTGATTGGTCAGTTAAAATCAGCTGTAGACTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCaa  <  1:251035/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                              
TGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAA  >  CP000730/1124930‑1125192

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