Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,827,816 C→A 7.1% 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 7.1% 396.0 / 16.4 113L148L (CTG→CTTastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base C (59/46);  new base A (5/3);  total (64/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 123 total aligned reads displayed.

GCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGC  >  NC_000913/1827668‑1827964
                                                                                                                                                    |                                                                                                                                                    
gCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAAc                                                                                                                                                      >  4:59789/1‑149 (MQ=255)
     tAGACTAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGa                                                                                                                                                 <  4:439715/149‑1 (MQ=255)
      aGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGa                                                                                                                                                 <  1:653386/148‑1 (MQ=255)
         cAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTc                                                                                                                                              <  1:531758/148‑1 (MQ=255)
          aaaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                             >  1:356962/1‑148 (MQ=255)
          aaaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAttttt                                                                                                                                            >  3:270006/1‑146 (MQ=255)
           aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                            <  2:846824/148‑1 (MQ=255)
           aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                            >  2:836416/1‑148 (MQ=255)
           aaaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAg                                                                                                                                            >  2:260323/1‑148 (MQ=255)
            aaaGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTGTTAAAAGCGAtttttaa                                                                                                                                          <  4:270006/149‑6 (MQ=255)
              aGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGa                                                                                                                                         <  2:210278/148‑1 (MQ=255)
               gTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                       <  4:110884/149‑1 (MQ=255)
                  gCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgc                                                                                                                                    <  3:366910/149‑1 (MQ=255)
                     tACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGgcgc                                                                                                                                  >  1:171543/1‑148 (MQ=255)
                       cccGGTTCACCGTAATGACCGCCTAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                >  1:592114/1‑148 (MQ=255)
                       cccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCCTTTAACAGCGATTCAGTGACGGGCGCtt                                                                                                                                <  1:724551/148‑1 (MQ=255)
                        ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                       <  3:460673/140‑1 (MQ=255)
                        ccGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                       >  4:460673/1‑140 (MQ=255)
                           gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAg                                                                                                                            <  1:531508/148‑1 (MQ=255)
                           gTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGc                                                                                                                           >  4:431778/1‑149 (MQ=255)
                                ccGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAAc                                                                                                                       >  1:877638/1‑148 (MQ=255)
                                 cGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACgg                                                                                                                     <  3:431778/149‑1 (MQ=255)
                                 cGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGCTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACg                                                                                                                      >  1:760026/1‑148 (MQ=255)
                                    aaTGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                             <  4:850938/122‑1 (MQ=255)
                                    aaTGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCa                                                                                                                                             >  3:850938/1‑122 (MQ=255)
                                      tGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                                >  4:794712/1‑149 (MQ=255)
                                       gACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaaa                                                                                                                >  2:625890/1‑148 (MQ=255)
                                         ccgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAtt                                                                                                              >  2:300564/1‑148 (MQ=255)
                                          cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAttt                                                                                                             <  2:910772/148‑1 (MQ=255)
                                          cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAAttt                                                                                                             <  2:592681/148‑1 (MQ=255)
                                          cgccgAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTa                                                                                                            <  4:83352/149‑1 (MQ=255)
                                           gccgAGACGATTGTGGTTTGCCGCCCCCTCACCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCACTGACGGGCGCTTCCAGCGAACGGTGAAGTTTa                                                                                                            <  1:260323/148‑1 (MQ=255)
                                               agaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                       <  1:860279/148‑1 (MQ=255)
                                                gaCGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGtt                                                                                                       >  2:782991/1‑148 (MQ=255)
                                                  cGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAg                                                                                                    <  3:708137/149‑1 (MQ=255)
                                                     ttGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                                 >  3:325357/1‑149 (MQ=255)
                                                      tgtgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                                 >  1:156563/1‑148 (MQ=255)
                                                      tgtgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                                 >  2:508825/1‑148 (MQ=255)
                                                        tgGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                                              >  4:135326/1‑149 (MQ=255)
                                                          gTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCa                                                                                             <  2:787633/148‑1 (MQ=255)
                                                           tttGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCaa                                                                                            >  1:93102/1‑148 (MQ=255)
                                                             tGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACgg                                                                                         <  3:135326/149‑1 (MQ=255)
                                                               ccgccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTAt                                                                                                           >  1:176443/1‑129 (MQ=255)
                                                                 gccCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                                      >  1:774734/1‑148 (MQ=255)
                                                                    cccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                                  >  4:713949/1‑149 (MQ=255)
                                                                     ccTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                                 >  4:723705/1‑149 (MQ=255)
                                                                      cTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGc                                                                                 <  1:782991/148‑1 (MQ=255)
                                                                       tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaa                                                                                                                 <  2:44640/115‑1 (MQ=255)
                                                                       tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGaa                                                                                                                 >  1:44640/1‑115 (MQ=255)
                                                                       tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCa                                                                                >  2:554121/1‑148 (MQ=255)
                                                                       tcatcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCAc                                                                               <  3:218889/149‑1 (MQ=255)
                                                                        catcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTGAATTTCTCTTCGACGTTAGAAATCGCTTTTGACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTGATTGTTCAGTTTGGCACCGGTGTGATGCACt                                                                              >  3:559119/1‑149 (MQ=255)
                                                                          tcaCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACttt                                                                            <  4:567457/149‑1 (MQ=255)
                                                                            aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTg                                                                           >  1:523289/1‑148 (MQ=255)
                                                                            aCCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCATCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGc                                                                          <  4:81033/149‑1 (MQ=255)
                                                                             ccGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                                                  >  2:658139/1‑140 (MQ=255)
                                                                                   aaCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                                                   >  4:481729/1‑149 (MQ=255)
                                                                                    aCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                                                   <  1:508825/148‑1 (MQ=255)
                                                                                    aCGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGc                                                                  >  3:266163/1‑149 (MQ=255)
                                                                                     cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                                 <  1:46545/117‑1 (MQ=255)
                                                                                     cGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGtt                                                                                                 >  2:46545/1‑117 (MQ=255)
                                                                                             gtgGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGc                                                          >  2:108465/1‑148 (MQ=255)
                                                                                                gCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGa                                                       >  2:242299/1‑148 (MQ=255)
                                                                                                 cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTg                                                                                       >  4:291846/1‑115 (MQ=255)
                                                                                                 cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTg                                                                                       <  3:291846/115‑1 (MQ=255)
                                                                                                 cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGAt                                                      <  1:836416/148‑1 (MQ=255)
                                                                                                 cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGAt                                                      >  1:812522/1‑148 (MQ=255)
                                                                                                 cAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGAt                                                      <  1:671885/148‑1 (MQ=255)
                                                                                                   aCGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgc                                                   >  4:491352/1‑149 (MQ=255)
                                                                                                      ccGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGgcgc                                                 >  1:399971/1‑148 (MQ=255)
                                                                                                                gCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACg                                                                                                                                       >  3:465262/1‑52 (MQ=255)
                                                                                                                  tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTg                                         >  3:495160/1‑144 (MQ=255)
                                                                                                                  tAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTg                                         <  4:495160/144‑1 (MQ=255)
                                                                                                                       ttCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTgg                                >  1:660786/1‑148 (MQ=255)
                                                                                                                               gtcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACCGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCcaca                       >  4:184616/1‑149 (MQ=255)
                                                                                                                                tcgtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCcaca                       <  1:902577/148‑1 (MQ=255)
                                                                                                                                  gtTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTg                    <  3:794712/149‑1 (MQ=255)
                                                                                                                                   ttAAAGATCGCTTTTAACAGCGATTCGGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATGCTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTgg                   >  4:574029/1‑149 (MQ=255)
                                                                                                                                   ttAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTgg                   >  3:312469/1‑149 (MQ=255)
                                                                                                                                       gaaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGa                >  1:605232/2‑148 (MQ=255)
                                                                                                                                        aaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGaa               >  1:102787/1‑148 (MQ=255)
                                                                                                                                        aaTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAg              >  3:672684/1‑149 (MQ=255)
                                                                                                                                         aTCGCTGTTAAAAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGc             <  4:728565/149‑1 (MQ=255)
                                                                                                                                          tCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGaa               <  4:246357/146‑1 (MQ=255)
                                                                                                                                              ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc        >  4:397029/1‑149 (MQ=255)
                                                                                                                                              ttttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc        >  3:229227/1‑149 (MQ=255)
                                                                                                                                               tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc        >  1:547541/1‑148 (MQ=255)
                                                                                                                                               tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc        >  2:262343/1‑148 (MQ=255)
                                                                                                                                               tttAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACt       >  3:6769/1‑149 (MQ=255)
                                                                                                                                                 tAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTg      <  2:124325/148‑1 (MQ=255)
                                                                                                                                                  aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGa     <  2:156563/148‑1 (MQ=255)
                                                                                                                                                   aCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCGCATTGGCGAAGCGGAACTGACg   >  3:513201/1‑149 (MQ=255)
                                                                                                                                                    cAGCGATTCGGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTg                                                                           <  2:849948/76‑1 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGATGGGCGCTTCCAGCGAATGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACGTTGCCATCCAGCGTATCGGCAGATGGCGCGATTGTTACCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGc  <  3:558154/149‑1 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTg                                                                           >  1:849948/1‑76 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTg                                                                           >  1:442883/1‑76 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                                    >  2:6599/1‑115 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCa                                    <  1:6599/115‑1 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGc  <  3:59789/149‑1 (MQ=255)
                                                                                                                                                    cAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGc  <  4:266163/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCCGTAGACAAAAAGTTGCATACCCGGTTCACCGTAATGACCGCCGAGACGATTGTGGTTTGCCGCCCCCTCATCACCGAGCAACGCTACCTGTGGCAACGCCGAATGGACGCTAAATTTCTCTTCGTCGTTAAAAATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGC  >  NC_000913/1827668‑1827964

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: