Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,294 A→T 9.2% L169Q (CTG→CAG)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2940AT 9.2% 236.8 / 8.9 89L169Q (CTG→CAG) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base T (4/5);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

GAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAG  >  NC_000913/4307146‑4307442
                                                                                                                                                    |                                                                                                                                                    
gAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                                                                                                     <  1:1645200/150‑1 (MQ=255)
 aaaCGGGCGTCTTCCGCAAGGACCAAAAGATCTCGCAATGTGTAATATCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGGTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                                                                                                    <  1:264530/150‑1 (MQ=255)
  aaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                   <  1:684657/150‑1 (MQ=255)
   aGGGGCGTTTTCCGCATGTACGAAAAGATCGCGCTATGGGTAATTTCGGGGCTGTTGTTCGTAGCAGCGTCTTCGCGCCATTCCAGAGCATTTTGTTTCCAGGCACCCAGGATTGCCACACGCGCTGTGTTGGGTCATATCTCCCTgggg                                                                                                                                                  <  2:241164/150‑1 (MQ=255)
   aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                  >  2:960140/1‑150 (MQ=255)
    ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                 <  2:256071/150‑1 (MQ=255)
         gTTTTCCGCATTGACCAAAAGATCGCGCAAGTGGTAATTTCGGGTCTGTTGTGCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTCGTGTCCAGGCACCCCGGCTCGCCACACGCGCAGTGTTGGGTCATAGCTCCCAGGGGGATATg                                                                                                                                            <  2:1247148/150‑1 (MQ=255)
         gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                            <  2:1285688/150‑1 (MQ=255)
           tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                          <  1:780618/150‑1 (MQ=255)
            ttCCGCATGGACGAAAAGAGCGCGCAATGGGTAATTTCGGGGCGGTTGTTCGGACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGAGTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                         <  1:581430/150‑1 (MQ=255)
             tCCGCATGGACGAAAAGATCGCGCAATGGGTAATATCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                        <  2:123447/150‑1 (MQ=255)
             gccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGGCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                        <  1:1340339/149‑1 (MQ=255)
              ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCa                                                                                                                                       <  2:1290108/150‑1 (MQ=255)
              ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCa                                                                                                                                       <  2:668542/150‑1 (MQ=255)
               cGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTCGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGAGATGCCCCAg                                                                                                                                      >  2:66265/1‑150 (MQ=255)
                  atgtgcGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTc                                                                                                                                   <  2:1444327/145‑1 (MQ=255)
                    ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                 <  1:363412/150‑1 (MQ=255)
                    ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                 <  2:1447402/150‑1 (MQ=255)
                    ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGATCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTTTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                 >  2:1065715/1‑150 (MQ=255)
                     gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                >  1:534054/1‑150 (MQ=255)
                     gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                          >  2:343622/1‑140 (MQ=255)
                     gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                          <  1:343622/140‑1 (MQ=255)
                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              >  2:1592220/1‑150 (MQ=255)
                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              >  2:1057844/1‑150 (MQ=255)
                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACACAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              >  2:865038/1‑150 (MQ=255)
                       cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              <  1:1358014/150‑1 (MQ=255)
                        gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                             <  2:87443/150‑1 (MQ=255)
                        gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATATACCAGGGGGATATGAACCAGATAGACTGc                                                                                                                             >  1:1669488/1‑150 (MQ=255)
                        caaaaGACCGCGCAATTGGTAATTTCGGGGCGGGTCTGCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                             <  2:111742/149‑1 (MQ=255)
                         aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                            >  1:1294959/1‑150 (MQ=255)
                            aGATCGCGAAATGGGTAATTTCGGTGCTGTTGTTCGTACAAGCGTCTTAGCGACATTCCAGAGCTATTGCTTGCCCGGCCCCCAGGATTGCCACACGCGCAGTGTTGGGTCGTATCTCCCAGGGGGATATGACCCCGTTCGCCTGCaacc                                                                                                                         >  2:954231/1‑149 (MQ=255)
                             gATCGCGCAATGGGTAATTTCGGTGCTGTTGTTCGTACCAGCTGCTTATCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                                                                        <  2:1059745/150‑1 (MQ=255)
                                 gcgcAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCCGTTCGCCTGCCACAccttt                                                                                                                    >  2:76192/1‑147 (MQ=255)
                                    cAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtga                                                                                                                 >  2:960334/1‑149 (MQ=255)
                                         ggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                                                                                                            >  1:532862/1‑150 (MQ=255)
                                               ttCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                      >  1:234440/1‑150 (MQ=255)
                                                  ggggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGGAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcg                                                                                                   >  2:113200/1‑150 (MQ=255)
                                                   gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                                  <  2:234440/150‑1 (MQ=255)
                                                     gCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGTTATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCCCAc                                                                                                >  1:808045/1‑150 (MQ=255)
                                                      cTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGGTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAcc                                                                                               <  1:244140/150‑1 (MQ=255)
                                                       tgttgtTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGACTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                              >  2:1073250/1‑150 (MQ=255)
                                                          tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                           <  1:1399646/150‑1 (MQ=255)
                                                            ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                         <  2:1240598/150‑1 (MQ=255)
                                                                  ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                   >  1:1382880/1‑150 (MQ=255)
                                                                  ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCACAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                   >  1:1376225/1‑150 (MQ=255)
                                                                         cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGCTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCAc                                                                            >  1:65813/1‑150 (MQ=255)
                                                                         cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                            <  1:1592220/150‑1 (MQ=255)
                                                                           tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                          <  1:1103028/150‑1 (MQ=255)
                                                                             gcgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCa                                                                        >  1:166185/1‑150 (MQ=255)
                                                                                ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                     >  2:313777/1‑150 (MQ=255)
                                                                                      cAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCCCCGTTCATCCACACTGCGACCCCCATCCCCGTAc                                                               >  1:700510/1‑150 (MQ=255)
                                                                                           cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                          <  1:960140/150‑1 (MQ=255)
                                                                                           cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGGACCgaga                                                          >  1:292120/1‑150 (MQ=255)
                                                                                                 gTTTCCAGGCACCCAGGATTTCCACACGCGCAGTGTTGGGTCATATCTCCCAGTGTGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCTAACCCCATCCCCGTACCGAGATAGgcc                                                    <  1:571175/150‑1 (MQ=255)
                                                                                                  tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                   <  2:1382880/150‑1 (MQ=255)
                                                                                                  tttCCAGGAACCCAGGATTTCCACACGCGCAGTGTTGGGTCATATCTCCAAGTGGGATATGACCAAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                   <  2:1376225/150‑1 (MQ=255)
                                                                                                      cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                               >  2:606762/1‑150 (MQ=255)
                                                                                                        ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                             >  1:729669/1‑150 (MQ=255)
                                                                                                         gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAAc                                            <  1:1561004/150‑1 (MQ=255)
                                                                                                           aCCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCAGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCa                                          >  1:852083/1‑150 (MQ=255)
                                                                                                             ccAGGATTGCCACACGCGCAGTTTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                        <  1:1065715/150‑1 (MQ=255)
                                                                                                              cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCAACGGCGCACCGTTCATCCACACTGCGACCCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAgtt                                       >  2:93124/1‑150 (MQ=255)
                                                                                                                 gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                    >  1:79287/1‑150 (MQ=255)
                                                                                                                  aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGttgttg                                   <  1:1487252/150‑1 (MQ=255)
                                                                                                                     gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACCCTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGGTGTTGCgg                                >  2:456377/1‑149 (MQ=255)
                                                                                                                     gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                >  2:427481/1‑150 (MQ=255)
                                                                                                                     gCCACACGCGCAGTGTTGGGTAATATCTCCCAGGGGGATATTACCCAGTTAGACTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTTAGAACCACAACCCCGTTGCGAGAAAGGCCTCCAGAAACAGTTGTTGCga                                >  2:767134/1‑149 (MQ=255)
                                                                                                                      ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCCCCAGAACCAGTTGTTGCGTa                               <  1:825604/150‑1 (MQ=255)
                                                                                                                      ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTACGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTa                               <  1:1125206/150‑1 (MQ=255)
                                                                                                                       cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                              >  2:1647173/1‑150 (MQ=255)
                                                                                                                           cgcgCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCCGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCACCCCCTCCCCGGACCGAGATTGGCCGCCAGCACCAGTTGGTGCGTAAGGCg                          >  1:1347491/1‑150 (MQ=255)
                                                                                                                           cgcgCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg                          <  1:867863/150‑1 (MQ=255)
                                                                                                                           cgcgCAGTGTTGGGGCATATCTCCCAGGGGTATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg                          <  1:1063490/150‑1 (MQ=255)
                                                                                                                                  tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCTTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct                   <  2:1349139/150‑1 (MQ=255)
                                                                                                                                   gtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACTTTTTTTTTTAATGTTAAGGGGTCAACGGTGATCTACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                  <  2:3019/150‑1 (MQ=255)
                                                                                                                                    ttGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAGCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                 <  1:93124/150‑1 (MQ=255)
                                                                                                                                     tGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact                <  2:166185/150‑1 (MQ=255)
                                                                                                                                      gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta               >  1:641647/1‑150 (MQ=255)
                                                                                                                                      gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCACATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTgttcta               >  1:221336/1‑150 (MQ=255)
                                                                                                                                       ggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCTCCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTctactcc              >  1:797507/1‑148 (MQ=255)
                                                                                                                                        gTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACGAGTTGTTGCGTAAGGCGGGTTTCTACTACg             >  1:895298/1‑150 (MQ=255)
                                                                                                                                         tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGt            >  1:767192/1‑150 (MQ=255)
                                                                                                                                            tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCCCGGCGCACCGTTCATCCACCCTGCGCACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGGTTTCTACTCCGTccc         >  1:704767/1‑150 (MQ=255)
                                                                                                                                            tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc         <  1:570408/150‑1 (MQ=255)
                                                                                                                                            tatCTCCCAGGGGGATATGACACAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGAGAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTAGTACGTccc         >  2:1404786/1‑150 (MQ=255)
                                                                                                                                               cccccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAACCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTATTCTACAACGTCCCAgg      >  2:1478748/3‑150 (MQ=255)
                                                                                                                                                tcCCAGGGGGATTTGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCACCGTTCATCCACCCTGCGACCCCCATCCCCGTACCGAGATAGGCCGCCGGAACCCGGTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAggg     >  2:1007411/1‑149 (MQ=255)
                                                                                                                                                  ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGGTTTCTACTACGTCCCAGgaga   >  1:364737/1‑150 (MQ=255)
                                                                                                                                                   cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag  <  2:532862/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAG  >  NC_000913/4307146‑4307442

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