Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,852 C→T 5.2% R236H (CGT→CAT)  btuD ← vitamin B12 ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,8520CT 5.2% 248.3 / 2.4 77R236H (CGT→CAT) btuDvitamin B12 ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base T (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCAT  >  NC_000913/1792705‑1792997
                                                                                                                                                   |                                                                                                                                                 
tgtgATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACg                                                                                                                                                  >  6:167032/1‑149 (MQ=255)
    aTCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCg                                                                                                                                               >  4:28672/1‑148 (MQ=255)
     tCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGa                                                                                                                                              >  2:398959/1‑148 (MQ=255)
       aaaaTCAGGCAGAATCGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGaaaa                                                                                                                                           <  2:385225/149‑1 (MQ=255)
        aaaTAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAAt                                                                                                                                          >  1:10858/1‑149 (MQ=255)
         aaTAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAAtt                                                                                                                                         <  4:132562/149‑1 (MQ=255)
          aTAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTc                                                                                                                                        <  5:167032/149‑1 (MQ=255)
             aGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCAt                                                                                                                                      <  7:197595/148‑1 (MQ=255)
                   aaaCGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAg                                                                                                                               >  8:202150/1‑149 (MQ=255)
                        cATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGCTGGCGAAAATTCATCCCATAGGCCTg                                                                                                                          >  6:308523/1‑149 (MQ=255)
                         aTTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTgc                                                                                                                         <  7:7973/149‑1 (MQ=255)
                         aTTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTgc                                                                                                                         <  1:398959/149‑1 (MQ=255)
                          ttCCGTAGTTTTCTTCTGTTTTTTATTCTTTCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATTCTAATGAAATCAGATGGTCGAAATCATCATTCTGTGACCTTCGATATCCATACTGCGAAAATTCATCCCATAGGCCTgc                                                                                                                         <  2:270602/148‑1 (MQ=255)
                             cGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCa                                                                                                                     <  8:279338/149‑1 (MQ=255)
                                 cGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGAtt                                                                                                                 >  3:327918/1‑149 (MQ=255)
                                            ttttttATTCTTGCATTTATTTAGCGGCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGt                                                                                                        <  3:303884/147‑1 (MQ=255)
                                            ttttttATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTg                                                                                                       <  3:286942/148‑1 (MQ=255)
                                                 tATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGCCCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCAc                                                                                                  <  8:36440/148‑1 (MQ=255)
                                                 tATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCAc                                                                                                  <  3:28672/148‑1 (MQ=255)
                                                   ttCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACctc                                                                                               <  6:370880/149‑1 (MQ=255)
                                                      ttGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTc                                                                                            >  4:359825/1‑149 (MQ=255)
                                                        gCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTccc                                                                                          >  2:119143/1‑149 (MQ=255)
                                                          aTTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTg                                                                                        >  4:226951/1‑149 (MQ=255)
                                                           ttAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGc                                                                                       >  3:240599/1‑149 (MQ=255)
                                                             aaTTTAGCGTCGTAATTCCCCGATTATCAAGATACTAATGAAATCAGATGGTCGAAAGCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATACGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCg                                                                                      <  1:386124/148‑1 (MQ=255)
                                                              aTTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGt                                                                                     <  6:137151/148‑1 (MQ=255)
                                                              aTTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTc                                                                                    >  3:165974/1‑149 (MQ=255)
                                                               tttAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTc                                                                                    >  6:100212/1‑148 (MQ=255)
                                                                   gCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTg                                                                               <  4:240599/149‑1 (MQ=255)
                                                                    cgtcgtAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTg                                                                               <  4:327918/148‑1 (MQ=255)
                                                                     gtcgtAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGc                                                                             >  1:166221/1‑149 (MQ=255)
                                                                        gtAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCa                                                                           <  5:277638/148‑1 (MQ=255)
                                                                            ttACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGATGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCAtt                                                                      >  6:228063/1‑149 (MQ=255)
                                                                              aCCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCAtttt                                                                    >  6:76699/1‑149 (MQ=255)
                                                                                        cAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCCTCCACTGGCCAGCATTTTTCCACCttt                                                           >  8:231895/1‑148 (MQ=255)
                                                                                         aaGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCttt                                                           >  6:269252/1‑147 (MQ=255)
                                                                                         aaGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTa                                                         <  2:23532/149‑1 (MQ=255)
                                                                                         aaGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTa                                                         <  5:269252/149‑1 (MQ=255)
                                                                                                     aaaTCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGccc                                             >  4:72528/1‑149 (MQ=255)
                                                                                                          aGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATgc                                        <  6:86385/149‑1 (MQ=255)
                                                                                                           gATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATgcg                                       <  5:100212/149‑1 (MQ=255)
                                                                                                            aTGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCttt                                                           <  6:276776/128‑1 (MQ=255)
                                                                                                            aTGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCttt                                                           >  5:276776/1‑128 (MQ=255)
                                                                                                             tGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCa                                     <  3:226951/149‑1 (MQ=255)
                                                                                                             tGGTCGAAAGCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATCGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCa                                     <  4:348099/149‑1 (MQ=255)
                                                                                                              ggTCGTATTCTGCATTCTGTGTCCTTTGATTTGCAGAGTTGGTAAATTCATCCCTTTGGCCTGCGTCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCa                                     <  8:386233/148‑1 (MQ=255)
                                                                                                              ggTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCAt                                    >  2:176300/1‑149 (MQ=255)
                                                                                                               gTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCA                                                      <  6:320276‑M1/130‑1 (MQ=255)
                                                                                                               gTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCA                                                      >  5:320276‑M1/1‑130 (MQ=255)
                                                                                                                tCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGa                                  >  8:261054/1‑149 (MQ=255)
                                                                                                                  gAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACg                                <  3:72528/149‑1 (MQ=255)
                                                                                                                  gAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACg                                >  2:238143/1‑149 (MQ=255)
                                                                                                                      tCAGCATTCTGTGACCTTCGATATCCAGATGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAAt                            <  5:308523/149‑1 (MQ=255)
                                                                                                                           aTTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTgtg                       >  1:119143‑M2/1‑146 (MQ=255)
                                                                                                                             tCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTc                                                                                    >  6:10268/1‑86 (MQ=255)
                                                                                                                             tCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTc                                                                                    <  5:10268/86‑1 (MQ=255)
                                                                                                                               tgtgACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCAt                                                                       <  5:78217/97‑1 (MQ=255)
                                                                                                                               tgtgACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCAt                                                                       >  6:78217/1‑97 (MQ=255)
                                                                                                                               tgtgACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTgtggtt                    <  1:205760‑M2/148‑7 (MQ=255)
                                                                                                                                 tgACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGTGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACacattgcgtcatgcgcatcgggcgtggttgc                  <  5:238554‑M1/148‑32 (MQ=255)
                                                                                                                                 tgACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTgag                 <  1:238143/149‑1 (MQ=255)
                                                                                                                                    ccTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTgtggtt                    <  8:23485‑M2/143‑7 (MQ=255)
                                                                                                                                      ttCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            >  5:118250/1‑149 (MQ=255)
                                                                                                                                      ttCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            >  7:101171/1‑149 (MQ=255)
                                                                                                                                        cGATATCCAGATGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            <  5:228063/147‑1 (MQ=255)
                                                                                                                                        cGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            <  6:118250/147‑1 (MQ=255)
                                                                                                                                        cGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            <  2:166221/147‑1 (MQ=255)
                                                                                                                                         gATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTACCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGAc         >  2:130000/1‑149 (MQ=255)
                                                                                                                                         gATATCCAGACGGCGAAAATGCATCCCATAGGCCTGCGCCAGATTTGGCGGCTTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGa          <  3:274149/148‑1 (MQ=255)
                                                                                                                                          atatCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGGGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACt        >  3:240720/1‑149 (MQ=255)
                                                                                                                                              ccAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            <  5:37149/141‑1 (MQ=255)
                                                                                                                                              ccAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGt            >  6:37149/1‑141 (MQ=255)
                                                                                                                                               cAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCtttt                                                          <  3:188522/94‑1 (MQ=255)
                                                                                                                                               cAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCtttt                                                          >  4:188522/1‑94 (MQ=255)
                                                                                                                                               cAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCtc    <  7:231895/148‑1 (MQ=255)
                                                                                                                                                aGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCtcat  >  7:7847/1‑149 (MQ=255)
                                                                                                                                                aGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCtcat  >  3:346048/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGTGATCAAAATAAGGCAGAAACGCATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCAT  >  NC_000913/1792705‑1792997

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