Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,816 C→A 8.9% L148L (CTG→CTT astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,8160CA 8.9% 477.1 / 29.1 135L148L (CTG→CTTastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base C (66/57);  new base A (7/5);  total (73/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 135 total aligned reads displayed.

CCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGA  >  NC_000913/1827669‑1827961
                                                                                                                                                   |                                                                                                                                                 
cccGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                  >  1:723018/1‑149 (MQ=255)
    tAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                              >  1:1743293/1‑149 (MQ=255)
       aCAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                           >  1:1137722/1‑149 (MQ=255)
        cAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                          <  1:497210/149‑1 (MQ=255)
        cAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                          <  2:1743293/149‑1 (MQ=255)
             aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGAc                                                                                                                                     >  3:1077225/1‑149 (MQ=255)
                      cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                                                                              >  1:1768540/1‑147 (MQ=255)
                      cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCt                                                                                                                              <  2:1768540/147‑1 (MQ=255)
                      cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTc                                                                                                                            >  1:870689/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)
                             cACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGaa                                                                                                                     >  3:1376824/1‑149 (MQ=255)
                                  tAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTg                                                                                                                <  1:1015481/149‑1 (MQ=255)
                                    aTGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                          >  3:768103/1‑121 (MQ=255)
                                       aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                  <  3:802308/110‑1 (MQ=255)
                                       aCCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACa                                                                                                                                                  >  4:802308/1‑110 (MQ=255)
                                       aCCGCCGAGACGATTGTGGTTTACCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTCCTCTTTGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                           >  2:772557/1‑149 (MQ=255)
                                          gccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                              <  3:174012/111‑1 (MQ=255)
                                          gccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                              >  4:174012/1‑111 (MQ=255)
                                           ccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                       <  1:201159/149‑1 (MQ=255)
                                                 cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAg                                                                                                 <  1:420166/149‑1 (MQ=255)
                                                  gATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACgg                                                                                                                  <  4:1820091/131‑1 (MQ=255)
                                                  gATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACgg                                                                                                                  >  3:1820091/1‑131 (MQ=255)
                                                    ttGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                              >  2:317389/1‑149 (MQ=255)
                                                      gtgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTgg                                                                                            >  3:657017/1‑149 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            <  4:264494/93‑1 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            >  3:1802592/1‑93 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            <  2:1265044/93‑1 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            >  3:264494/1‑93 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            >  1:1265044/1‑93 (MQ=255)
                                                              ccgccCCCTCATCACCGAGCAACGCTACCTGAGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGAtt                                                                                                                                            <  4:1802592/93‑1 (MQ=255)
                                                               cgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                   <  4:1869946/149‑1 (MQ=255)
                                                                 cccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                                 <  4:1077225/149‑1 (MQ=255)
                                                                 cccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTgaga                                                                                 >  1:253399/1‑149 (MQ=255)
                                                                  ccccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGAt                                                                                >  4:1753399/1‑149 (MQ=255)
                                                                   cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                               <  4:657017/149‑1 (MQ=255)
                                                                   cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                               >  4:1269425/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)
                                                                         tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACttt                                                                         <  1:1254416/149‑1 (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)
                                                                                  aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                                                <  3:1753399/149‑1 (MQ=255)
                                                                                      cTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTa                                                            >  1:1483579/1‑149 (MQ=255)
                                                                                          cTGTAGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCgg                                                        <  3:1765461/149‑1 (MQ=255)
                                                                                           tgtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc                                                       <  2:1483579/149‑1 (MQ=255)
                                                                                               gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGAt                                                   <  2:253399/149‑1 (MQ=255)
                                                                                                cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                                >  2:741923/1‑103 (MQ=255)
                                                                                                cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAgg                                                                                                <  1:741923/103‑1 (MQ=255)
                                                                                                 aaCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAtt                                                                                                       <  4:565279/95‑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)
                                                                                                       gAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGAt                                           <  2:870689/149‑1 (MQ=255)
                                                                                                           ggACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATGCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGtt                                       <  1:772557/149‑1 (MQ=255)
                                                                                                               gCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                              >  4:1003180/1‑90 (MQ=255)
                                                                                                               gCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCg                                   >  1:686143/1‑149 (MQ=255)
                                                                                                                 tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                                 >  3:124585/1‑149 (MQ=255)
                                                                                                                   aaTTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCAtt                               <  2:1651376/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)
                                                                                                                      ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                                                                       <  1:918759/106‑1 (MQ=255)
                                                                                                                      ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGc                            <  2:957258/149‑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)
                                                                                                                                  ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTgg                <  1:505855/149‑1 (MQ=255)
                                                                                                                                  ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACACACATTg                 <  2:1620173/148‑1 (MQ=255)
                                                                                                                                    aaaaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc                                                               <  1:428449/100‑1 (MQ=255)
                                                                                                                                    aaaaaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc                                                               >  2:428449/1‑100 (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)
                                                                                                                                         tCGCTTTTAACAGCGCTACAGTGAGGGGCGTTTCCCGCGAACGGTGAAATTTATTGTTCAGGTTGGCAGCGGTGAGATGCACTTTGCCAACCAGCGTATCGGCAGATGGCGCTATCGTTACCTGATTGGCTACCCACATTGGCGAAGCg         >  4:1173925/1‑149 (MQ=255)
                                                                                                                                         tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCa                                                      >  2:716769/1‑104 (MQ=255)
                                                                                                                                         tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg         >  4:1042735/1‑149 (MQ=255)
                                                                                                                                         tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATg                                                  >  1:555157/1‑108 (MQ=255)
                                                                                                                                         tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg         >  2:773553/1‑149 (MQ=255)
                                                                                                                                         tCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg         >  3:67041/1‑149 (MQ=255)
                                                                                                                                            ctgttAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATc                                                          >  2:674812/4‑97 (MQ=255)
                                                                                                                                            ctgttAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATc                                                          <  1:674812/94‑1 (MQ=255)
                                                                                                                                            cTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGaa      >  2:1877873/1‑149 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc                                  >  3:249945/1‑120 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc                                  <  4:249945/120‑1 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc     <  4:1373312/149‑1 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc     >  2:611328/1‑149 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc     >  4:565627/1‑149 (MQ=255)
                                                                                                                                             ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc     >  1:1268859/1‑149 (MQ=255)
                                                                                                                                             ttttAACAGCAATTCAGTGACGGGCGCTTCCGGCGAACGGTGAAATTTATTGTTCAGGTTGGCAATGGTGAGATGCACTTTGCCAACCAGCGTATCGGCAGATGGCGCGATCGTTGCCGGATTGGCTGCCCACATTGGCGAAGCGGAAc     >  3:1168662/1‑149 (MQ=255)
                                                                                                                                                tAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGa  <  3:1850812/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGA  >  NC_000913/1827669‑1827961

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