Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,351 C→T 9.3% R150H (CGC→CAC) ‡ hybE ← hydrogenase 2‑specific chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,3510CT 9.3% 253.2 / 11.3 86R150H (CGC→CAC) ‡hybEhydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base C (43/35);  new base T (3/5);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGA  >  NC_000913/3140203‑3140493
                                                                                                                                                    |                                                                                                                                              
aCGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                                                                                               <  1:823633/150‑1 (MQ=255)
 cGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGa                                                                                                                                              >  1:447329/1‑150 (MQ=255)
     ggAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                                                                                          <  1:933886/150‑1 (MQ=255)
      gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGccc                                                                                                                                         <  1:537958/150‑1 (MQ=255)
      gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGccc                                                                                                                                         <  1:763341/150‑1 (MQ=255)
        gagCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTg                                                                                                                                       <  1:1112443/150‑1 (MQ=255)
          gcgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGct                                                                                                                                     <  1:728761/150‑2 (MQ=255)
           cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgt                                                                                                                                    >  1:1302191/1‑150 (MQ=255)
            gcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgt                                                                                                                                   <  1:768012/150‑4 (MQ=255)
            gcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgt                                                                                                                                   <  1:550489/150‑4 (MQ=255)
             cgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgtt                                                                                                                                  <  1:1345641/150‑5 (MQ=255)
             cgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt                                                                                                                                  >  2:114416/1‑150 (MQ=255)
               ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCTTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGTCTACTCATGCATTTTCGCCACTCCTTCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                                <  2:1366376/150‑1 (MQ=255)
               ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgtttg                                                                                                                                <  1:1306493/150‑7 (MQ=255)
               ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                                >  1:1489923/1‑150 (MQ=255)
                cAATTTCCCGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCCCCGCCTGCCGCTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGt                                                                                                                               <  1:586008/150‑1 (MQ=255)
                cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGt                                                                                                                               >  1:765327/1‑150 (MQ=255)
                      ccAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATcc                                                                                                                         >  2:1126746/1‑150 (MQ=255)
                                     tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGc                                                                                                          <  2:754017/150‑1 (MQ=255)
                                     tGACGCGCTTAACAACCTGCTGCTCCGCCTGCCCTAGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGc                                                                                                          >  2:1436946/1‑150 (MQ=255)
                                      gACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGTCCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCa                                                                                                         <  2:1293668/150‑1 (MQ=255)
                                      gACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCa                                                                                                         >  1:100366/1‑150 (MQ=255)
                                         gcgcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCg                                                                                                      <  1:1154515/150‑1 (MQ=255)
                                         gcgcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCg                                                                                                      >  2:43388/1‑150 (MQ=255)
                                             ttAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAg                                                                                                  <  1:329955/150‑1 (MQ=255)
                                              tAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAgg                                                                                                 >  2:13046/1‑150 (MQ=255)
                                               aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGa                                                                                                >  2:589907/1‑150 (MQ=255)
                                               aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGa                                                                                                <  1:1023598/150‑1 (MQ=255)
                                                aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                               >  1:145091/1‑150 (MQ=255)
                                                 cATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATc                                                                                              <  2:177598/150‑1 (MQ=255)
                                                  aTCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGCCAACTCATGCATTTTCGCCACTCCTGCGCCCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCa                                                                                             >  1:206597/1‑150 (MQ=255)
                                                  aTCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCa                                                                                             <  1:958273/150‑1 (MQ=255)
                                                        ccgctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTg                                                                                       >  2:813151/3‑150 (MQ=255)
                                                         tgctCCGCCTGCCGTTGGGTAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATGCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgc                                                                                      >  2:299833/1‑150 (MQ=255)
                                                           ctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                                                                                    <  2:491823/150‑1 (MQ=255)
                                                             ccGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAg                                                                                  >  1:371962/1‑150 (MQ=255)
                                                                   gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATcc                                                                            >  2:11999/1‑150 (MQ=255)
                                                                   gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATcc                                                                            >  2:973017/1‑150 (MQ=255)
                                                                   gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATcc                                                                            <  2:1489923/150‑1 (MQ=255)
                                                                   gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATcc                                                                            >  2:861283/1‑150 (MQ=255)
                                                                    ccGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCg                                                                           >  1:1354763/1‑150 (MQ=255)
                                                                     cGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGt                                                                          <  1:232977/150‑1 (MQ=255)
                                                                     cGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGt                                                                          >  1:582764/1‑150 (MQ=255)
                                                                           aTAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCg                                                                    <  1:1436946/150‑1 (MQ=255)
                                                                            tAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGt                                                                   >  2:335817/1‑150 (MQ=255)
                                                                              aTTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTg                                                                 >  1:175388/1‑150 (MQ=255)
                                                                               tttCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTgg                                                                <  1:241280/150‑1 (MQ=255)
                                                                               tttCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGAACAAACAGCAGGGCGCGACGCCCAGCGTGTGGTACATCCGGATTCGTGCCTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTgg                                                                >  2:42380/1‑150 (MQ=255)
                                                                                  cAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGccc                                                             <  2:100366/150‑1 (MQ=255)
                                                                                   aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCct                                                            <  2:447329/150‑1 (MQ=255)
                                                                                     cGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCctct                                                          <  1:553918/150‑1 (MQ=255)
                                                                                      ggCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCtt                                                         <  1:502985/150‑1 (MQ=255)
                                                                                      ggCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCtt                                                         <  1:1371958/150‑1 (MQ=255)
                                                                                      ggCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCtt                                                         <  2:145091/150‑1 (MQ=255)
                                                                                       gcgcTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTc                                                        <  2:582764/150‑1 (MQ=255)
                                                                                       gcgcTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTc                                                        <  1:235943/150‑1 (MQ=255)
                                                                                           tctGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATc                                                    <  2:1302191/150‑1 (MQ=255)
                                                                                            ctGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCg                                                   >  2:204739/1‑150 (MQ=255)
                                                                                             tGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAAt                                                     <  2:1130908/147‑1 (MQ=255)
                                                                                             tGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAAt                                                     >  1:1130908/1‑147 (MQ=255)
                                                                                              ggCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGAc                                                 <  2:947819/150‑1 (MQ=255)
                                                                                               ggcaaGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACa                                                >  2:990638/4‑150 (MQ=255)
                                                                                                        cAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCg                                       <  1:335817/150‑1 (MQ=255)
                                                                                                           cTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGaaa                                    <  1:861283/150‑1 (MQ=255)
                                                                                                           cTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGaaa                                    >  2:257333/1‑150 (MQ=255)
                                                                                                             cATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgc                                  >  1:596819/1‑150 (MQ=255)
                                                                                                             cATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgc                                  >  1:370072/1‑150 (MQ=255)
                                                                                                              aTGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcg                                 >  1:449955/1‑150 (MQ=255)
                                                                                                               tGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcgg                                >  1:464617/1‑150 (MQ=255)
                                                                                                                 cATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGGAGTCATCCGTCATGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggag                              >  2:1407069/1‑148 (MQ=255)
                                                                                                                 cATTTTCGCAACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggcg                              >  2:1451641/1‑150 (MQ=255)
                                                                                                                  aTTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATACGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGa                             >  1:1087632/1‑150 (MQ=255)
                                                                                                                   ttttCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGAc                            <  1:973017/150‑1 (MQ=255)
                                                                                                                        gCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCCTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAg                       >  2:551128/1‑150 (MQ=255)
                                                                                                                         ccACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGc                      >  2:428052/1‑150 (MQ=255)
                                                                                                                          cACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                     <  1:834335/150‑1 (MQ=255)
                                                                                                                          cACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                     >  2:44930/1‑150 (MQ=255)
                                                                                                                              ccTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACa                 >  1:1366279/1‑150 (MQ=255)
                                                                                                                               cTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAg                >  1:652788/1‑150 (MQ=255)
                                                                                                                                      cAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGat         <  1:872869/150‑1 (MQ=255)
                                                                                                                                       aaaCAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGata        >  2:1242577/1‑150 (MQ=255)
                                                                                                                                          acgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTg     >  2:790066/3‑150 (MQ=255)
                                                                                                                                           cgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCCGTCCTCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAGAGCGGCGACATCAGCGAACAGGAGAGATATTGc    >  1:649724/2‑150 (MQ=255)
                                                                                                                                           cgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGc    >  2:1115162/2‑150 (MQ=255)
                                                                                                                                           agcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCACTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGc    <  2:1488151/150‑1 (MQ=255)
                                                                                                                                            gcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGCCATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCg   >  1:1388300/1‑150 (MQ=255)
                                                                                                                                             cagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGa  >  2:1521559/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGA  >  NC_000913/3140203‑3140493

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: