Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,046 G→T 100% P164T (CCG→ACG)  guaA ← GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,0460GT97.8% 301.2 / ‑0.9 89P164T (CCG→ACG) guaAGMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (45/42);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-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.

CGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGTTTA  >  NC_000913/2631903‑2632189
                                                                                                                                               |                                                                                                                                               
cGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCa                                                                                                                                              <  8:52115/147‑1 (MQ=255)
 gCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCCAACGTGCAGc                                                                                                                                            >  7:89812/1‑148 (MQ=255)
 gCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAg                                                                                                                                             <  7:117603/147‑1 (MQ=255)
         ccACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTg                                                                                                                                    >  8:314175/1‑148 (MQ=255)
         ccACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGACATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTg                                                                                                                                    >  8:308794/1‑148 (MQ=255)
          cacaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGc                                                                                                                                   <  7:328712/148‑1 (MQ=255)
          cacaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGc                                                                                                                                   <  6:223597/148‑1 (MQ=255)
           acaGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  <  7:46951/148‑1 (MQ=255)
           acaGGGCTTCACACTAGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  <  8:136907/148‑1 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  <  3:178621/147‑1 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  >  4:178621/1‑147 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTg                                                                                                                                 >  5:124072/1‑148 (MQ=255)
            caGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTg                                                                                                                                 >  8:364807/1‑148 (MQ=255)
               ggCTTCAGACTGGCAGATATCTCGCACAAAACGCTCGAGCATGCGCATACCCTGGCGGGTATGAGTCAGTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCt                                                                                                                              <  4:90608/148‑1 (MQ=255)
               ggCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCt                                                                                                                              <  8:83750/148‑1 (MQ=255)
                    caccagggCAGATAGCACGCACACACCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCAATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACg                                                                                                                           <  2:264918/140‑1 (MQ=255)
                    cacaCTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGt                                                                                                                         >  4:192480/1‑148 (MQ=255)
                    cacaCTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGt                                                                                                                         >  1:183307/1‑148 (MQ=255)
                    cacaCTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGt                                                                                                                         >  5:64657/1‑148 (MQ=255)
                     acacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtg                                                                                                                        <  2:183307/148‑1 (MQ=255)
                      cacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  <  4:102762/137‑1 (MQ=255)
                      cacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                                  >  3:102762/1‑137 (MQ=255)
                      cacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtga                                                                                                                       >  8:306098/1‑148 (MQ=255)
                        cTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtgat                                                                                                                      >  8:10179/1‑147 (MQ=255)
                        cTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtgat                                                                                                                      >  2:133299/1‑147 (MQ=255)
                          ggCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtgatg                                                                                                                     >  5:54339/1‑146 (MQ=255)
                           gCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtgatga                                                                                                                    <  4:23785/146‑1 (MQ=255)
                              gATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcg                                                                                                               <  6:124072/148‑1 (MQ=255)
                               atatCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcgg                                                                                                              <  2:3227/148‑1 (MQ=255)
                                   cACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                                                                                                           >  3:4863/1‑147 (MQ=255)
                                   cACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGGGATGAAGTcggacg                                                                                                           <  8:89812/147‑1 (MQ=255)
                                       cacaAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATa                                                                                                      <  7:306098/148‑1 (MQ=255)
                                         caAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc                                                                                                    <  3:192480/148‑1 (MQ=255)
                                         caAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc                                                                                                    >  7:107852/1‑148 (MQ=255)
                                          aaaaCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCg                                                                                                   <  1:90159/148‑1 (MQ=255)
                                            aaCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCgg                                                                                                  >  7:296304/1‑147 (MQ=255)
                                             aCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTa                                                                                                <  1:212/148‑1 (MQ=255)
                                                     gCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATc                                                                                        >  5:322669/1‑148 (MQ=255)
                                                       aTGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCga                                                                                      >  1:204236/1‑147 (MQ=255)
                                                            cATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGGTACGGTGATGAAGTCGGACGGATTAGCGGTAACTTTATCGCCGTgg                                                                                 >  5:55953/1‑148 (MQ=255)
                                                              tACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCt                                                                               <  4:210569/148‑1 (MQ=255)
                                                               aCCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCCCCGTGGCt                                                                               >  4:96624/1‑147 (MQ=255)
                                                                  cTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATc                                                                           >  4:205740/1‑148 (MQ=255)
                                                                     gCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAg                                                                        <  3:114165/148‑1 (MQ=255)
                                                                     gCGGGAATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAg                                                                        <  3:293450/148‑1 (MQ=255)
                                                                      cGGGTATGAGTATCTTCCGTGTGGATCTGTACGCCATAGAAGCTTTTTTCTTCGTTAGCCATAATGGCTAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGa                                                                       <  6:55953/148‑1 (MQ=255)
                                                                       gggTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGAc                                                                      <  8:107852/148‑1 (MQ=255)
                                                                            tGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCg                                                                  <  6:64657/147‑1 (MQ=255)
                                                                            tGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGa                                                                 <  3:66568/148‑1 (MQ=255)
                                                                             gAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGag                                                                >  1:225656/1‑148 (MQ=255)
                                                                                  aCTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCg                                                           <  4:4863/148‑1 (MQ=255)
                                                                                   cTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGGGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCgg                                                          <  6:7571/148‑1 (MQ=255)
                                                                                         ggTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACAt                                                                    >  4:30828/1‑132 (MQ=255)
                                                                                         ggTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACAt                                                                    <  3:30828/132‑1 (MQ=255)
                                                                                               aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGc                                              >  7:250355/1‑148 (MQ=255)
                                                                                               aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAAAGGCAAACGTGAAGCTTTCGGTGCAGGCTACGCTGATGAAGTACGACGGAATAGAGGTAACTTTATCGCCTTGGCACAACCAGACAACGACCAGCGCATAACCGACTGc                                              >  8:258419/1‑148 (MQ=255)
                                                                                                  gTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGt                                           >  1:882/1‑148 (MQ=255)
                                                                                                     cGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAg                                        <  8:235214/148‑1 (MQ=255)
                                                                                                         aTAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCgg                                            <  5:24325/140‑1 (MQ=255)
                                                                                                         aTAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCgg                                            >  6:24325/1‑140 (MQ=255)
                                                                                                         aTAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCa                                    >  8:123865/1‑148 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTTATGGCTACTGTGATGAAGTCGGACGGAAAAGCGGTAACTTTATCGCCGTGGCTCATCCCGACATCGAGCAGCGGTTTACCGTATGCGGTCAGCGCATc                                  >  7:264791/1‑148 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCgg                                                          >  7:342416/1‑124 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCgg                                                          <  8:342416/124‑1 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATc                                  >  1:136602/1‑148 (MQ=255)
                                                                                                           aGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATc                                  >  1:9223/1‑148 (MQ=255)
                                                                                                            gacgCGTGTGTCTTCGTTCGCCCGACGGGCACACGTGCAGCGTGCGGTGCGGGCTACGGGGAGGACGGCGGCCGGTATCGCGGTACCGTGCTCGCCGTGGCGCATCCCGACCTCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCctct                                 <  3:5252/145‑1 (MQ=255)
                                                                                                              aGCGTTTTTCTTCGTTAGCCCTAATGGCCAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGTAGAGCGGTAACTGTCTCGCCGTGGCGCATCCCGCCATCGAGCCGCGGGTTACCGGCGGCGGTCCGCGCATCtg                                >  4:5252/1‑146 (MQ=255)
                                                                                                              aGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTc                               <  2:88974/148‑1 (MQ=255)
                                                                                                              aGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTc                               >  2:179216/1‑148 (MQ=255)
                                                                                                               gCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCg                              <  3:35809/148‑1 (MQ=255)
                                                                                                               gCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCg                              <  7:308794/148‑1 (MQ=255)
                                                                                                               gCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCg                              <  7:314175/148‑1 (MQ=255)
                                                                                                               gCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCg                              <  2:373560/148‑1 (MQ=255)
                                                                                                                cGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATGGCGGAACGTTTAACGCGGGGGCGCAGCAGGCCAGCGCGCGGCGTTTTACCGTGTGCGGACGCCGCAGGTTCGa                             >  6:81076/1‑148 (MQ=255)
                                                                                                                      tcttcGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACcg                        <  1:179216/147‑1 (MQ=255)
                                                                                                                        ttcGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGa                     <  2:225656/148‑1 (MQ=255)
                                                                                                                         tcGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGaa                    >  3:145636/1‑148 (MQ=255)
                                                                                                                                 cATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTgcg            >  2:214438/1‑148 (MQ=255)
                                                                                                                                   tAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGAAGCGGTTTACAGTCAGCGGTCAGCGCATCTTCGATACCGCGAACCAGtcagct          >  5:311135/1‑143 (MQ=255)
                                                                                                                                    aaTGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTg         >  6:321302/1‑148 (MQ=255)
                                                                                                                                    aaTGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTg         >  5:249072/1‑148 (MQ=255)
                                                                                                                                    aaTGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGAGCGCTg         >  7:237671/1‑148 (MQ=255)
                                                                                                                                    aaTGGCAAACGTGCAGCTTTCGGTGCTGGATACGGTGATGAACTCGGACGTAATAGCGGTAACTTTATCGCCGTGGCGCATCCAGACATCGAGAAGCGGATAACCGTCTTCGGTCAGCGCAACTTCGAAACCTCGAACCAGTGCCCtc         >  3:25071/1‑147 (MQ=255)
                                                                                                                                      tGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTc       >  3:220277/1‑148 (MQ=255)
                                                                                                                                         cAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGtt    <  7:123865/148‑1 (MQ=255)
                                                                                                                                          aaaCGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGttt   <  7:364807/148‑1 (MQ=255)
                                                                                                                                          aaaCGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGttt   <  7:10179/148‑1 (MQ=255)
                                                                                                                                           aaCGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGCCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGTTTa  >  5:33850/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CGCTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTCGTTTA  >  NC_000913/2631903‑2632189

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