Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,488,989 G→T 7.1% H398N (CAC→AAC)  argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,488,9890GT 7.1% 212.5 / 4.3 85H398N (CAC→AAC) argDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (32/46);  new base T (2/5);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACTGTG  >  NC_000913/3488840‑3489138
                                                                                                                                                     |                                                                                                                                                     
cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                       >  2:1270685/1‑150 (MQ=255)
  aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTggg                                                                                                                                                     >  1:572606/1‑150 (MQ=255)
    cGCATCACGGGATTTGGCTGAGTCTCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCAACGGATAAGGCATTCACGCCGCATCCGACAGCTTTTGCCGGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCg                                                                                                                                                   <  1:150398/150‑1 (MQ=255)
    cGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCg                                                                                                                                                   <  2:136007/150‑1 (MQ=255)
    cGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCAATCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCg                                                                                                                                                   <  2:1166908/150‑1 (MQ=255)
        tCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaagg                                                                                                                                               <  1:1048589/150‑3 (MQ=255)
        tCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTGTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaagg                                                                                                                                               <  1:667750/150‑3 (MQ=255)
         cACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCCTCCGACATCTTTTGCCTGATGCTTCGGATTTACGCCCCAGCCACCTTCGCCACCGCGTGGGCGAAACg                                                                                                                                              >  1:237107/1‑150 (MQ=255)
         cACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggt                                                                                                                                              <  1:178944/150‑4 (MQ=255)
          aCGGTATTTGGCTGAGTCGCAAGTTGCGGGATTCGAATGCGTAATGCATGCTAAACCATCGGATAAGTCATTCACGCCGCATACGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtg                                                                                                                                             <  1:653899/150‑5 (MQ=255)
          aCGGGATTTGGCTGAGTCGCAAGGTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGt                                                                                                                                             <  1:506129/150‑1 (MQ=255)
           cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtgg                                                                                                                                            <  1:1031760/150‑6 (MQ=255)
           cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCCAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGACTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                            >  2:16898/1‑150 (MQ=255)
             ggAATTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTACGCCACCGCGTGGGCGAAACGTTGc                                                                                                                                          <  2:1065391/150‑1 (MQ=255)
                tttGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATc                                                                                                                                       <  1:1288458/150‑1 (MQ=255)
                 ttGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATcc                                                                                                                                      <  1:1192571/150‑1 (MQ=255)
                 ttGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATcc                                                                                                                                      <  2:780505/150‑1 (MQ=255)
                     ctgcggcGCACGGTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTc                                                                                                                                  <  1:788876/144‑1 (MQ=255)
                     cTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTc                                                                                                                                  >  1:1190843/1‑150 (MQ=255)
                      tGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCa                                                                                                                                 <  2:58788/150‑1 (MQ=255)
                       gAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCAt                                                                                                                                <  1:972953/150‑1 (MQ=255)
                        aGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATc                                                                                                                               <  1:329838/150‑1 (MQ=255)
                          tCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGa                                                                                                                             <  1:919894/150‑1 (MQ=255)
                           cGCAAGTAGCGGGATACAAACGAGAAATGCAGGCAAAACCATCCGATAACGCATACACGCCGCATCCGACATCTTTTGCCAGATGCTTCGTATTTACGCACCAAACACCATCGCCCCCGCGTGGGCAAAACGTTGCATCCCTTCATCGat                                                                                                                            >  2:1131480/1‑150 (MQ=255)
                             cAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGatat                                                                                                                          <  1:802888/150‑1 (MQ=255)
                             cAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGatat                                                                                                                          <  1:450800/150‑1 (MQ=255)
                             cAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGACTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGCAACGTTGCATCCCTTCATCGatat                                                                                                                          >  2:1046510/1‑150 (MQ=255)
                               aGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATcc                                                                                                                        <  1:40558/150‑1 (MQ=255)
                               aGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATcc                                                                                                                        >  1:556369/1‑150 (MQ=255)
                                       gATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCATATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTc                                                                                                                <  2:509038/150‑1 (MQ=255)
                                              aGGCGTAATTCAGGCTAAACCTTCGGATAATGCATTCACGCCGCATCCGACATCTTTTGCCTGATGGTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTccaccacc                                                                                                         <  2:368290/150‑1 (MQ=255)
                                                gCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAg                                                                                                       <  1:1247061/150‑1 (MQ=255)
                                                     aTGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGCCg                                                                                                  >  2:470493/1‑150 (MQ=255)
                                                       gCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGt                                                                                                >  2:277299/1‑150 (MQ=255)
                                                         aGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGc                                                                                              >  1:9567/1‑150 (MQ=255)
                                                            cTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaaa                                                                                           >  1:109160/1‑150 (MQ=255)
                                                              aaaCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCACCCCCCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATATTCCACCACCAGCGACGGTGCAAAAc                                                                                         >  1:1042231/1‑150 (MQ=255)
                                                                gccATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGCAGGTGCAAAACGc                                                                                       >  1:1082802/2‑150 (MQ=255)
                                                                 ccATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCa                                                                                      <  1:1152639/150‑1 (MQ=255)
                                                                 ccATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCa                                                                                      <  2:572476/150‑1 (MQ=255)
                                                                    tCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCCCCACCAGCGACGGTGCAAAACGCATca                                                                                   >  1:31689/1‑150 (MQ=255)
                                                                     cGGATAAGGCATTCACGCCGCATCCGACATCTTATGCCTGATGCTAATTATTTACTCCCCAACAACCTTCGCCACCGCGTGGGCGAAACGTTGCATACCTTCATCGATATCCGCATCTTCCACCACCACAGACCGTGCAAAACGCATccc                                                                                  >  1:647135/1‑148 (MQ=255)
                                                                        aTAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtc                                                                               <  2:340736/150‑1 (MQ=255)
                                                                         tAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCATATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtcc                                                                              <  1:508836/150‑1 (MQ=255)
                                                                         tAAGGCAATCACGCCGCATCCGACATCTTTAGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGAAGGTGCAAAACGCATCACAtcc                                                                              <  2:569118/150‑1 (MQ=255)
                                                                          aaGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccg                                                                             <  1:1285473/150‑1 (MQ=255)
                                                                          aaGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccg                                                                             <  1:396343/150‑1 (MQ=255)
                                                                           aGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccgg                                                                            >  1:607823/1‑150 (MQ=255)
                                                                           aGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccgg                                                                            <  1:804942/150‑1 (MQ=255)
                                                                           aGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccgg                                                                            >  2:971381/1‑150 (MQ=255)
                                                                                 tCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGc                                                                      <  1:468830/150‑1 (MQ=255)
                                                                                  cACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCa                                                                     >  1:522881/1‑150 (MQ=255)
                                                                                   aCGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCAt                                                                    <  2:9567/150‑1 (MQ=255)
                                                                                     gccgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTc                                                                  <  1:828403/150‑1 (MQ=255)
                                                                                       cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCat                                                                <  1:16898/150‑2 (MQ=255)
                                                                                       cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAg                                                                <  1:1038853/150‑1 (MQ=255)
                                                                                       cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACAGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCCCCACCAGCGACGGTGCAAAACGCCACACATCCAGTCCGGCATTCAg                                                                >  2:84883/1‑150 (MQ=255)
                                                                                            ccGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCa                                                           <  1:666009/150‑1 (MQ=255)
                                                                                             cGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCAt                                                          >  1:1107908/1‑150 (MQ=255)
                                                                                              gACATCTTTTGCCTGATGCTTCGTATTTACGCACCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGAAAAACGCATCAAATCCGGTCCGGCATTCAGCACCAtt                                                         >  1:1258229/1‑150 (MQ=255)
                                                                                                         ccTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctc                                              >  2:745782/1‑150 (MQ=255)
                                                                                                         ccTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctc                                              <  1:1074736/150‑1 (MQ=255)
                                                                                                         ccTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCcta                                              >  2:1260129/1‑149 (MQ=255)
                                                                                                           tGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTg                                            <  2:1094614/150‑1 (MQ=255)
                                                                                                             aTGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTgcg                                          <  2:1313506/150‑1 (MQ=255)
                                                                                                               gCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAACCGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGcc                                        >  1:587829/1‑150 (MQ=255)
                                                                                                                cTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCGACATTACGCCAGCCTCTGCGccc                                       >  2:1279391/1‑150 (MQ=255)
                                                                                                                 ttCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCg                                      <  1:285115/150‑1 (MQ=255)
                                                                                                                      aTTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCCTTACGCCAGCCTCTGCGCCCGCATAc                                 >  1:1050034/1‑150 (MQ=255)
                                                                                                                       tttACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACa                                >  1:665516/1‑150 (MQ=255)
                                                                                                                         tACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAgg                              >  1:293062/1‑150 (MQ=255)
                                                                                                                             ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCCCCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                          >  2:1019554/1‑150 (MQ=255)
                                                                                                                             ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCCCCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                          >  2:1271435/1‑150 (MQ=255)
                                                                                                                               ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                        <  1:959771/150‑1 (MQ=255)
                                                                                                                                cAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAc                       >  2:406393/1‑150 (MQ=255)
                                                                                                                                cAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCAACGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAc                       <  2:1176908/150‑1 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATTAcg                      <  2:572606/150‑1 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg                      <  2:109160/150‑1 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg                      <  2:556369/150‑1 (MQ=255)
                                                                                                                                    caccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACgcgc                   <  1:363869/150‑1 (MQ=255)
                                                                                                                                    acgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAGTCACgcgc                   >  1:1054731/5‑150 (MQ=255)
                                                                                                                                     cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAATTCACGCGcc                  >  2:975182/4‑150 (MQ=255)
                                                                                                                                         tCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCAGCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGAc              <  1:1260115/150‑1 (MQ=255)
                                                                                                                                            ccACCGCGTGGGCGAAACGTTGTATCCCTTCATCGAAATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCtt           <  2:660882/150‑1 (MQ=255)
                                                                                                                                            ccACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGTAATCACGCGCCCGAccgt           >  1:53634/1‑148 (MQ=255)
                                                                                                                                               ccGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGCCCTTTGt        >  2:1123617/1‑150 (MQ=255)
                                                                                                                                                  cgTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACt     <  1:229738/150‑1 (MQ=255)
                                                                                                                                                     gggCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACtgtg  <  1:969008/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACTGTG  >  NC_000913/3488840‑3489138

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