Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,533,039 T→G 5.8% R74R (CGT→CGG) ‡ pck → phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,533,0390TG 5.8% 439.1 / 15.4 120R74R (CGT→CGG) ‡pckphosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base T (54/59);  new base G (4/3);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 123 total aligned reads displayed.

AACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGTTGTCGACGCT  >  NC_000913/3532893‑3533186
                                                                                                                                                  |                                                                                                                                                   
aaCCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTg                                                                                                                                                    <  3:770943/148‑1 (MQ=255)
   ccAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACg                                                                                                                                                 >  4:108096/1‑148 (MQ=255)
       gCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATat                                                                                                                                             <  2:893151/148‑2 (MQ=255)
                    gTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATat                                                                                                                                             >  2:1009359/1‑134 (MQ=255)
                      aTCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACttt                                                                                                                              <  2:736749/148‑1 (MQ=255)
                            aaGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggtg                                                                                                                        <  3:108096/148‑1 (MQ=255)
                                  tCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTACTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          <  2:924462/140‑1 (MQ=255)
                                  tCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          >  1:924462/1‑140 (MQ=255)
                                  tCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATata                                                                                                                                            >  4:526169/1‑120 (MQ=255)
                                        cGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                            >  4:262234/1‑148 (MQ=255)
                                        cGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                            <  3:262234/148‑1 (MQ=255)
                                           gCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggt                                                                                                                         <  4:211737/132‑1 (MQ=255)
                                           gCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggt                                                                                                                         >  3:211737/1‑132 (MQ=255)
                                               gACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:645226/1‑139 (MQ=255)
                                                   ggTTATGAGCGCGGGGTGTTAGCCAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                            <  2:264010/137‑1 (MQ=255)
                                                      tATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAg                                                                                                  <  3:437337/144‑1 (MQ=255)
                                                         gAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                           <  1:896807/148‑1 (MQ=255)
                                                            cgcgGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGGTACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  2:235165/1‑126 (MQ=255)
                                                            cgcgGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGGTACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              <  1:235165/126‑1 (MQ=255)
                                                            cgcgGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATat                                                                                                                                             >  1:34510/1‑94 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  1:638825/1‑116 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGt                                                                                                     <  2:283277/125‑1 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGt                                                                                                     >  1:283277/1‑125 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGaa                                                                                                <  1:595963/130‑1 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGaa                                                                                                >  2:595963/1‑130 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTCACGATACCACTCGCTATACTTTCTGGTGGGCAGAcaaa                                                                                                              <  2:638825/116‑1 (MQ=255)
                                                                        aCTAATCTGGGTGCCATTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCa                                                                                                                                           <  1:16748/85‑1 (MQ=255)
                                                                        aCTAATCTGGGTGCCATTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCa                                                                                                                                           >  2:16748/1‑85 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:510423/1‑113 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              <  4:510423/113‑1 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTcc                                                                           >  2:593107/1‑148 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTcc                                                                           >  2:472756/1‑148 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATat                                                                                                                                             <  3:952561/82‑2 (MQ=255)
                                                                          tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAccaa                                                                                                              <  4:946081/112‑1 (MQ=255)
                                                                          tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:946081/1‑112 (MQ=255)
                                                                          tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:993951/1‑112 (MQ=255)
                                                                           aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          >  3:176448/1‑99 (MQ=255)
                                                                           aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          <  4:176448/99‑1 (MQ=255)
                                                                             tCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaa                                                                       <  1:262101/148‑1 (MQ=255)
                                                                                  gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       <  3:29837/127‑1 (MQ=255)
                                                                                  gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       >  4:29837/1‑127 (MQ=255)
                                                                                      cgttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgcag                                                              >  1:475776/1‑148 (MQ=255)
                                                                                         tgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                            >  2:822715/1‑131 (MQ=255)
                                                                                                 aTACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAgg                                                   <  4:1030756/148‑1 (MQ=255)
                                                                                                   aCCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGcc                                                 <  1:593107/148‑1 (MQ=255)
                                                                                                   aCCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGcc                                                 <  1:472756/148‑1 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       >  4:434541/1‑107 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       <  3:434541/107‑1 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                                      <  3:1044459/124‑1 (MQ=255)
                                                                                                        gATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:235359/1‑82 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                                        <  4:111635/85‑1 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                                        >  3:111635/1‑85 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                                 <  1:881290/124‑1 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                                 >  2:881290/1‑124 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATa                                                                                                                                              <  4:303499/45‑1 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATa                                                                                                                                              >  3:303499/1‑45 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                             >  2:643413/1‑78 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                             <  3:147293/78‑1 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                             <  1:643413/78‑1 (MQ=255)
                                                                                                                ccGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGTGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                                      >  4:661172/1‑114 (MQ=255)
                                                                                                                ccGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                                      <  3:661172/114‑1 (MQ=255)
                                                                                                                  ggTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGc                                     >  1:601557/1‑145 (MQ=255)
                                                                                                                  ggTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGc                                     <  2:601557/145‑1 (MQ=255)
                                                                                                                    tCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              >  3:639719/1‑70 (MQ=255)
                                                                                                                    tCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                              <  4:639719/70‑1 (MQ=255)
                                                                                                                       ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTCAAAGGCCTGGTGACCAGGCa                                    >  1:510491/1‑141 (MQ=255)
                                                                                                                       ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTCAAAGGCCTGGTGACCAGGCa                                    <  2:510491/141‑1 (MQ=255)
                                                                                                                        tCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                                        >  4:860188/1‑72 (MQ=255)
                                                                                                                         cACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGaa                                                      <  1:791002/121‑1 (MQ=255)
                                                                                                                           ccAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                            <  2:339192/65‑1 (MQ=255)
                                                                                                                           ccAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                            >  1:339192/1‑65 (MQ=255)
                                                                                                                             aaaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          <  1:218734/49‑1 (MQ=255)
                                                                                                                             aaaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                          >  2:218734/1‑49 (MQ=255)
                                                                                                                             aaaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAg                                   >  4:108540/1‑136 (MQ=255)
                                                                                                                             aaaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAg                                   <  3:108540/136‑1 (MQ=255)
                                                                                                                              aaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCAt                                                           >  1:45243/1‑111 (MQ=255)
                                                                                                                              aaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCAt                                                           <  2:45243/111‑1 (MQ=255)
                                                                                                                                 gATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgc                                                                >  1:520166/1‑103 (MQ=255)
                                                                                                                                 gATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgc                                                                <  2:520166/103‑1 (MQ=255)
                                                                                                                                      tggtatCGTCACGGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTc              >  3:421419/4‑148 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                           >  3:189432/1‑70 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                           <  4:189432/70‑1 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                                  >  4:581682/1‑95 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                                  >  4:210221/1‑95 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCg             >  2:598553/1‑148 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAgg                                      <  1:123898/122‑1 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGt            <  1:598553/148‑1 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGt            <  2:475776/148‑1 (MQ=255)
                                                                                                                                            cGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       <  1:286390/69‑1 (MQ=255)
                                                                                                                                            cGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                                       >  2:286390/1‑69 (MQ=255)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACaaa                                                                                      >  2:853347/1‑68 (MQ=255)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACaaa                                                                                      <  1:853347/68‑1 (MQ=255)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCa                                        <  3:109902/114‑1 (MQ=255)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCa                                        >  4:109902/1‑114 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTaa                                                                                                   >  4:316701/1‑54 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                            <  1:765124/77‑1 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                            >  2:765124/1‑77 (MQ=255)
                                                                                                                                                  tGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTgg                                              <  4:887/104‑1 (MQ=255)
                                                                                                                                                  tGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTgg                                              >  3:887/1‑104 (MQ=255)
                                                                                                                                                  tGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGTTGTCGACGCt  >  2:196682/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGTTGTCGACGCT  >  NC_000913/3532893‑3533186

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