Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,624,673 T→A 100% K192I (AAA→ATA)  malA ← alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,624,6730TA95.7% 201.0 / 0.3 70K192I (AAA→ATA) malAalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (38/28);  total (41/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 = 8.53e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCC  >  CP000730/1624536‑1624809
                                                                                                                                         |                                                                                                                                        
aaGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTatat                                                                                                                                       >  2:372727/1‑141 (MQ=255)
  gTTTGTATACCTGGCTGATTAATATCTACATCAAATGCTGGGGCATATGTTTTAGCCTCAGGTACAGGTAAGTCACCCGGTTCAAACGTCTTCTTAATATGCATACTTGCATGTACTCTAAATCCATCAATGCCTGTATCa                                                                                                                                     >  1:430291/1‑141 (MQ=255)
    ttGTATACCTGGCTGATTCATATCTACACCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAACCGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCaaa                                                                                                                                   <  1:529871/141‑1 (MQ=255)
      gTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAacc                                                                                                                                 <  1:509771/141‑1 (MQ=255)
        ataCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAaccac                                                                                                                               >  1:204633/1‑141 (MQ=255)
        ataCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGGGGTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATTCCTATATCAAaccac                                                                                                                               >  1:520473/1‑141 (MQ=255)
         taCCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAacaacc                                                                                                                              >  1:247527/1‑141 (MQ=255)
                  gATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTtcatca                                                                                                                     <  1:372727/141‑1 (MQ=255)
                   aTTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTGAATCCATCAATGCCTATATGATACCACCAGTTcatcat                                                                                                                    >  1:21032/1‑141 (MQ=255)
                   aTTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTcatcat                                                                                                                    <  1:312112/141‑1 (MQ=255)
                    ttCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCAtt                                                                                                                   >  2:499058/1‑141 (MQ=255)
                     tCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCAttt                                                                                                                  >  2:91602/1‑141 (MQ=255)
                        tatCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCaa                                                                                                               <  2:611640/141‑1 (MQ=255)
                          tCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAAt                                                                                                             >  1:264686/1‑141 (MQ=255)
                          tCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAAt                                                                                                             >  1:665114/1‑141 (MQ=255)
                          tCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGGACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAAt                                                                                                             >  1:440217/1‑141 (MQ=255)
                              cATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAg                                                                                                         <  2:169154/141‑1 (MQ=255)
                              cATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAg                                                                                                         <  1:343081/141‑1 (MQ=255)
                                     tGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACCGccacact                                                                                                  >  2:495673/1‑135 (MQ=255)
                                      gCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCTTCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTa                                                                                                 >  2:2607/1‑141 (MQ=255)
                                        tGGGGCATATGTTTTACCCTCAGGTACAGGGAAGTCACCCGCTTCAAACGTATTCTTACTATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAAc                                                                                               >  1:454242/1‑141 (MQ=255)
                                         ggggCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACt                                                                                              >  1:28402/1‑141 (MQ=255)
                                          gggCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACtt                                                                                             <  1:444318/141‑1 (MQ=255)
                                             cATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCg                                                                                          <  1:273844/141‑1 (MQ=255)
                                               tatGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCACCCCTAACTTCCGGa                                                                                        >  2:325325/1‑141 (MQ=255)
                                                  gTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATAGGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTa                                                                                     <  1:506164/141‑1 (MQ=255)
                                                   ttttACCCTCAGGTACAGGTAAGTCACCAGCGTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAc                                                                                    >  1:659586/1‑141 (MQ=255)
                                                      tACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTCCTCTAAATCCATCAATGCCTATATCAAACCACAAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAcccc                                                                                 <  2:641687/141‑1 (MQ=255)
                                                      tACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTAcccc                                                                                 >  1:581173/1‑141 (MQ=255)
                                                          cTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAATCCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAAtt                                                                             <  1:431912/141‑1 (MQ=255)
                                                          cTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCCAtt                                                                             >  1:545/1‑141 (MQ=255)
                                                          cTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAAATTCCGGATTACCCCAAtt                                                                             >  1:392101/1‑141 (MQ=255)
                                                          cTCAGGTACAGGTAAGTCACCCGCGTCAAACGTCTTCTTAAGATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAAtt                                                                             <  2:101487/141‑1 (MQ=255)
                                                            caggtacaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCa                                                                           >  1:443927/1‑141 (MQ=255)
                                                              ggtacaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACACCAATTCaaa                                                                         >  2:31569/1‑141 (MQ=255)
                                                                 acaggtaAGTCACCTGCTTCAAAGGTCTTCTTAATATGCGTAATTGCATCTACTTTAAATCCATCAATGCCTATATCAAACCACAAGTTCATCATTTCAAATACAGCATCTTTAACTTCTGGATTACTCGATTTAAAATcg                                                                      >  2:295416/1‑140 (MQ=255)
                                                                 acaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCa                                                                      >  1:449557/1‑141 (MQ=255)
                                                                 acaggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATACATCAATGCCTATATCAAACCACCAGTTCATCATTTCACATACAGCATCTCTAACTTCCGGATTACTCCAATTCAAATCa                                                                      >  1:402034/1‑141 (MQ=255)
                                                                    ggtaAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGgt                                                                   >  2:238592/1‑141 (MQ=255)
                                                                              ccGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                         <  2:633188/141‑1 (MQ=255)
                                                                              ccGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                         >  1:386472/1‑141 (MQ=255)
                                                                              ccGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                         <  2:647564/141‑1 (MQ=255)
                                                                              ccGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                         >  1:53889/1‑141 (MQ=255)
                                                                              ccGCGTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTAc                                                         <  2:193566/141‑1 (MQ=255)
                                                                               cGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACt                                                        >  1:599118/1‑141 (MQ=255)
                                                                                ccTTAAAACGTCTTCTTAATTTGCTTAAGTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATGCAAATCAGGTTGTTTTTTACTg                                                       <  2:101527/140‑1 (MQ=255)
                                                                                 cTTCAAACGTCTTCATAATCTGCGTAATTGCATCTACTCTAAATCCATCATTGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGa                                                      <  1:145376/141‑1 (MQ=255)
                                                                                      aaCGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATaaa                                                 <  1:278525/141‑1 (MQ=255)
                                                                                        cGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATg                                               >  1:573774/1‑141 (MQ=255)
                                                                                        cGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATg                                               <  1:380887/141‑1 (MQ=255)
                                                                                        cGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATg                                               >  2:254679/1‑141 (MQ=255)
                                                                                        cGTCTTCTTAATATGCGTAATTGCAGCTACTCTAAATCCATCATTGCCGATATCAAACCACCAGTTCATAATTTCCAATACAGCAACTCTAACTCCCCGATTTCCCCAATTAAAATCAGGTTGTTTTTTACTGCATAAATg                                               >  2:290969/1‑141 (MQ=255)
                                                                                         ttcttctCAATACGCACAACTGCATCTGTTCTCACCCCCTCAATGCCCATTCCCAACCCCCCGTTCGCCAGTTCACATTCACCCTTTCTAACTGCCGGATTACCCCAATTCAAATCACGTTGTTCTGTACTGAATacttgc                                              <  2:616781/140‑6 (MQ=255)
                                                                                                tgagatGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAACCCACCAGTTCAGCATTTCAAATACAGCATCTCTAACTTCCGGAGTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAaataat                                       <  2:214459/137‑1 (MQ=255)
                                                                                                  atatGCGTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAatat                                     >  2:364892/1‑141 (MQ=255)
                                                                                                       ggTAATTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTc                                <  2:454242/140‑1 (MQ=255)
                                                                                                       cGTAATTGCACCTACTCTAACTCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTc                                <  2:382946/141‑1 (MQ=255)
                                                                                                          aaTTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGTTCAGt                             >  2:426306/1‑141 (MQ=255)
                                                                                                          aaTTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCAGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGt                             >  2:553701/1‑141 (MQ=255)
                                                                                                           aTTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCACTTCTAATCAGGTTCTTTTTTACTGCATAAATGGGACTAATATTGCTCAGTa                            >  2:293666/1‑141 (MQ=255)
                                                                                                           aTTGCATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGtt                            >  1:176249/1‑140 (MQ=255)
                                                                                                              gCATCTACTCTAAATCCATCAATGCCTATATCACAACACCAGTTCATCACTTAAAATAAAGCATCTCTAACTTCCGGATTACCCCAACACAAATCAGGTGGTGTGTTACTTTAGAACTGGAAATAATATCGCTCAGGATTa                         >  2:554942/1‑141 (MQ=255)
                                                                                                               cATCTACTCTAAATCCATCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAg                        >  2:1795/1‑141 (MQ=255)
                                                                                                                aTCTACTCTAAATCCCTCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATGTTGCTCAGTATTAGc                       <  1:494347/141‑1 (MQ=255)
                                                                                                                      cccaactCCATGAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACCGCATCTCTAACTTCCGGATTACCCCACTTCAAATACGGGGTCTGTTGACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCata                 <  2:640458/135‑1 (MQ=255)
                                                                                                                               aTCAATGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGt        <  1:429311/141‑1 (MQ=255)
                                                                                                                                tCAAGGCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTTCTCAGTATTAGCATCATATTCCCATGTa       <  1:429168/141‑1 (MQ=255)
                                                                                                                                 cgacgCCTATATCAAACCACCAGTTCTTCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATCCCTATGTAg      <  1:488830/137‑1 (MQ=255)
                                                                                                                                     gCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCACTTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATcc  >  2:275562/1‑141 (MQ=255)
                                                                                                                                     gCCTATATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATcc  >  1:422943/1‑141 (MQ=255)
                                                                                                                                     gCCTATATCAAACCACCAGTTCATCATTTCAAATAAAGCATCTGTAAATTCCGGATTACCCCAATTCAAATCAGATTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCGCTGGTAGAGcc  >  1:505022/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                        
AAGTTTGTATACCTGGCTGATTCATATCTACATCAAATGCTGGGGCATATGTTTTACCCTCAGGTACAGGTAAGTCACCCGCTTCAAACGTCTTCTTAATATGCGTAATTGCATCTACTCTAAATCCATCAATGCCTTTATCAAACCACCAGTTCATCATTTCAAATACAGCATCTCTAACTTCCGGATTACCCCAATTCAAATCAGGTTGTTTTTTACTGAATAAATGGAAATAATATTGCTCAGTATTAGCATCATATTCCCATGTAGATCC  >  CP000730/1624536‑1624809

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: