Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,811 T→G 5.7% K150T (AAA→ACA)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,8110TG 5.7% 315.4 / 9.8 87K150T (AAA→ACA) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base G (3/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

TTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAC  >  NC_000913/1827663‑1827958
                                                                                                                                                    |                                                                                                                                                   
ttctCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                     <  1:121408/149‑1 (MQ=255)
ttctCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                     <  4:749690/149‑1 (MQ=255)
ttctCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                     <  4:286325/149‑1 (MQ=255)
ttctCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                     >  2:121408/1‑149 (MQ=255)
ttctCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCtt                                                                                                                                                     >  2:240702/1‑149 (MQ=255)
   tcGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTa                                                                                                                                                  >  3:1344184/1‑149 (MQ=255)
      cccGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                               <  3:981038/149‑1 (MQ=255)
       ccGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAg                                                                                                                                              <  1:240702/149‑1 (MQ=255)
         gTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCg                                                                                                                                            >  2:1365963/1‑149 (MQ=255)
            gACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtt                                                                                                                                         <  4:1344184/149‑1 (MQ=255)
             aCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAAc                                                                                                                                                >  4:183377/1‑141 (MQ=255)
             aCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAAc                                                                                                                                                <  3:183377/141‑1 (MQ=255)
               aaaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                      >  3:995369/1‑149 (MQ=255)
                     ttGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTaa                                                                                                                                                 >  3:208260/1‑132 (MQ=255)
                     ttGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAACTCCTCTTCGTCGTTAAAAATCGCTTTTaa                                                                                                                                                 <  4:208260/132‑1 (MQ=255)
                       gCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgc                                                                                                                              <  3:1430431/149‑1 (MQ=255)
                         aTACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcgc                                                                                                                            >  4:397647/1‑149 (MQ=255)
                          tACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                                                                           >  3:1395134/1‑149 (MQ=255)
                            cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                         <  3:1305833/149‑1 (MQ=255)
                            cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                         <  3:397647/149‑1 (MQ=255)
                             ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                        >  3:755341/1‑149 (MQ=255)
                                gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAt                                                                                                                                          >  3:278888/1‑128 (MQ=255)
                                gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAt                                                                                                                                          <  4:278888/128‑1 (MQ=255)
                                  tCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGa                                                                                                                   <  4:189068/149‑1 (MQ=255)
                                   cACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGaa                                                                                                                  >  4:654665/1‑149 (MQ=255)
                                    aCCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCttt                                                                                                                                                    >  2:5970/1‑114 (MQ=255)
                                    aCCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCttt                                                                                                                                                    <  1:5970/114‑1 (MQ=255)
                                             aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                        >  4:1180771/1‑149 (MQ=255)
                                                gccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                                                                     <  4:1395134/149‑1 (MQ=255)
                                                 ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  4:755341/149‑1 (MQ=255)
                                                       cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAg                                                                                              <  1:1164772/149‑1 (MQ=255)
                                                                       cccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                              <  4:1292437/149‑1 (MQ=255)
                                                                        ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                        <  3:571642/90‑1 (MQ=255)
                                                                        ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                        >  4:571642/1‑90 (MQ=255)
                                                                        ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGAt                                                                             <  1:1365963/149‑1 (MQ=255)
                                                                          ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTc                                                                                                >  4:816556/1‑128 (MQ=255)
                                                                          ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTc                                                                                                <  3:816556/128‑1 (MQ=255)
                                                                          ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                                        <  2:1378726/136‑1 (MQ=255)
                                                                          ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                                        >  1:1378726/1‑136 (MQ=255)
                                                                               tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    >  2:253300/1‑119 (MQ=255)
                                                                               tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  1:253300/119‑1 (MQ=255)
                                                                                caCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGt                                                                                  >  2:959414/1‑136 (MQ=255)
                                                                                caCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGt                                                                                  <  1:959414/136‑1 (MQ=255)
                                                                                      gCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATcc                                                               >  4:568320/1‑149 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcgc                                                                                                                            <  2:141649/86‑1 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcgc                                                                                                                            >  1:141649/1‑86 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 <  3:503622/113‑1 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 >  4:503622/1‑113 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGt                                                                                  <  2:1367094/128‑1 (MQ=255)
                                                                                        aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGt                                                                                  >  1:1367094/1‑128 (MQ=255)
                                                                                                  gtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                          <  3:238642/94‑1 (MQ=255)
                                                                                                  gtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                          >  4:238642/1‑94 (MQ=255)
                                                                                                     gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                          <  2:616018/91‑1 (MQ=255)
                                                                                                     gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                          >  1:616018/1‑91 (MQ=255)
                                                                                                      cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTa                                                                                                      <  1:903828/94‑1 (MQ=255)
                                                                                                      cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTa                                                                                                      >  2:903828/1‑94 (MQ=255)
                                                                                                       aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATgg                                              >  3:1345311/1‑149 (MQ=255)
                                                                                                       aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATgg                                              >  3:1224880/1‑149 (MQ=255)
                                                                                                        aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 <  3:1407287/97‑1 (MQ=255)
                                                                                                        aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 >  4:1407287/1‑97 (MQ=255)
                                                                                                         cGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgcg                                            <  3:654665/149‑1 (MQ=255)
                                                                                                            cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                        <  3:308955/86‑1 (MQ=255)
                                                                                                            cGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                        >  4:308955/1‑86 (MQ=255)
                                                                                                                 ggACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGtt                                    <  4:995369/149‑1 (MQ=255)
                                                                                                                  gACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                      <  2:1256391/50‑1 (MQ=255)
                                                                                                                  gACGCTAAAGTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                      >  1:1256391/1‑50 (MQ=255)
                                                                                                                       tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgcgc                                           >  1:775737/1‑136 (MQ=255)
                                                                                                                       tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgcgc                                           <  2:775737/136‑1 (MQ=255)
                                                                                                                       tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                              >  3:64941/1‑149 (MQ=255)
                                                                                                                          aTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTg                           >  1:174186/1‑149 (MQ=255)
                                                                                                                           tttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTgg                          >  2:81304/1‑149 (MQ=255)
                                                                                                                               tctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 <  1:105439/74‑1 (MQ=255)
                                                                                                                               tctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                 >  2:105439/1‑74 (MQ=255)
                                                                                                                                ctTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAcc                     <  3:1180771/149‑1 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTg                                                                                 <  3:218333/88‑1 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTg                                                                                 >  4:218333/1‑88 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAtt                            >  2:1362613/1‑141 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAtt                            <  1:1362613/141‑1 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAccc                    <  4:64941/149‑1 (MQ=255)
                                                                                                                                 ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAccc                    >  4:618269/1‑149 (MQ=255)
                                                                                                                                              aTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                                                                          <  3:422745/82‑1 (MQ=255)
                                                                                                                                              aTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                                                                          >  4:422745/1‑82 (MQ=255)
                                                                                                                                              aTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGc       >  2:950484/1‑149 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATc                                                                <  2:295474/91‑1 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATc                                                                >  1:295474/1‑91 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc                                                            >  4:1400890/1‑95 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc                                                            <  3:1400890/95‑1 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGtt                                    <  1:854505/119‑1 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGtt                                    >  2:854505/1‑119 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt               >  4:856070/1‑140 (MQ=255)
                                                                                                                                               tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAtt               <  3:856070/140‑1 (MQ=255)
                                                                                                                                                  cTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa   >  2:115136/1‑149 (MQ=255)
                                                                                                                                                  cTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATCCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa   <  1:115136/149‑1 (MQ=255)
                                                                                                                                                  cTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACCGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa   <  4:1224880/149‑1 (MQ=255)
                                                                                                                                                   ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  >  3:731767/1‑149 (MQ=255)
                                                                                                                                                   ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  >  2:598075/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCTCGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAC  >  NC_000913/1827663‑1827958

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