Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,152 (C)5→6 coding (816/2463 nt) thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,1471.C97.0% 232.6 / ‑0.3 67coding (811/2463 nt)thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (30/35);  total (31/36)
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.99e-01

GACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACA  >  NC_000913/1018‑1279
                                                                                                                                  |                                                                                                                                    
gacgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTAccccccatcg                                                                                                                            <  1:673665/140‑6 (MQ=255)
gacgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTAccccccatcg                                                                                                                            <  2:420617/140‑6 (MQ=255)
  cgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGcc                                                                                                                           <  2:640607/140‑1 (MQ=255)
  cgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGcc                                                                                                                           <  2:365889/140‑1 (MQ=255)
  cgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGcc                                                                                                                           <  2:445403/140‑1 (MQ=255)
  cgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGcc                                                                                                                           <  1:367722/140‑1 (MQ=255)
   gTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGccc                                                                                                                          <  1:45725/140‑1 (MQ=255)
    ttGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACTATTACCCCCCATCGCCCa                                                                                                                         <  1:550942/140‑1 (MQ=255)
       aCGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGtt                                                                                                                      <  2:451203/140‑1 (MQ=255)
          gggTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCa                                                                                                                   <  1:1036220/140‑1 (MQ=255)
           ggTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAg                                                                                                                  >  2:482397/1‑140 (MQ=255)
             tCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGAt                                                                                                                >  2:430324/1‑140 (MQ=255)
                     tGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCa                                                                                                                   >  2:601935/1‑129 (MQ=255)
                     tGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCa                                                                                                                   <  1:601935/129‑1 (MQ=255)
                     tGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGcc                                                                                                        <  1:285241/140‑1 (MQ=255)
                      gCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCt                                                                                                       >  2:536347/1‑140 (MQ=255)
                      gCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCt                                                                                                       >  2:653620/1‑140 (MQ=255)
                       cGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTg                                                                                                      >  2:981882/1‑140 (MQ=255)
                       cGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTg                                                                                                      >  1:1075634/1‑140 (MQ=255)
                           ccGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTcc                                                                                                                    <  1:102160/122‑1 (MQ=255)
                           ccGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTcc                                                                                                                    >  2:102160/1‑122 (MQ=255)
                                    gTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCg                                                                                         >  2:42980/1‑140 (MQ=255)
                                           aTGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATcc                                                                                  >  2:648466/1‑140 (MQ=255)
                                            tGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCt                                                                                 <  1:981882/140‑1 (MQ=255)
                                              cGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCa                                                                               <  1:425469/140‑1 (MQ=255)
                                              cGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCa                                                                               <  2:638760/140‑1 (MQ=255)
                                                   ttgttgAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCAc                                                                          <  2:597098/140‑1 (MQ=255)
                                                    tgttgAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCAcc                                                                         <  1:482397/140‑1 (MQ=255)
                                                        gAAGTNGATGTCCTACCAGGAAGCGATGGAGCTTTCCTNNNNCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGt                                                                     <  1:323217/140‑1 (MQ=255)
                                                        gAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGt                                                                     <  1:775502/140‑1 (MQ=255)
                                                           gTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACg                                                                  >  1:776098/1‑140 (MQ=255)
                                                             cGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCt                                                                >  2:317982/1‑140 (MQ=255)
                                                                tGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCAt                                                             >  1:961748/1‑140 (MQ=255)
                                                                     tACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCAtt                                                            <  1:651183/136‑1 (MQ=255)
                                                                     tACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCAtt                                                            >  2:651183/1‑136 (MQ=255)
                                                                      aCCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGc                                                       >  2:866828/1‑140 (MQ=255)
                                                                         aGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCGCCAGGTACGCTCATTGGTGCCAg                                                    <  2:61513/140‑1 (MQ=255)
                                                                            aaGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCg                                                 >  2:1001982/1‑140 (MQ=255)
                                                                               cGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtga                                              >  2:834598/1‑140 (MQ=255)
                                                                               cGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtga                                              >  2:430232/1‑140 (MQ=255)
                                                                                 aTGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtgatg                                            >  1:532351/1‑140 (MQ=255)
                                                                                 aTGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtgatg                                            >  1:831589/1‑140 (MQ=255)
                                                                                  tGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCTCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATACTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGtgatga                                           >  2:957846/1‑140 (MQ=255)
                                                                                  tGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtgatga                                           >  1:617953/1‑140 (MQ=255)
                                                                                  tGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGtgatga                                           <  1:508253/140‑1 (MQ=255)
                                                                                            cTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTa                                 <  2:776098/140‑1 (MQ=255)
                                                                                            cTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTAAGCTCATTGGTGCCAGCCGTGATGAAGACGAATTa                                 <  2:961748/140‑1 (MQ=255)
                                                                                             tACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTAc                                <  1:430324/140‑1 (MQ=255)
                                                                                               cTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCg                              >  2:1067319/1‑140 (MQ=255)
                                                                                                  cGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTAc                                <  2:130216/134‑1 (MQ=255)
                                                                                                  cGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTAc                                >  1:130216/1‑134 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGaa                                          <  1:147293/123‑1 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGaa                                          >  2:147293/1‑123 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCaa                         >  2:782540/1‑140 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCaa                         >  2:924031/1‑140 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCaa                         >  1:776898/1‑140 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGACGTGATGAAGACGAAt                                   >  2:701270/1‑130 (MQ=255)
                                                                                                    gcgcTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGACGTGATGAAGACGAAt                                   <  1:701270/130‑1 (MQ=255)
                                                                                                       cTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAggg                      <  2:1075634/140‑1 (MQ=255)
                                                                                                       cTAAAGTTCTTCACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAggg                      <  1:834598/140‑1 (MQ=255)
                                                                                                               cttcACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCa              >  1:128335/1‑140 (MQ=255)
                                                                                                               cttcACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCa              <  1:782540/140‑1 (MQ=255)
                                                                                                                ttcACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCaa             <  2:316434/140‑1 (MQ=255)
                                                                                                                 tcACCCCCGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAAt            <  2:532351/140‑1 (MQ=255)
                                                                                                                      cccGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGaa       <  1:648466/140‑1 (MQ=255)
                                                                                                                      cccGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGaa       <  1:434841/140‑1 (MQ=255)
                                                                                                                        cGCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATa     >  2:817020/1‑140 (MQ=255)
                                                                                                                         gCACCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATaa    <  1:653620/140‑1 (MQ=255)
                                                                                                                           aCCATTACCCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACa  <  1:866828/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACA  >  NC_000913/1018‑1279

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: