Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,488,983 C→A 6.9% V400L (GTG→TTG)  argD ← N‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,488,9830CA 6.9% 303.5 / 2.2 93V400L (GTG→TTG) argDN‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base A (6/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTG  >  NC_000913/3488836‑3489131
                                                                                                                                                   |                                                                                                                                                    
tCAGCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTATGCCCCAACCACCTTCTCCACcg                                                                                                                                                    <  2:706236/150‑1 (MQ=255)
tCAGCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACcg                                                                                                                                                    <  2:894115/150‑1 (MQ=255)
    cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                <  1:198086/150‑1 (MQ=255)
    cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                <  1:1281765/150‑1 (MQ=255)
     gACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                               >  1:785404/1‑150 (MQ=255)
      aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTggg                                                                                                                                              <  1:1428323/150‑1 (MQ=255)
      aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgcg                                                                                                                                              <  2:545682/150‑3 (MQ=255)
        cGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCg                                                                                                                                            <  2:310937/150‑1 (MQ=255)
         gCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGa                                                                                                                                           >  2:1544800/1‑150 (MQ=255)
          cATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGaa                                                                                                                                          >  1:1022940/1‑150 (MQ=255)
             cACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggt                                                                                                                                       <  2:774061/150‑4 (MQ=255)
             cACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACg                                                                                                                                       <  2:1518737/150‑1 (MQ=255)
               cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                     >  2:1200636/1‑150 (MQ=255)
               cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                     >  2:1318466/1‑150 (MQ=255)
               cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                     >  1:742424/1‑150 (MQ=255)
                gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtggt                                                                                                                                    <  1:974122/150‑7 (MQ=255)
                gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                    <  1:1544800/150‑1 (MQ=255)
                  gATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCa                                                                                                                                  <  1:787373/150‑1 (MQ=255)
                    tttGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATc                                                                                                                                <  2:1246158/150‑1 (MQ=255)
                      tGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATccc                                                                                                                              >  2:408688/1‑150 (MQ=255)
                      tGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATccc                                                                                                                              >  1:306002/1‑150 (MQ=255)
                        gCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCtt                                                                                                                            <  2:491841/150‑1 (MQ=255)
                         cTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTc                                                                                                                           >  1:1328320/1‑150 (MQ=255)
                             gTCGCAAGTTGCGGGATTCGAACGCGTAATGCAGGCTAAACCATCGGATAAGGCAATCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTCGGCGAAACGATGCATCCCTACATCg                                                                                                                       >  1:112926/1‑150 (MQ=255)
                              tCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGa                                                                                                                      >  2:265701/1‑150 (MQ=255)
                                gCAAGTTGCGGGATTCGAAGGCGTACTGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGCCCTCTTTTGCCTGATGCTTCGTATTTACGCCCCACCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCCTCGata                                                                                                                    >  2:1058732/1‑150 (MQ=255)
                                gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGata                                                                                                                    <  2:981277/150‑1 (MQ=255)
                                gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGata                                                                                                                    <  2:975341/150‑1 (MQ=255)
                                gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGata                                                                                                                    <  1:163167/150‑1 (MQ=255)
                                gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGata                                                                                                                    >  1:52791/1‑150 (MQ=255)
                                      tGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCa                                                                                                              <  2:640899/150‑1 (MQ=255)
                                           gATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTc                                                                                                         <  1:1200636/150‑1 (MQ=255)
                                            aTTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTcg                                                                                                        <  2:1506441/150‑2 (MQ=255)
                                            aTTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTcc                                                                                                        <  1:726839/150‑1 (MQ=255)
                                               cGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTccacc                                                                                                     >  2:305674/1‑150 (MQ=255)
                                                      gTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAAACACCTTCGCCACCGCGTGGGCGAAACGTTGCCTCCCTTCATCGATATCCGCATCTTCCCCCACCAGCg                                                                                              >  2:892945/1‑150 (MQ=255)
                                                          tGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACgg                                                                                          >  2:314280/1‑150 (MQ=255)
                                                                cTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaaa                                                                                    >  1:1427694/1‑150 (MQ=255)
                                                                  aaaCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACGGCGTTTGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAAc                                                                                  <  2:654857/150‑1 (MQ=255)
                                                                    aCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGc                                                                                <  1:903979/150‑1 (MQ=255)
                                                                     ccATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCa                                                                               >  2:1566986/1‑150 (MQ=255)
                                                                       aTCGGATAAGGGATTCAGGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATc                                                                             <  2:112926/150‑1 (MQ=255)
                                                                          ggATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATcaca                                                                          <  1:592480/150‑1 (MQ=255)
                                                                          ggATAAGGCATTCACGCCGCATCCCACAGCTTTTGCCTGATGCTTCGTATTTACGCCCCCACCACCTTCGCCACCGCGAGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATATTCCAACACCACCGACGGTGCAAAACGCATcaca                                                                          >  2:116760/1‑150 (MQ=255)
                                                                             tAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtcc                                                                       <  1:352612/150‑1 (MQ=255)
                                                                                ggCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccggt                                                                    >  1:1499081/1‑150 (MQ=255)
                                                                                     tCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGc                                                               <  1:314280/150‑1 (MQ=255)
                                                                                     tCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCAACATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGc                                                               >  2:713510/1‑150 (MQ=255)
                                                                                         gccgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTc                                                           <  2:107064/150‑1 (MQ=255)
                                                                                           cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAg                                                         <  1:509265/150‑1 (MQ=255)
                                                                                                 cGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCAt                                                   <  1:168946/150‑1 (MQ=255)
                                                                                                   aCATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTa                                                 <  1:305674/150‑1 (MQ=255)
                                                                                                              cTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctct                                      >  2:1046424/1‑150 (MQ=255)
                                                                                                               tGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTg                                     <  2:1231680/150‑1 (MQ=255)
                                                                                                                  tGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTgcgc                                  >  2:887458/1‑150 (MQ=255)
                                                                                                                  tGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTgcgc                                  <  1:675550/150‑1 (MQ=255)
                                                                                                                   gCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGcc                                 >  1:1374730/1‑150 (MQ=255)
                                                                                                                    cTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGccc                                >  2:391518/1‑150 (MQ=255)
                                                                                                                     ttCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCg                               <  1:1333082/150‑1 (MQ=255)
                                                                                                                      tCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGc                              >  2:118716/1‑150 (MQ=255)
                                                                                                                       cGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCa                             >  2:1525906/1‑150 (MQ=255)
                                                                                                                        gTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCAt                            <  1:24349/150‑1 (MQ=255)
                                                                                                                         tATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaa                                                                                     <  2:817956/92‑1 (MQ=255)
                                                                                                                         tATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaa                                                                                     >  1:817956/1‑92 (MQ=255)
                                                                                                                           tttACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGCCGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACa                         >  1:542781/1‑150 (MQ=255)
                                                                                                                               cGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaa                     >  1:577777/1‑150 (MQ=255)
                                                                                                                               cGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaa                     <  2:923172/150‑1 (MQ=255)
                                                                                                                               cGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaa                     >  1:725319/1‑150 (MQ=255)
                                                                                                                                 ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGa                      >  1:1234479/1‑147 (MQ=255)
                                                                                                                                 ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGa                      <  2:1234479/147‑1 (MQ=255)
                                                                                                                                 ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                   <  2:785404/150‑1 (MQ=255)
                                                                                                                                   ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 <  2:742424/150‑1 (MQ=255)
                                                                                                                                   ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 <  2:955032/150‑1 (MQ=255)
                                                                                                                                   ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 <  2:262952/150‑1 (MQ=255)
                                                                                                                                   ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCACATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 <  1:713510/150‑1 (MQ=255)
                                                                                                                                     gcaacgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg               >  2:1487781/8‑150 (MQ=255)
                                                                                                                                     aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg               >  1:1493080/1‑150 (MQ=255)
                                                                                                                                      accaccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgc              >  2:1561308/1‑150 (MQ=255)
                                                                                                                                      accaccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgc              >  1:268689/1‑150 (MQ=255)
                                                                                                                                       ccaccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgcg             <  2:52791/150‑1 (MQ=255)
                                                                                                                                        acgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAACCCGCCTCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGACATCACgcgc            >  1:639051/5‑150 (MQ=255)
                                                                                                                                        acgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCTGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACgcgc            >  2:1316753/5‑150 (MQ=255)
                                                                                                                                         cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCTCCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc           >  1:1140190/4‑150 (MQ=255)
                                                                                                                                         cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc           >  1:324175/4‑150 (MQ=255)
                                                                                                                                         accTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc           >  1:881654/1‑150 (MQ=255)
                                                                                                                                         accTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc           >  2:1362365/1‑150 (MQ=255)
                                                                                                                                          gtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGccc          >  1:1607847/3‑150 (MQ=255)
                                                                                                                                            ttCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGa        <  2:1022940/150‑1 (MQ=255)
                                                                                                                                            ttCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGa        >  1:1159104/1‑150 (MQ=255)
                                                                                                                                               gCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCt     >  1:979266/1‑150 (MQ=255)
                                                                                                                                                 cACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCttt   >  2:795595/1‑150 (MQ=255)
                                                                                                                                                  aCCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCTTTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCGCGACCTTTg  <  2:1328320/150‑1 (MQ=255)
                                                                                                                                                  aCCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTg  >  1:938584/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAGCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTG  >  NC_000913/3488836‑3489131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: