Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,387,553 T→A 5.1% R212R (CGA→CGT aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,387,5530TA 5.1% 288.8 / 5.0 78R212R (CGA→CGTaaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base A (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GCGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATCCG  >  NC_000913/3387405‑3387701
                                                                                                                                                    |                                                                                                                                                    
gcgTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGt                                                                                                                                                      <  4:520938/149‑1 (MQ=255)
               gataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGccc                                                                                                                                       >  4:869552/1‑149 (MQ=255)
                  aagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAg                                                                                                                                    >  4:1135284/1‑149 (MQ=255)
                         gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                             <  2:1013998/149‑1 (MQ=255)
                          cAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                            <  4:269449/149‑1 (MQ=255)
                             gATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                              <  1:473411/128‑1 (MQ=255)
                             gATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                              >  2:473411/1‑128 (MQ=255)
                                  ggTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACtt                                                                                                                    >  3:1397262/1‑149 (MQ=255)
                                     aaTCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCAc                                                                                                                 <  3:1135284/149‑1 (MQ=255)
                                       tcagcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCa                                                                                                               >  3:810298/1‑149 (MQ=255)
                                              cagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCAt                                                                                                        <  3:869552/149‑1 (MQ=255)
                                                    tAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTCCACCATCGCCATGCTTGa                                                                                                  >  1:391591/1‑149 (MQ=255)
                                                       cGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATAc                                                                                               <  2:391591/149‑1 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                              <  1:343693/101‑1 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                              >  2:343693/1‑101 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGa                                                                                                  <  3:1009850/145‑1 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATAATCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGa                                                                                                  >  4:1009850/1‑145 (MQ=255)
                                                           ttGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGc                                                                                                                                                   <  1:168440/93‑1 (MQ=255)
                                                           ttGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGc                                                                                                                                                   >  2:168440/1‑93 (MQ=255)
                                                              aTCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                        >  1:1208490/1‑149 (MQ=255)
                                                              aTCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                        >  1:876474/1‑149 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCGGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                               >  3:266244/1‑93 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCGGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                               <  4:266244/93‑1 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACtt                                                                                                                    <  2:246966/120‑1 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACtt                                                                                                                    >  1:246966/1‑120 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGt                                                                                          >  4:513492/1‑146 (MQ=255)
                                                               tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGt                                                                                          <  3:513492/146‑1 (MQ=255)
                                                                        aaCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGa                                                                              <  4:810298/149‑1 (MQ=255)
                                                                        aaCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGa                                                                              <  2:1208490/149‑1 (MQ=255)
                                                                           gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCa                                                                                                               <  1:975209/113‑1 (MQ=255)
                                                                           gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCa                                                                                                               >  2:975209/1‑113 (MQ=255)
                                                                               aTTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcg                                                                       >  2:419560/1‑149 (MQ=255)
                                                                                     ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                     >  3:452282/1‑113 (MQ=255)
                                                                                     ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                                     <  4:452282/113‑1 (MQ=255)
                                                                                      cgcgCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAg                                                                >  1:89206/1‑149 (MQ=255)
                                                                                          cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGATAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            <  4:1397262/149‑1 (MQ=255)
                                                                                          cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            <  2:876474/149‑1 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                             >  1:597505/1‑83 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                             <  2:597505/83‑1 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAtt                                                                                    >  2:1152528/1‑124 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAtt                                                                                    <  1:1152528/124‑1 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                                                           >  3:483380/1‑149 (MQ=255)
                                                                                               cGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGCGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       <  2:830980/149‑1 (MQ=255)
                                                                                                  gAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAg                                                    >  4:537122/1‑149 (MQ=255)
                                                                                                     gATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       >  1:679220/1‑143 (MQ=255)
                                                                                                     gATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       <  2:679220/143‑1 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                        <  1:1134691/124‑1 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                        >  2:1134691/1‑124 (MQ=255)
                                                                                                        tCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCg                                              >  1:1054064/1‑149 (MQ=255)
                                                                                                         ccATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGa                                             <  1:419560/149‑1 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGa                                                                              <  4:88223/115‑1 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGa                                                                              >  3:88223/1‑115 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt                                            >  1:1384455/1‑149 (MQ=255)
                                                                                                           atatTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       >  1:639288/1‑137 (MQ=255)
                                                                                                           atatTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       <  2:639288/137‑1 (MQ=255)
                                                                                                              ttCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGGTTGATACa                                                                                              <  4:123151/95‑1 (MQ=255)
                                                                                                              ttCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACa                                                                                              >  3:123151/1‑95 (MQ=255)
                                                                                                              ttCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       >  4:1219441/1‑134 (MQ=255)
                                                                                                              ttCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                                                       <  3:1219441/134‑1 (MQ=255)
                                                                                                                   gTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtctt                                   <  2:1384455/149‑1 (MQ=255)
                                                                                                                      gctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGtt                                <  4:483380/149‑1 (MQ=255)
                                                                                                                          cgcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgat                            <  2:89206/149‑1 (MQ=255)
                                                                                                                           gcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatc                           <  3:537122/149‑1 (MQ=255)
                                                                                                                              tGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGAt                                                                                                 <  1:1145027/76‑1 (MQ=255)
                                                                                                                              tGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGAt                                                                                                 >  2:1145027/1‑76 (MQ=255)
                                                                                                                               gCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatcgatc                       >  3:541559/1‑149 (MQ=255)
                                                                                                                                  ttGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgcg                    >  1:1361643/1‑149 (MQ=255)
                                                                                                                                   tGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                                                <  1:18846/88‑1 (MQ=255)
                                                                                                                                   tGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                                                >  2:18846/1‑88 (MQ=255)
                                                                                                                                      aGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgcg                    <  1:56436/145‑1 (MQ=255)
                                                                                                                                      aGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgcg                    >  2:56436/1‑145 (MQ=255)
                                                                                                                                           cGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga           <  2:1054064/149‑1 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            <  4:1058748/98‑1 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            <  4:580426/98‑1 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            >  3:580426/2‑99 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                            >  3:1058748/2‑99 (MQ=255)
                                                                                                                                             tggtgCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGc         <  2:1361643/149‑1 (MQ=255)
                                                                                                                                              ggtgCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCa        <  4:1103167/149‑1 (MQ=255)
                                                                                                                                              ggtgCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCa        >  3:1103167/1‑149 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                                                <  3:559037/71‑1 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                                                >  4:559037/1‑71 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTg                                  <  4:137627/117‑1 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTg                                  >  3:137627/1‑117 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga           >  4:1117855/1‑140 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga           <  3:1117855/140‑1 (MQ=255)
                                                                                                                                                    tCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATCCg  >  3:446378/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAGCAAATCCG  >  NC_000913/3387405‑3387701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: