Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,568,986 G→A 5.3% R115C (CGC→TGC)  glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,568,9860GA 5.3% 243.3 / 2.2 75R115C (CGC→TGC) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/41);  new base A (2/2);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTGGGTGATCACGCCCATACGACGGA  >  NC_000913/3568841‑3569132
                                                                                                                                                 |                                                                                                                                                  
cAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGAt                                                                                                                                                  <  2:133065/148‑1 (MQ=255)
 aaCGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATAc                                                                                                                                                >  4:189209/1‑149 (MQ=255)
       cgcCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGttt                                                                                                                                          >  4:340198/1‑149 (MQ=255)
        gcCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGtttt                                                                                                                                         <  5:227776/149‑1 (MQ=255)
          cTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTcc                                                                                                                                       >  1:204812/1‑149 (MQ=255)
             tttCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTccc                                                                                                                                      <  3:6467/147‑1 (MQ=255)
                      gATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATtctc                                                                                                                           >  6:103863/1‑149 (MQ=255)
                      gATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATtctc                                                                                                                           >  7:380186/1‑149 (MQ=255)
                       aTCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATtctc                                                                                                                           <  7:371979/148‑1 (MQ=255)
                       aTCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATtctc                                                                                                                           <  7:377088/148‑1 (MQ=255)
                           aTAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGAGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTctgctg                                                                                                                      >  5:73338/1‑149 (MQ=255)
                              aGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCtgctgtgc                                                                                                                   <  8:65969/149‑1 (MQ=255)
                               gCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCGTCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCtgctgtgct                                                                                                                  >  3:204027/1‑149 (MQ=255)
                                     gcgAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgcag                                                                                                             >  1:98388/1‑148 (MQ=255)
                                     gcgAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgcag                                                                                                             >  4:38723/1‑148 (MQ=255)
                                             tcctgcttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGAc                                                                                                    <  6:46160/146‑1 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTgg                                                                                                                <  2:73820/131‑1 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgca                                                                                                              >  4:137301/1‑133 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgca                                                                                                              <  3:137301/133‑1 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTc                                                                                              >  4:350201/1‑149 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTc                                                                                              >  4:44593/1‑149 (MQ=255)
                                                   ttgTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTc                                                                                              >  5:214735/1‑149 (MQ=255)
                                                      tAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGtt                                                                                           <  3:340198/149‑1 (MQ=255)
                                                      tAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGt                                                                                            <  3:350201/148‑1 (MQ=255)
                                                        gATATGGTCGCCCGCCAGGATCACCACGTATTCCTCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCa                                                                                         >  2:341375/1‑149 (MQ=255)
                                                        gATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGa                                                                                                         >  6:214073/1‑133 (MQ=255)
                                                        gATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGa                                                                                                         <  5:214073/133‑1 (MQ=255)
                                                        gATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCAATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCa                                                                                         >  7:84603/1‑149 (MQ=255)
                                                           atGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCAATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCAttt                                                                                      >  1:85464/1‑149 (MQ=255)
                                                            tGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATttc                                                                                     <  6:73338/149‑1 (MQ=255)
                                                            tGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGATTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACACACTCGTTCATttc                                                                                     >  1:154902/1‑149 (MQ=255)
                                                              gTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATttctt                                                                                   >  8:100395/1‑149 (MQ=255)
                                                              gTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATttct                                                                                    <  1:353304/148‑1 (MQ=255)
                                                                cgcccgccAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCa                                                                                 <  2:98388/149‑1 (MQ=255)
                                                                 gcccgccAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCAt                                                                                >  2:338667/1‑149 (MQ=255)
                                                                   ccgccAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTg                                                                              <  3:189209/149‑1 (MQ=255)
                                                                     gccAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTg                                                                              >  2:178602/1‑147 (MQ=255)
                                                                      ccAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTgaag                                                                           <  2:99338/149‑1 (MQ=255)
                                                                        aGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTgaagaa                                                                         <  5:103863/149‑1 (MQ=255)
                                                                         ggATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAAt                                                                        <  5:35902/149‑1 (MQ=255)
                                                                         ggATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAAt                                                                        <  3:38723/149‑1 (MQ=255)
                                                                          gATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAtt                                                                                                                              >  5:69918/1‑94 (MQ=255)
                                                                          gATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAtt                                                                                                                              <  6:69918/94‑1 (MQ=255)
                                                                           aTCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTACTTCATTGAAGAATGa                                                                      >  4:282295/1‑149 (MQ=255)
                                                                               ccacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGt                                                                                            >  7:196330/1‑123 (MQ=255)
                                                                               ccacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGt                                                                                            <  8:196330/123‑1 (MQ=255)
                                                                               ccacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAg                                                                  >  8:249463/1‑149 (MQ=255)
                                                                               ccacGTATTCCGCTTAATAACGGCGCATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGGGCCGCGATACCAGTTTTCCCCTTTCATTCTCAGCTGTGCTGGCAGCACATCCACAAACACGTTCATTTCTTCATTTAAGAATCACCAg                                                                  >  8:299730/1‑149 (MQ=255)
                                                                                cacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGc                                                                 >  8:39677/1‑149 (MQ=255)
                                                                                cacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAGTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCAtt                                                                               >  2:243537/1‑135 (MQ=255)
                                                                                cacGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCACTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCAtt                                                                               <  1:243537/135‑1 (MQ=255)
                                                                                         cGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGaa                                                        <  1:178602/149‑1 (MQ=255)
                                                                                         cGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGaa                                                        <  4:26720/149‑1 (MQ=255)
                                                                                          gCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAAt                                                       <  2:146066/149‑1 (MQ=255)
                                                                                            tttATAACGTCGGATAATGTCGAGGTTTTGGGTGACCGCATCTCCCGTGCCGCGATACCAGTTAACCCCTTTCAATCTCTGCTGTGCTGGCAGCAGATAGACAAACTCGTTCATTTCTTCATTGAACAATGACCAGCCGCGCTGAAtgt                                                     >  5:358063/1‑149 (MQ=255)
                                                                                                  aCGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTCCTTCATTGAAGAATGACCAGCCGCGCTGAATGtgctgc                                                <  3:282295/148‑1 (MQ=255)
                                                                                                   cggcggATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCAc                                              <  8:380186/149‑1 (MQ=255)
                                                                                                     gcggATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCa                                            <  7:249463/149‑1 (MQ=255)
                                                                                                     gcggATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCa                                            <  7:298136/149‑1 (MQ=255)
                                                                                                       ggATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCaga                                          <  6:214735/149‑1 (MQ=255)
                                                                                                            aTGTGGAGGTTTTGGGTTACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCTTTCATTTCTTCATTGAAGAATGACCAGCCGCTCTGAATGTGCTGCACAAGAGtgtg                                     <  6:358063/149‑1 (MQ=255)
                                                                                                            aTGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCAATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGtgtg                                     <  8:84603/149‑1 (MQ=255)
                                                                                                             tGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTgg                                    <  7:39677/149‑1 (MQ=255)
                                                                                                             tGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCAATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTAAATGTGCTGCACCAGAGTGTgg                                    <  2:85464/149‑1 (MQ=255)
                                                                                                                   ggTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTgg                              <  1:338667/149‑1 (MQ=255)
                                                                                                                          ggTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCACCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTgg                        <  4:204027/148‑1 (MQ=255)
                                                                                                                           gTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTg                               >  5:214062/1‑140 (MQ=255)
                                                                                                                           gTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTg                               <  6:214062/140‑1 (MQ=255)
                                                                                                                             gACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCt                                                           <  2:332938/110‑1 (MQ=255)
                                                                                                                             gACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCt                                                           >  1:332938/1‑110 (MQ=255)
                                                                                                                               ccGCATCTGCGGTGCCGCGATTCCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGGAGATCGACTAACTCGTTCATTTCTTCATTGAATAATGACGAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGAACTGGTACTGGGTTATc                  <  8:132924/149‑1 (MQ=255)
                                                                                                                               ccGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTGGGTGATc                  <  3:44593/149‑1 (MQ=255)
                                                                                                                                   aTCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTGGGTGATCACGc              <  7:299730/149‑1 (MQ=255)
                                                                                                                                       gCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCaga                                          <  8:170790/117‑1 (MQ=255)
                                                                                                                                       gCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCaga                                          >  7:170790/1‑117 (MQ=255)
                                                                                                                                               gcgATACCCGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTGGGTGATCACGCCCATACGACGGa  <  2:204812/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCATTGAAGAATGACCAGCCGCGCTGAATGTGCTGCACCAGAGTGTGGGACTGGTACTGGGTGATCACGCCCATACGACGGA  >  NC_000913/3568841‑3569132

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: