Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 532,380 A→G 23.3% intergenic (‑159/‑71) allS ← / → allA allD operon transcriptional activator/ureidoglycolate lyase, releasing urea

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913532,3800AG23.3% 252.2 / 86.4 155intergenic (‑159/‑71)allS/allAallD operon transcriptional activator/ureidoglycolate lyase, releasing urea
Reads supporting (aligned to +/- strand):  ref base A (49/68);  new base G (17/19);  total (68/87)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

Only 100 of 156 total aligned reads displayed.

AGTATTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAGT  >  NC_000913/532281‑532477
                                                                                                   |                                                                                                 
   aTTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTg                                                                                                <  1:2863010/100‑1 (MQ=255)
         aCATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTg                                                                                            <  1:3525826/98‑1 (MQ=255)
         aCATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGt                                                                                           <  2:82264/99‑1 (MQ=255)
          cATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTa                                                                                         <  1:1924815/100‑1 (MQ=255)
          cATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTa                                                                                         <  1:661697/100‑1 (MQ=255)
           aTGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTAt                                                                                        <  2:2603430/100‑1 (MQ=255)
           aTGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTa                                                                                         <  1:3811978/99‑1 (MQ=255)
           aTGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTAt                                                                                        <  1:3596318/100‑1 (MQ=255)
            tggtAACGCTGTTCAATTTGTGGAGTAAATCTGGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTAt                                                                                        <  2:2199055/96‑1 (MQ=255)
             gATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTg                                                                                      <  1:1024522/100‑1 (MQ=255)
             gATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTg                                                                                      <  1:254822/100‑1 (MQ=255)
              aTAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGa                                                                                     <  2:4002564/100‑1 (MQ=255)
                aaCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGaaa                                                                                   <  2:1782041/100‑1 (MQ=255)
                aaCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGaaa                                                                                   <  2:3451104/100‑1 (MQ=255)
                 aCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTTTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  >  1:663115/1‑100 (MQ=255)
                 aCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  >  1:1546572/1‑100 (MQ=255)
                 aCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  >  1:2611758/1‑100 (MQ=255)
                 aCCCTGTTCAATTTGTGGACTAAAAATAGATATAAAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  <  2:535378/100‑1 (MQ=255)
                 aCCCTGTTCAATCTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  >  1:2711453/1‑100 (MQ=255)
                 aCCCTGTTAAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  <  2:1934017/100‑1 (MQ=255)
                  cccTGTTTAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAcc                                                                                 <  2:3262502/100‑1 (MQ=255)
                  cccTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAAcc                                                                                 >  2:1536377/1‑100 (MQ=255)
                   ccTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCa                                                                                >  1:3092827/1‑100 (MQ=255)
                   ccTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCa                                                                                <  2:2633479/100‑1 (MQ=255)
                    cTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCATCTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGaaa                                                                                   <  2:3688714/96‑1 (MQ=255)
                     tggtCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGAGGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAAc                                                                                  <  2:1468828/93‑1 (MQ=255)
                     tGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAg                                                                              <  2:2536290/100‑1 (MQ=255)
                      gTTCAATTTGTGGACTAAATCTAGTTTTGGCAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGc                                                                             <  2:4126036/100‑1 (MQ=255)
                      gTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAgn                                                                             <  2:2136126/100‑2 (MQ=255)
                      gTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGc                                                                             >  1:3998652/1‑100 (MQ=255)
                       ttCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCt                                                                            >  1:2347710/1‑100 (MQ=255)
                         cAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGa                                                                          >  2:1082571/1‑100 (MQ=255)
                          aaTTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTTATGTTTTTCTCTTAAGGTTTTAGATTGTCTGGTATTGAAACCAATCTGAc                                                                         >  2:1023840/1‑100 (MQ=255)
                           aTTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGAcc                                                                        <  1:535116/100‑1 (MQ=255)
                            tttGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCg                                                                       >  2:2218678/1‑100 (MQ=255)
                               gtgGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTc                                                                    <  1:1694241/100‑1 (MQ=255)
                               gtgGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTACTCTTAAGTGTTTAGCTTGCCTGTGATTGAAACAAAGCTGACCGGt                                                                     >  2:2791654/1‑99 (MQ=255)
                                 ggACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcgg                                                                  >  1:2840751/1‑100 (MQ=255)
                                 ggACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTcg                                                                   >  2:4007182/1‑99 (MQ=255)
                                  gACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTCTTGATGTTGTTCTCTTAGGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTcggc                                                                 >  1:1606072/1‑100 (MQ=255)
                                  gACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAGGCTGACCGGTcggc                                                                 <  2:488320/100‑1 (MQ=255)
                                   aCTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTcggcg                                                                <  1:103305/100‑1 (MQ=255)
                                    cTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcggcgg                                                               <  2:1164425/100‑1 (MQ=255)
                                    cTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTcggcgg                                                               >  2:1227350/1‑100 (MQ=255)
                                    cTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTcggcgg                                                               <  1:1169118/100‑1 (MQ=255)
                                     tAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcggcgg                                                               <  2:953957/99‑1 (MQ=255)
                                       aaTCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCggtgg                                                            <  1:3026380/100‑1 (MQ=255)
                                        aTCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCggtggt                                                           <  1:2892576/100‑1 (MQ=255)
                                         tCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGtt                                                          <  1:1470617/100‑1 (MQ=255)
                                         tCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGtt                                                          <  1:3432152/100‑1 (MQ=255)
                                          cTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTg                                                         <  1:2923416/100‑1 (MQ=255)
                                          cTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTg                                                         >  1:2148751/1‑100 (MQ=255)
                                            aGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGaa                                                       >  1:2362281/1‑100 (MQ=255)
                                              ttttGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACg                                                     <  1:1625606/100‑1 (MQ=255)
                                              ttttGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACg                                                     <  2:1871704/100‑1 (MQ=255)
                                                 tgaaaaaaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACg                                                     <  2:1852616/94‑1 (MQ=255)
                                                 tGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGaa                                                  <  2:3762940/100‑1 (MQ=255)
                                                  ggAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTTAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAt                                                 >  2:2511042/1‑100 (MQ=255)
                                                   gAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAtt                                                >  1:2394307/1‑100 (MQ=255)
                                                    aaaaaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTa                                               <  2:2362268/100‑1 (MQ=255)
                                                      aaaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCTGCGGTGGTTGAACGGAATTATg                                             <  2:2397284/100‑1 (MQ=255)
                                                       aaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGt                                            <  1:2440410/100‑1 (MQ=255)
                                                        atatTCTAACTTTGGTATTGTTGTTGTTCTCTTAAGGTTTTATATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTa                                               <  2:3618978/96‑1 (MQ=255)
                                                        atatTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTTACCGGTCGGCGGTGGTTGAACGGAATTATGtt                                           <  2:1669404/100‑1 (MQ=255)
                                                         tatTTCATCTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGTTTGTTATTGTAACCAAGCTGACCGGTCGGCGGTGGTTGAACCGAACAATGt                                            <  2:305469/98‑1 (MQ=255)
                                                            tctaaattttGTAGTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACaa                                       <  2:584281/94‑1 (MQ=255)
                                                            tCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACaa                                       >  2:3248280/1‑100 (MQ=255)
                                                            tCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACaa                                       >  1:3648039/1‑100 (MQ=255)
                                                             ccAACTTTTGTCTTGATGTTGTTCTCTTAAGGCTTCAGATTGCCTGTTATTGCAACCACTCTGACCGGTCGGCGCCGGTTGACCGGAATTATGCTACaa                                       >  2:467708/1‑99 (MQ=255)
                                                             ccAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                      >  1:1445410/1‑100 (MQ=255)
                                                               aaCTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATGATGTTACAAGGa                                    <  1:2708759/100‑1 (MQ=255)
                                                                 cTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACa                                  >  1:1526672/1‑100 (MQ=255)
                                                                   tttGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACaaa                                >  2:638399/1‑100 (MQ=255)
                                                                       tATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGa                            >  1:3781723/1‑100 (MQ=255)
                                                                       tATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGa                            >  1:3119054/1‑100 (MQ=255)
                                                                          tGATGTTGTTCTCTTAAGTTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATgn                         <  2:1596998/100‑2 (MQ=255)
                                                                          tGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGa                         <  2:4129782/100‑1 (MQ=255)
                                                                           gATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGaa                        <  2:1289085/100‑1 (MQ=255)
                                                                            aTGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGaaa                       >  1:2297966/1‑100 (MQ=255)
                                                                             tgttgtTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAAc                      <  1:1184919/100‑1 (MQ=255)
                                                                             tgttgtTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAAc                      <  2:2755442/100‑1 (MQ=255)
                                                                                 gttCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGCAACttga                  >  2:3623268/1‑98 (MQ=255)
                                                                                 gttCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCa                  >  1:3914306/1‑100 (MQ=255)
                                                                                  ttCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAg                 >  1:3417223/1‑100 (MQ=255)
                                                                                   tctttttacgTTTTAGATTTCTTGTTATTGAAACAAATCGGACTGGTCGGCGGTGGTTGAACCGAATTATGTTACAAGGTCAAAAAGATGAAACTTCAg                 <  2:1242843/90‑1 (MQ=255)
                                                                                   tctcTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAgg                <  1:2187280/100‑1 (MQ=255)
                                                                                   tctcTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAgg                <  1:3276376/100‑1 (MQ=255)
                                                                                    ctctTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGt               >  1:1528873/1‑100 (MQ=255)
                                                                                    ctctTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGt               >  1:1749754/1‑100 (MQ=255)
                                                                                     tctTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTa              >  1:3250288/1‑100 (MQ=255)
                                                                                       ttAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTAtt            >  1:3906690/1‑100 (MQ=255)
                                                                                        tAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTa           >  1:2597199/1‑100 (MQ=255)
                                                                                        tAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTa           >  2:883602/1‑100 (MQ=255)
                                                                                        tAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTa           >  2:1651083/1‑100 (MQ=255)
                                                                                          aGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTAcc         <  1:2277808/100‑1 (MQ=255)
                                                                                          aGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATGATGTTATAAGGACAAAAAGATGAAACTTCAGGTATTAcc         >  2:2902638/1‑100 (MQ=255)
                                                                                            gTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGt       >  1:3939276/1‑100 (MQ=255)
                                                                                              tttAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTa     >  1:1693896/1‑100 (MQ=255)
                                                                                              tttAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTa     <  1:974394/100‑1 (MQ=255)
                                                                                                 aGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAGt  >  1:1812526/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
AGTATTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAGT  >  NC_000913/532281‑532477

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