Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,836,188 G→A 5.2% A413V (GCG→GTG)  hypF ← carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,836,1880GA 5.2% 244.6 / 2.2 77A413V (GCG→GTG) hypFcarbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base A (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCTTTCGCG  >  NC_000913/2836040‑2836333
                                                                                                                                                    |                                                                                                                                                 
ccACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCg                                                                                                                                                   >  6:398829/1‑149 (MQ=255)
         ggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCAc                                                                                                                                          <  5:398829/149‑1 (MQ=255)
         ggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCAc                                                                                                                                          >  5:172015/1‑149 (MQ=255)
         ggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCAc                                                                                                                                          >  7:14503/1‑149 (MQ=255)
               aTAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCa                                                                                                                                    >  2:356790/1‑149 (MQ=255)
               aTAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCa                                                                                                                                    >  3:334354/1‑149 (MQ=255)
                tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCAc                                                                                                                                   >  7:331869/1‑149 (MQ=255)
                   cccGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAg                                                                                                                                <  7:160817/149‑1 (MQ=255)
                     cGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAgg                                                                                                                               >  6:61106/1‑148 (MQ=255)
                           cgcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAg                                                                                                                         <  2:210959/148‑1 (MQ=255)
                           cgcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAg                                                                                                                         <  5:33340/148‑1 (MQ=255)
                                  tGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTAttttt                                                                                                                 <  1:356790/149‑1 (MQ=255)
                                  tGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTAttttt                                                                                                                 <  8:14503/149‑1 (MQ=255)
                                  tGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTAttttt                                                                                                                 <  6:172015/149‑1 (MQ=255)
                                   gCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGtt                                                                                                                                             <  3:4292/120‑1 (MQ=255)
                                   gCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGtt                                                                                                                                             >  4:4292/1‑120 (MQ=255)
                                    cacaACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTc                                                                                                                <  4:108746/148‑1 (MQ=255)
                                    cacaACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCa                                                                                                               <  3:239963/149‑1 (MQ=255)
                                     acaACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAg                                                                                                              <  4:248600/149‑1 (MQ=255)
                                       aaCGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTAttttt                                                                                                                 >  2:332742/1‑144 (MQ=255)
                                         cGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCCCCAGGCAGAAGGTATTTTTCAGATc                                                                                                           >  5:406489/1‑148 (MQ=255)
                                         cGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATc                                                                                                           <  6:406489/148‑1 (MQ=255)
                                           tATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCcgc                                                                                                        <  5:61106/149‑1 (MQ=255)
                                                     tAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGaca                                                                                               >  2:201188/1‑148 (MQ=255)
                                                       aaaTCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGAcacag                                                                                            >  1:17391/1‑149 (MQ=255)
                                                       aaaTCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGAcaca                                                                                             <  8:331869/148‑1 (MQ=255)
                                                         aTCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGAcacagca                                                                                          <  2:17391/149‑1 (MQ=255)
                                                             tAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGAcacagcacag                                                                                       <  7:385773/148‑1 (MQ=255)
                                                             tAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGAcacagcacag                                                                                       >  8:385773/1‑148 (MQ=255)
                                                             tAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAgg                                                                                      >  6:369143/1‑149 (MQ=255)
                                                              aGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGAGGCTGACTCAACACCGCTTGTTCACACCGCCCCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGc                                                                                     >  3:302964/1‑149 (MQ=255)
                                                                        cATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCga                                                                                                   <  2:210848/125‑1 (MQ=255)
                                                                        cATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCga                                                                                                   >  1:210848/1‑125 (MQ=255)
                                                                        cATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAtttt                                                                           <  1:201188/149‑1 (MQ=255)
                                                                        cATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAttt                                                                            <  1:67569/148‑1 (MQ=255)
                                                                         aTCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAtttt                                                                           >  1:209205/1‑148 (MQ=255)
                                                                          tCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCgcgc                                                                                                                                     <  3:52767/89‑1 (MQ=255)
                                                                          tCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCgcgc                                                                                                                                     >  4:52767/1‑89 (MQ=255)
                                                                          tCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAttttt                                                                          >  1:241977/1‑148 (MQ=255)
                                                                          tCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAttttt                                                                          >  7:429597/1‑148 (MQ=255)
                                                                          tCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTa                                                                         <  2:209205/149‑1 (MQ=255)
                                                                            aGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTaaa                                                                       <  5:369143/149‑1 (MQ=255)
                                                                                 tAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCg                                                                  >  1:111626/1‑149 (MQ=255)
                                                                                 tAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCACTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTaaa                                                                       <  6:232545/144‑1 (MQ=255)
                                                                                 tAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCACTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTaaa                                                                       >  5:105316/1‑144 (MQ=255)
                                                                                 tAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCACTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTaaa                                                                       >  5:232545/1‑144 (MQ=255)
                                                                                 tAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCACTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTaaa                                                                       <  6:105316/144‑1 (MQ=255)
                                                                                    cGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAg                                                               <  8:323549/149‑1 (MQ=255)
                                                                                          gcgcCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCaaa                                                          <  7:8946/148‑1 (MQ=255)
                                                                                          gcgcCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAg                                                         >  4:50139/1‑149 (MQ=255)
                                                                                           cgcCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGc                                                        <  3:262899/149‑1 (MQ=255)
                                                                                             ccACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCa                                                      <  2:241977/149‑1 (MQ=255)
                                                                                                cTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACCCCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAgcg                                                   >  1:221539/1‑149 (MQ=255)
                                                                                                 tGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAgcg                                                   >  1:137670/1‑148 (MQ=255)
                                                                                                       tGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCgg                                            >  4:169328/1‑149 (MQ=255)
                                                                                                        ggATTCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCTCACCAGGCa                                                                                                                             <  5:381554/67‑1 (MQ=255)
                                                                                                        ggATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCTCACCAGGCa                                                                                                                             >  6:381554/1‑67 (MQ=255)
                                                                                                         gATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCgg                                            <  3:169328/147‑1 (MQ=255)
                                                                                                         gATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGca                                          >  1:251059/1‑149 (MQ=255)
                                                                                                           tGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGcaca                                        <  2:230337/149‑1 (MQ=255)
                                                                                                              cATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATAc                                     >  7:370069/1‑149 (MQ=255)
                                                                                                                tCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGATCACCAGGCAGAAGGTATTTTTCAAATCCGCGCCGAGAcaca                                                                                             >  5:262157/1‑91 (MQ=255)
                                                                                                                     aCTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCTAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATAccccccc                               >  6:410926/1‑146 (MQ=255)
                                                                                                                       ttAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGaa                            >  7:294392/1‑149 (MQ=255)
                                                                                                                       ttAAATCGCCCAGATGCTGACTCAACACCGCTTGTTAACCGCGCACCAGCCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTATAAAGCCCGGAGGCAAAGCCAGCGCATCAGGCACAAACCCCCGCGaa                            >  2:149542/1‑149 (MQ=255)
                                                                                                                          aaTCGCCCAGATGCTGACTCACCACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACgg                         <  4:334354/149‑1 (MQ=255)
                                                                                                                                ccAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCgcag                    >  4:374392/1‑148 (MQ=255)
                                                                                                                                ccAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCgcag                    >  2:187905/1‑148 (MQ=255)
                                                                                                                                  aGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCAt                 >  1:409213/1‑149 (MQ=255)
                                                                                                                                           cTCAACACCGCTTGTTCACCGCGCACCAGGCAGAATGTATTTTTCAGCTCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCt        <  2:221539/149‑1 (MQ=255)
                                                                                                                                            taaaCACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCtt       >  2:27153/3‑149 (MQ=255)
                                                                                                                                            tCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCgcagca                  >  4:114962/1‑138 (MQ=255)
                                                                                                                                            tCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCgcagca                  <  3:114962/138‑1 (MQ=255)
                                                                                                                                            tCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCtt       >  4:72484/1‑149 (MQ=255)
                                                                                                                                            tCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCtt       >  5:144422/1‑149 (MQ=255)
                                                                                                                                              aaCACCGCTTGTTCACAGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCttt      <  8:85021/148‑1 (MQ=255)
                                                                                                                                                cacCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCTTTcgc   >  4:227305/1‑149 (MQ=255)
                                                                                                                                                  ccGCTTGTTCACCGCGCCCCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGAGCATCCGGCTCATACCCCCGCGAACGGCGCAGCATTTCGCCGCTTTcgcg  >  5:270977/1‑148 (MQ=255)
                                                                                                                                                    gCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAt                                                                              >  6:24546/1‑70 (MQ=255)
                                                                                                                                                    gCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACAt                                                                              <  5:24546/70‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCGCGCCACTGCATCTGGATGCCATCGTCACTTAAATCGCCCAGATGCTGACTCAACACCGCTTGTTCACCGCGCACCAGGCAGAAGGTATTTTTCAGATCCGCGCCGAGACACAGCACAGGCGGAACATTTTTAAAGCCCGGAGGCAAAGCCAGCGCATCCGGCACATACCCCCGCGAACGGCGCAGCATTTCGCCGCTTTCGCG  >  NC_000913/2836040‑2836333

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