Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,533,038 G→C 5.8% R74P (CGT→CCT) ‡ pck → phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,533,0380GC 5.8% 432.7 / 12.8 120R74P (CGT→CCT) ‡pckphosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base G (54/59);  new base C (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 = 9.36e-01

Only 100 of 122 total aligned reads displayed.

AACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGT  >  NC_000913/3532893‑3533176
                                                                                                                                                 |                                                                                                                                          
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)
                            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)
                                           gCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggt                                                                                                               >  3:211737/1‑132 (MQ=255)
                                            ccTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGGCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                              <  1:448037/148‑1 (MQ=255)
                                               gACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    <  4:645226/139‑1 (MQ=255)
                                               gACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    >  3:645226/1‑139 (MQ=255)
                                                 cAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                   >  2:527214/1‑122 (MQ=255)
                                                 cAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTaa                                                                                         >  2:262101/1‑148 (MQ=255)
                                                   ggTTATGAGCGCGGGGTGTTAGCCAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaagg                                                                                                  <  2:264010/137‑1 (MQ=255)
                                                      tATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAg                                                                                        >  4:437337/1‑144 (MQ=255)
                                                         gAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                 <  3:1109126/148‑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)
                                                            cgcgGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATat                                                                                                                                   <  2:34510/95‑2 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGt                                                                                           >  1:283277/1‑125 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGaa                                                                                      <  1:595963/130‑1 (MQ=255)
                                                                      tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTCACGATACCACTCGCTATACTTTCTGGTGGGCAGAcaaa                                                                                                    <  2:638825/116‑1 (MQ=255)
                                                                        aCTAATCTGGGTGCCATTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCa                                                                                                                                 <  1:16748/85‑1 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    <  4:510423/113‑1 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    >  3:510423/1‑113 (MQ=255)
                                                                         cTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTcc                                                                 >  3:1030756/1‑148 (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:993951/1‑112 (MQ=255)
                                                                          tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    >  3:946081/1‑112 (MQ=255)
                                                                           aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtgg                                                                                                                <  4:176448/99‑1 (MQ=255)
                                                                             tCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaa                                                             <  1:262101/148‑1 (MQ=255)
                                                                                  gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                             >  4:29837/1‑127 (MQ=255)
                                                                                  gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                             <  3:29837/127‑1 (MQ=255)
                                                                                      cgttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgcag                                                    >  1:475776/1‑148 (MQ=255)
                                                                                         tgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                  >  2:822715/1‑131 (MQ=255)
                                                                                         tgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                  <  1:822715/131‑1 (MQ=255)
                                                                                             gtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCg                                                                                                                                            <  4:116880/53‑1 (MQ=255)
                                                                                             gtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCg                                                                                                                                            >  3:116880/1‑53 (MQ=255)
                                                                                                   aCCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGcc                                       <  1:593107/148‑1 (MQ=255)
                                                                                                   aCCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGcc                                       <  1:472756/148‑1 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGacaa                                                                                <  1:63556/104‑1 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGacaa                                                                                >  2:63556/1‑104 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                             <  3:434541/107‑1 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGAcaacaa                                                                             >  4:434541/1‑107 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                            >  4:1044459/1‑124 (MQ=255)
                                                                                                      gggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                            <  3:1044459/124‑1 (MQ=255)
                                                                                                        gATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    >  3:235359/1‑82 (MQ=255)
                                                                                                        gATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    <  4:235359/82‑1 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                              >  3:111635/1‑85 (MQ=255)
                                                                                                           cTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                       <  1:881290/124‑1 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATa                                                                                                                                    <  4:303499/45‑1 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATa                                                                                                                                    >  3:303499/1‑45 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                   <  1:643413/78‑1 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                   >  2:643413/1‑78 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                   >  4:147293/1‑78 (MQ=255)
                                                                                                             tCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                   <  3:147293/78‑1 (MQ=255)
                                                                                                                ccGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGTGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaaa                                                            >  4:661172/1‑114 (MQ=255)
                                                                                                                  ggTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGc                           >  1:601557/1‑145 (MQ=255)
                                                                                                                  ggTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGc                           <  2:601557/145‑1 (MQ=255)
                                                                                                                    tCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                    >  3:639719/1‑70 (MQ=255)
                                                                                                                       ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTCAAAGGCCTGGTGACCAGGCa                          >  1:510491/1‑141 (MQ=255)
                                                                                                                       ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTCAAAGGCCTGGTGACCAGGCa                          <  2:510491/141‑1 (MQ=255)
                                                                                                                        tCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                              >  4:860188/1‑72 (MQ=255)
                                                                                                                        tCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaa                                                                                              <  3:860188/72‑1 (MQ=255)
                                                                                                                         cACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGaa                                            <  1:791002/121‑1 (MQ=255)
                                                                                                                         cACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGaa                                            >  2:791002/1‑121 (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)
                                                                                                                              aaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCAt                                                 <  2:45243/111‑1 (MQ=255)
                                                                                                                                 gATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgc                                                      <  2:520166/103‑1 (MQ=255)
                                                                                                                                 gATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGgc                                                      >  1:520166/1‑103 (MQ=255)
                                                                                                                                      tggtatCGTCACGGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTc    >  3:421419/4‑148 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                 <  4:189432/70‑1 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGaca                                                                                 >  3:189432/1‑70 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                        <  3:581682/95‑1 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                        >  4:210221/1‑95 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                        >  4:581682/1‑95 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                        <  3:210221/95‑1 (MQ=255)
                                                                                                                                       tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCg   >  2:598553/1‑148 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAgg                            <  1:123898/122‑1 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAgg                            >  2:123898/1‑122 (MQ=255)
                                                                                                                                        atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGt  <  1:598553/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)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCa                              >  4:109902/1‑114 (MQ=255)
                                                                                                                                              tCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCa                              <  3:109902/114‑1 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTaa                                                                                         <  3:316701/54‑1 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                  <  1:765124/77‑1 (MQ=255)
                                                                                                                                               ccGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                  >  2:765124/1‑77 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
AACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGT  >  NC_000913/3532893‑3533176

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