Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,811 T→G 7.9% 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 7.9% 439.0 / 26.7 127K150T (AAA→ACA) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (62/55);  new base G (6/4);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 141 total aligned reads displayed.

CGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAC  >  NC_000913/1827667‑1827958
                                                                                                                                                |                                                                                                                                                   
cGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTaa                                                                                                                                                 <  3:157058/149‑1 (MQ=255)
  cccGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                               >  1:723018/1‑149 (MQ=255)
      tAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                           >  1:1743293/1‑149 (MQ=255)
         aCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                        >  1:1137722/1‑149 (MQ=255)
         aCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                        >  4:503688/1‑149 (MQ=255)
          cAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                       <  1:497210/149‑1 (MQ=255)
               aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTaa                                                                                                                                                 <  1:1057882/134‑1 (MQ=255)
               aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTaa                                                                                                                                                 >  2:1057882/1‑134 (MQ=255)
               aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGAc                                                                                                                                  >  3:1077225/1‑149 (MQ=255)
                        cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                                                                           <  2:1768540/147‑1 (MQ=255)
                        cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                                                                           >  1:1768540/1‑147 (MQ=255)
                        cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                         >  1:870689/1‑149 (MQ=255)
                         ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCttt                                                                                                                                                    <  4:1331371/121‑1 (MQ=255)
                         ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                        >  2:201159/1‑149 (MQ=255)
                         ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTcc                                                                                                                        >  2:552023/1‑149 (MQ=255)
                         ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaac                                                                                                                                   >  2:1852210/1‑132 (MQ=255)
                            gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGc                                                                                                                     >  4:1343683/1‑149 (MQ=255)
                                aCCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAAc                                                                                                                 >  4:1850812/1‑149 (MQ=255)
                                    tAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTg                                                                                                             <  1:1015481/149‑1 (MQ=255)
                                      aTGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                       >  3:768103/1‑121 (MQ=255)
                                         aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                               >  4:802308/1‑110 (MQ=255)
                                         aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                        >  2:420166/1‑149 (MQ=255)
                                         aCCGCCGAGACGATTGTGGTTTACCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTCCTCTTTGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                        >  2:772557/1‑149 (MQ=255)
                                           cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTa                                                                                                      <  3:1343683/149‑1 (MQ=255)
                                            gccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                           <  3:174012/111‑1 (MQ=255)
                                             ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  3:503688/149‑1 (MQ=255)
                                             ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  1:201159/149‑1 (MQ=255)
                                             ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  2:1137722/149‑1 (MQ=255)
                                                   cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAg                                                                                              <  1:420166/149‑1 (MQ=255)
                                                    gATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCttt                                                                                                                                                    <  4:1668045/94‑1 (MQ=255)
                                                    gATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACgg                                                                                                               >  3:1820091/1‑131 (MQ=255)
                                                    gATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACgg                                                                                                               <  4:1820091/131‑1 (MQ=255)
                                                      ttGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                           <  1:317389/149‑1 (MQ=255)
                                                        gtgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                         >  3:657017/1‑149 (MQ=255)
                                                               gccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGt                                                                                  <  2:723018/149‑1 (MQ=255)
                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                         >  3:1802592/1‑93 (MQ=255)
                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                         <  4:264494/93‑1 (MQ=255)
                                                                ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                         >  1:1265044/1‑93 (MQ=255)
                                                                   cccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                              <  4:1077225/149‑1 (MQ=255)
                                                                   cccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                              <  1:552023/149‑1 (MQ=255)
                                                                    ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGAt                                                                             >  4:1753399/1‑149 (MQ=255)
                                                                     cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                            <  4:657017/149‑1 (MQ=255)
                                                                        tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCAc                                                                         >  3:1768002/1‑149 (MQ=255)
                                                                           tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    <  2:321163/119‑1 (MQ=255)
                                                                           tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                    >  1:321163/1‑119 (MQ=255)
                                                                           tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACttt                                                                      <  4:1768002/149‑1 (MQ=255)
                                                                             aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATGGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                         >  3:1592497/1‑96 (MQ=255)
                                                                             aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                                                                    >  4:1765461/1‑149 (MQ=255)
                                                                             aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACACTAAATTCCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                         <  4:1592497/96‑1 (MQ=255)
                                                                                  gCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATcc                                                               >  1:1651376/1‑149 (MQ=255)
                                                                                   cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                         >  2:1735785/1‑122 (MQ=255)
                                                                                   cAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                         <  1:1735785/122‑1 (MQ=255)
                                                                                    aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCa                                                                                                                       <  1:861476/91‑1 (MQ=255)
                                                                                    aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCa                                                                                                                       >  2:861476/1‑91 (MQ=255)
                                                                                    aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                                             <  3:1753399/149‑1 (MQ=255)
                                                                                             tgtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc                                                    <  2:1483579/149‑1 (MQ=255)
                                                                                                 gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGAt                                                <  2:253399/149‑1 (MQ=255)
                                                                                                  cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                             <  1:741923/103‑1 (MQ=255)
                                                                                                   aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc                                                    <  1:1251042/143‑1 (MQ=255)
                                                                                                   aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc                                                    >  2:1251042/1‑143 (MQ=255)
                                                                                                   aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGGGAAATTTAtt                                                                                                    >  3:565279/1‑95 (MQ=255)
                                                                                                    aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAt                                                                                                         <  4:1189915/89‑1 (MQ=255)
                                                                                                      gCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgcgc                                           >  1:240839/1‑149 (MQ=255)
                                                                                                       ccGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCTTATCGGCAGATGGcgcg                                          >  1:491020/1‑149 (MQ=255)
                                                                                                                 gCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                           >  4:1003180/1‑90 (MQ=255)
                                                                                                                  cTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATc                                                                >  3:911638/1‑116 (MQ=255)
                                                                                                                   tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                              >  2:1843756/1‑149 (MQ=255)
                                                                                                                   tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                              >  3:124585/1‑149 (MQ=255)
                                                                                                                    aaaTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAt                             <  4:124585/149‑1 (MQ=255)
                                                                                                                      aTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTg                           >  1:957258/1‑149 (MQ=255)
                                                                                                                       tttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTgg                          >  2:505855/1‑149 (MQ=255)
                                                                                                                        ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                                                                    >  2:918759/1‑106 (MQ=255)
                                                                                                                        ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGc                         <  2:957258/149‑1 (MQ=255)
                                                                                                                             ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                           <  4:1052183/94‑1 (MQ=255)
                                                                                                                             ttCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                           >  3:1052183/1‑94 (MQ=255)
                                                                                                                                gtcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCcaca                 >  3:154004/1‑149 (MQ=255)
                                                                                                                                 tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAtt                            >  2:1226819/1‑137 (MQ=255)
                                                                                                                                 tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAtt                            <  1:1226819/137‑1 (MQ=255)
                                                                                                                                   gtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTg              >  3:1373312/1‑149 (MQ=255)
                                                                                                                                    ttATAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTg              >  1:1620173/1‑148 (MQ=255)
                                                                                                                                    ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACACACATTg              <  2:1620173/148‑1 (MQ=255)
                                                                                                                                       tgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTgg                          >  4:83851/3‑133 (MQ=255)
                                                                                                                                       tgaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTgg                          <  3:83851/131‑1 (MQ=255)
                                                                                                                                        aaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGaa         <  3:1269425/149‑1 (MQ=255)
                                                                                                                                         aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAg        >  3:196308/1‑149 (MQ=255)
                                                                                                                                          aTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGc       <  3:743453/149‑1 (MQ=255)
                                                                                                                                          aTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGc                                  <  2:1711142/122‑1 (MQ=255)
                                                                                                                                           tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCa                                                   <  1:716769/104‑1 (MQ=255)
                                                                                                                                           tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCa                                                   >  2:716769/1‑104 (MQ=255)
                                                                                                                                           tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg      >  1:1466597/1‑149 (MQ=255)
                                                                                                                                           tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATg                                               <  2:555157/108‑1 (MQ=255)
                                                                                                                                           tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATg                                               >  1:555157/1‑108 (MQ=255)
                                                                                                                                           tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg      >  2:773553/1‑149 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc                               <  4:249945/120‑1 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc                               >  3:249945/1‑120 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  <  2:1466597/149‑1 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  >  4:565627/1‑149 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  <  4:1373312/149‑1 (MQ=255)
                                                                                                                                               ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc  >  1:1268859/1‑149 (MQ=255)
                                                                                                                                               ttttAACAGCAATTCAGTGACGGGCGCTTCCGGCGAACGGTGAAATTTATTGTTCAGGTTGGCAATGGTGAGATGCACTTTGCCAACCAGCGTATCGGCAGATGGCGCGATCGTTGCCGGATTGGCTGCCCACATTGGCGAAGCGGAAc  >  3:1168662/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAC  >  NC_000913/1827667‑1827958

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