Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,091 G→C 9.0% intergenic (+58/+170) glxK → / ← allE glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,0910GC 9.0% 170.7 / 9.5 65intergenic (+58/+170)glxK/allEglycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base C (4/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

GCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAA  >  NC_000913/542949‑543227
                                                                                                                                              |                                                                                                                                         
gCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCgg                                                                                                                                           <  2:993778/143‑1 (MQ=255)
 cTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGa                                                                                                                                          >  1:715935/1‑143 (MQ=255)
  tCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGAt                                                                                                                                         >  1:458968/1‑143 (MQ=255)
    gCCAGCTGTGAAGCCAATCTCTTCAACTGCGCGCGTAATTTTGCCTGCGTTATTAAAATAGGTCAGGGAATTAAAAACTACCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                                       <  2:1119834/143‑1 (MQ=255)
         cGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                  <  2:661827/143‑1 (MQ=255)
            tGAAACCAATCTCGTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAg                                                                                                                               <  1:756806/143‑1 (MQ=255)
               aaCAAATCTTGTCAACAGCGCGGGAAATATGGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAc                                                                                                                           <  1:203163/143‑1 (MQ=255)
                atctacctctCCAACAGCGCGTAAAATATTTCCTTCTCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGAGGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                                          <  2:604338/137‑1 (MQ=255)
                acgaaactctTGTACAGCGCGCGAAATTTTGTCCACGCCATTAAATGATGTGAGTATGTTTTAAACTGACCCTTTCCTTTAAAGCGGATGCGATGTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAct                                                                                                                          <  2:737986/137‑2 (MQ=255)
                   aaTCCCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTActgca                                                                                                                       <  1:543736/143‑5 (MQ=255)
                     tctcTTCAACATCTCGCTAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                                     <  1:296477/143‑1 (MQ=255)
                      ttctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcc                                                                                                                    <  2:1108149/142‑1 (MQ=255)
                        ctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGCTTTCTAGCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt                                                                                                                  >  1:518651/1‑143 (MQ=255)
                         ttCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAATACTAACCCTTACCTTTAAAGCGGATGCGTTTTATATCGCATAAGAGTGCAGTATTCATGCAGGATGCGGCATGAATACAatataatac                                                                                                                 >  1:589062/1‑138 (MQ=255)
                         ttCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttg                                                                                                                 >  2:62512/1‑142 (MQ=255)
                          tCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                <  1:216165/143‑1 (MQ=255)
                            aaCAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTg                                                                                                              >  1:982222/1‑143 (MQ=255)
                                  gcgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                        <  1:352209/143‑1 (MQ=255)
                                  gcgccAAATATTGCCTTCGCCATTAAAATAGGCCAGGGAATTAAAAACTAACGCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                        <  2:1061069/143‑1 (MQ=255)
                                    gctaaaTATTGCCTGCGCCATTAAAATAGGTAAGGGAAGTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgc                                                                                                      <  1:897425/140‑1 (MQ=255)
                                      gAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCCTGAGTACCTTTTCCTTCTTGAAAATCgcgc                                                                                                    >  1:234818/1‑143 (MQ=255)
                                      gAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCgcgc                                                                                                    >  1:1046609/1‑143 (MQ=255)
                                          tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCacct                                                                                                >  1:421097/1‑140 (MQ=255)
                                                   cgcCATTCAAATAGGTCAGGGAATTAAAAACTAAACCTTCCATTTTACGCGGATGAGATTTATATAGAATATGTGTGCAGTACTCATGCCGGATGCGGCATGTGTTAGCTTCTTTGCCTGAAAATCGCGCAAATTCTatatat                                                                                       <  2:553192/143‑1 (MQ=255)
                                                    gcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATAtt                                                                                      >  1:1038108/1‑143 (MQ=255)
                                                     ccATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTg                                                                                     <  2:741976/143‑1 (MQ=255)
                                                               aGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                           <  2:381418/143‑1 (MQ=255)
                                                                ggTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGTGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGTTGTTTGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTAGATATTGCAGAGATCATg                                                                          <  1:235047/143‑1 (MQ=255)
                                                                   cAGGGAATTAAAACCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATTCTCCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAg                                                                       <  1:1018599/143‑1 (MQ=255)
                                                                    aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTCTATATTGCAGAGATCATGTAgg                                                                      >  1:133193/1‑143 (MQ=255)
                                                                    aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGCGATCATGTAgg                                                                      >  1:1022539/1‑143 (MQ=255)
                                                                    aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAgg                                                                      >  1:550533/1‑143 (MQ=255)
                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATACGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCATTCCTGAAACTCGCGAAAATTCTATATATTACAGAGATCATGTAGGc                                                                     >  2:903325/1‑143 (MQ=255)
                                                                     gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGc                                                                     >  1:1025823/1‑143 (MQ=255)
                                                                            aaaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                                              >  2:193944/1‑143 (MQ=255)
                                                                                aCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgc                                                          >  1:563403/1‑142 (MQ=255)
                                                                                 cTAACCCTTACCTTTAAAGTGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAGCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaa                                                         <  2:482868/143‑1 (MQ=255)
                                                                                 cTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCCCGCAATTTCTATATATTGCAGAGATCCTGTAGGCCTGATaagcgaa                                                         >  1:940813/1‑143 (MQ=255)
                                                                                   aaCCCTGACCTGTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaagc                                                       <  2:935216/143‑1 (MQ=255)
                                                                                         tACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCGTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCa                                                 <  1:427419/143‑1 (MQ=255)
                                                                                          aCCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAg                                                <  2:715935/143‑1 (MQ=255)
                                                                                            cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGc                                              >  1:902867/1‑143 (MQ=255)
                                                                                                   gCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGCTAAGCGAAGCGCATCCGGCACTGTTa                                       >  2:1079166/1‑143 (MQ=255)
                                                                                                           caTTTATAACGCATAAGAGTGCAGTACTCAGGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                               <  1:623783/142‑1 (MQ=255)
                                                                                                            aTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg                              >  1:1114191/1‑143 (MQ=255)
                                                                                                             tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGGGGCATGAGTACCATATCCTTCATGAAAATCGCGCAAATTCGATATATTGCAGAGATCATGTATGCCTGATCAGAGAAGCGCCACAGGCAATGTTACAAAAAAAgg                             >  2:9843/1‑142 (MQ=38)
                                                                                                             tggatatCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCAGATCCATCCTGAAAATAGCGCAAATTCTATATATTGCAGAGATCATGTAGACCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGc                             <  2:279787/140‑1 (MQ=38)
                                                                                                               tataTCGAATAAGAGTGCAGTACTCAGGCCGGATGCGGCATGAGTACCATATCCTTCCGGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCa                           <  1:135927/143‑1 (MQ=255)
                                                                                                                 tatCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACg                         <  1:789976/143‑1 (MQ=255)
                                                                                                                     gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                     >  1:536577/1‑143 (MQ=255)
                                                                                                                       aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaa                   >  1:997188/1‑143 (MQ=255)
                                                                                                                       aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaa                   >  2:921209/1‑143 (MQ=255)
                                                                                                                            agTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCACTGTTACCAAAAAAGCCACGGTATAAACCGt              >  1:1030137/1‑143 (MQ=255)
                                                                                                                             gcgCAGTACGCCTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg             <  1:62512/141‑1 (MQ=255)
                                                                                                                             gTGGAGTAGTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg             <  2:804448/143‑1 (MQ=255)
                                                                                                                                cAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTGGct          <  1:401758/143‑2 (MQ=255)
                                                                                                                                  tgacgcATGCCGGAGCCGGCATTAGTACCCTATCCTTCCTGAACACCGCGCAAAGTCTATATACTCCAGCGATCAGGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaa        <  2:94113/138‑1 (MQ=255)
                                                                                                                                  gTACTCATGCCGCATCCGGCATGCGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATCTCTTGCCGAGATCATGTAGGCCTGATAAGCGACGCGCATCAGCCCATGTTACAAAAAAAGCCACCGCATAAACCGTGGccac        >  2:256153/1‑140 (MQ=37)
                                                                                                                                  gTACTCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaa        >  1:171643/1‑143 (MQ=255)
                                                                                                                                   tACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaaa       >  1:1026152/1‑143 (MQ=255)
                                                                                                                                    aCTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt      <  1:678808/143‑1 (MQ=255)
                                                                                                                                     cTCATGCCGCATCCGGCATGAGTACAATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCAACAGGCAATGTTACAAAAAAAGAAACGGGATAAACCGAGGCAAAATc     >  2:289515/1‑143 (MQ=37)
                                                                                                                                      tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCCTGTAGGCCTGCTACGCGACGCGCATCAGGCAATGTTACAAAAAAAGCCACGATATAAACCGTGACAAAATca    >  1:301326/1‑142 (MQ=255)
                                                                                                                                      tCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATCAGCGAAGCGCATCAGGCACTGTTACAAAAAACGCCACGGTATACACCGTGGCAAAAAcc    >  1:742192/1‑143 (MQ=37)
                                                                                                                                        gtGCTGGATGCGGCATGTGTACCGTGTCCTTCATTACAATCGTGCAAATTCTATATGTTGCCGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa  <  2:785227/142‑1 (MQ=37)
                                                                                                                                        aTGGCGGATGCGGCATGAGTACCCTATCCGTCCTGAGAATCGCGCTAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa  <  2:156599/143‑1 (MQ=38)
                                                                                                                                              |                                                                                                                                         
GCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAA  >  NC_000913/542949‑543227

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