Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,454,594 T→C 5.4% intergenic (+16/+17) yjfF → / ← fbp putative sugar ABC transporter permease/fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,454,5940TC 5.4% 405.2 / 9.5 129intergenic (+16/+17)yjfF/fbpputative sugar ABC transporter permease/fructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base T (57/65);  new base C (2/5);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 137 total aligned reads displayed.

AGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCT  >  NC_000913/4454447‑4454739
                                                                                                                                                   |                                                                                                                                                 
aGTTCCTGGTGGACGAATATCGCCATCGGTATGTTGTTGATTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCCGTGACAAGAGTCAACATTGCTCAGCAATAAAAACGCCTCTCCgtgt                                                                                                                                                   <  1:165250/148‑1 (MQ=255)
aGTTCCTGGTGGACGAAAATCGCCATCGTTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtgt                                                                                                                                                   <  1:312157/148‑1 (MQ=255)
aGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGGGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGtgtg                                                                                                                                                  <  3:20388/149‑1 (MQ=255)
aGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGAGTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCTAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtgt                                                                                                                                                   <  6:26047/148‑1 (MQ=255)
 gTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGtgtg                                                                                                                                                  <  5:362314/148‑1 (MQ=255)
 gTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTgg                                                                                                                                                 <  7:331592/149‑1 (MQ=255)
  ttCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                <  8:185553/149‑1 (MQ=255)
  ttCCTGGGGGACGAAAATCGCCAGCGGTAGTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                <  2:316901/149‑1 (MQ=255)
         tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                         <  4:292721/149‑2 (MQ=255)
           gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                <  8:301810/140‑1 (MQ=255)
           gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                >  7:301810/1‑140 (MQ=255)
                   tCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTAc                                                                                                                               <  1:334996/149‑1 (MQ=255)
                       cATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGt                                                                                                                           >  4:143923/1‑149 (MQ=255)
                         tCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGtgtg                                                                                                                                                  <  7:113360/124‑1 (MQ=255)
                         tCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGtgtg                                                                                                                                                  >  8:113360/1‑124 (MQ=255)
                           ggTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCgg                                                                                                                       <  3:139400/149‑1 (MQ=255)
                                tttGTTGTTTATTTTTTTAGCATTACAGCGTGGAGTATCGGTGCTGTGGGAGAATCGTCAGAGTACGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAAc                                                                                                                   <  2:98455/148‑1 (MQ=255)
                                 ttgttgTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACt                                                                                                                  <  4:166132/148‑1 (MQ=255)
                                   gttgttTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCa                                                                                                                <  7:68816/148‑1 (MQ=255)
                                      gttTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                <  6:220449/113‑1 (MQ=255)
                                      gttTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGga                                                                                                                                                >  5:220449/1‑113 (MQ=255)
                                          aTTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATa                                                                                                          >  1:257822/1‑147 (MQ=255)
                                          aTTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgc                                                                                                                                           <  4:210325/114‑1 (MQ=255)
                                           tttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                        <  5:198447/148‑1 (MQ=255)
                                           tttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAg                                                                                                       >  5:113817/1‑149 (MQ=255)
                                           tttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAg                                                                                                       <  3:143923/149‑1 (MQ=255)
                                                ataGAATTACAGCGTGTATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCGTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGtt                                                                                                   <  2:187440/148‑1 (MQ=255)
                                                   gCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGa                                                                                                <  3:347354/148‑1 (MQ=255)
                                                    cATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACa                                                                                              <  3:184128/149‑1 (MQ=255)
                                                          aGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTc                                                                                         >  5:199374/1‑148 (MQ=255)
                                                           gCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCaa                                                                                       <  6:113817/149‑1 (MQ=255)
                                                            cGTGGATTAACGGTGCGGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAAc                                                                                      <  2:199425/149‑1 (MQ=255)
                                                               ggATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTc                                                                                                                          <  6:183020/110‑1 (MQ=255)
                                                                gATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCata                                                                                  >  5:219695/1‑149 (MQ=255)
                                                                 aTTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCatat                                                                                 <  6:132888/149‑1 (MQ=255)
                                                                  ttAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATTAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGTAACTCTCGGATAAAGCGTTCGACATCGTCAACCTTATg                                                                                >  2:261288/1‑149 (MQ=255)
                                                                  ttAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                                                >  6:348912/1‑149 (MQ=255)
                                                                     aCGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCgg                                                                                                                       <  5:196913/107‑1 (MQ=255)
                                                                       ggTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGt                                                                           <  4:172017/149‑1 (MQ=255)
                                                                        gTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGtt                                                                          <  6:161116/149‑1 (MQ=255)
                                                                             gtgGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGa                                                                     >  7:72435/1‑149 (MQ=255)
                                                                                 gagaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACaaa                                                                 <  6:219695/149‑1 (MQ=255)
                                                                                    aaTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                              >  8:193338/1‑149 (MQ=255)
                                                                                        gTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGaa                                                          >  3:217406/1‑149 (MQ=255)
                                                                                         tCAGATTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGTAGAGGGGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAAc                                                         <  8:274140/149‑1 (MQ=255)
                                                                                         tCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAAc                                                         >  3:214309/1‑149 (MQ=255)
                                                                                         tCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAAc                                                         >  8:355871/1‑149 (MQ=255)
                                                                                         tCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAAc                                                         >  2:69512/1‑149 (MQ=255)
                                                                                            gagTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgc                                                      >  1:117859/1‑149 (MQ=255)
                                                                                              gTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgcgc                                                    <  1:347836/149‑1 (MQ=255)
                                                                                              gTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgcgc                                                    <  7:355871/149‑1 (MQ=255)
                                                                                               ttCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCt                                                   <  2:117859/149‑1 (MQ=255)
                                                                                                   ccAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTg                                               >  1:99567/1‑149 (MQ=255)
                                                                                                      gTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGcc                                                                       <  1:129068/122‑1 (MQ=255)
                                                                                                      gTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGcc                                                                       >  2:129068/1‑122 (MQ=255)
                                                                                                       tGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           <  1:69512/149‑1 (MQ=255)
                                                                                                       tGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           <  6:199374/149‑1 (MQ=255)
                                                                                                       tGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           <  1:189461/149‑1 (MQ=255)
                                                                                                       tGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           <  5:70579/149‑1 (MQ=255)
                                                                                                         aCAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACgg                                                       <  4:231579/135‑1 (MQ=255)
                                                                                                         aCAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGt                                         <  8:72435/149‑1 (MQ=255)
                                                                                                         aCAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGt                                         >  6:150566/1‑149 (MQ=255)
                                                                                                          cAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGtt                                        >  6:340620/1‑149 (MQ=255)
                                                                                                           aaGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACaaa                                                                 >  4:197395/1‑123 (MQ=255)
                                                                                                           aaGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACaaa                                                                 <  3:197395/123‑1 (MQ=255)
                                                                                                           aaGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGttt                                       >  2:385739/1‑149 (MQ=255)
                                                                                                           aaGAGTCAACATTGCGCAGCAATAAAAAAGCCTCTCCGTGTGGAGAGGAGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGttt                                       <  1:261288/149‑1 (MQ=255)
                                                                                                               gTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCgg                                   >  3:88890/1‑149 (MQ=255)
                                                                                                                tCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTc                                      >  3:219773/1‑145 (MQ=255)
                                                                                                                tCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTc                                      <  4:219773/145‑1 (MQ=255)
                                                                                                                  aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                                                <  1:351154/101‑1 (MQ=255)
                                                                                                                  aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                                                >  2:351154/1‑101 (MQ=255)
                                                                                                                  aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                >  2:346945/1‑149 (MQ=255)
                                                                                                                  aaCATTGCGCAGCAATAAAAACGCCGCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                <  7:75259/149‑1 (MQ=255)
                                                                                                                   aCATTGCGCAGCAATAAATACGCCTCTGCGTGTGTAGAGGCGCATGAGATTACGCGTCCGGGAACTCACGGTTAAAGCGTTCGACATGTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           <  4:332827/137‑1 (MQ=255)
                                                                                                                   aCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                                           >  3:332827/1‑137 (MQ=255)
                                                                                                                    cATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                <  5:264722/147‑1 (MQ=255)
                                                                                                                    cATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                >  6:264722/1‑147 (MQ=255)
                                                                                                                      ttGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGata                            <  5:150566/149‑1 (MQ=255)
                                                                                                                         cgcAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATc                          <  4:88890/148‑1 (MQ=255)
                                                                                                                          gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTgg                                                 <  6:174858/124‑1 (MQ=255)
                                                                                                                               aaTAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATa                   >  6:337204/1‑149 (MQ=255)
                                                                                                                               aaTAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATGACGCGTCCGGGACCTCCCGGAGAACGCGTTCGCCATCTTCACCCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGGGCAGGGTTTCCGGGATGATTTCCAGAATa                   >  1:83214/1‑149 (MQ=255)
                                                                                                                                   aaaaCGCCTCTCCGTGTTGAGAGGCTGAGGAGATTACGCGTCCGTGAACTCAGGGATAAAGCGTTCGACATGTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGc                <  1:81466/148‑1 (MQ=255)
                                                                                                                                    aaaCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTGCGGGATGATATCCAGAATACGctc              >  5:100757/1‑149 (MQ=255)
                                                                                                                                     tgcGCCTCTCCACACGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTc                                                                                                  <  6:77402/62‑1 (MQ=25)
                                                                                                                                     tgcGCCTCTCCACACGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTc                                                                                                  >  5:77402/3‑64 (MQ=25)
                                                                                                                                       cGCCTCTCCACACGGAGAGGCGTAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCttt           <  4:217891/149‑1 (MQ=255)
                                                                                                                                          ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTc                                                                             <  8:7844/80‑1 (MQ=255)
                                                                                                                                          ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTc                                                                             >  6:291929/1‑80 (MQ=255)
                                                                                                                                           tctcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                                                >  3:287079/1‑76 (MQ=255)
                                                                                                                                           tctcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCa       >  5:40133/1‑149 (MQ=255)
                                                                                                                                            ctcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATATTCAACCATATGGTCGTTTCCGACAAAGAATGAACGGCGCTCGTGCAGGGTTTCCGCGAAGATATCCAGAAAACGCTCTTTGCCa       >  4:318249/1‑148 (MQ=255)
                                                                                                                                             tcCGTGTGGAGTGGCGCAGGAGATTACGCGTCCGGGATCTCACGGATTAAGCGTTCGACATCTTCAACCATTTGGTCGTTTCCGACTAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCttt           <  3:71804/143‑1 (MQ=255)
                                                                                                                                             tcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCa       <  2:99567/147‑1 (MQ=255)
                                                                                                                                             tcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtc     >  8:213818/1‑149 (MQ=255)
                                                                                                                                             tcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtc     >  1:31011/1‑149 (MQ=255)
                                                                                                                                             tcCGTGTGGAGAGGCGCAGGAGATTACGCGTACGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGAAATCCAGAATACGcagttt           >  4:71804/1‑138 (MQ=255)
                                                                                                                                               cGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtc     <  2:370666/147‑1 (MQ=255)
                                                                                                                                                gtgtGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtcgct  >  8:108315/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCT  >  NC_000913/4454447‑4454739

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