Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,187,352 G→A L297L (CTG→TTG)  roxA ← 50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,3520GA100.0% 235.5 / NA 61L297L (CTG→TTG) roxA50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/34);  total (27/34)

GCAGTCAGCGCAATGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAGCGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCAGTTTAT  >  NC_000913/1187210‑1187498
                                                                                                                                              |                                                                                                                                                  
gCAGTCAGCGCAATGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGt                                                                                                                                            <  1:822793/151‑1 (MQ=255)
        cgcAATGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCGTCGTAGATTTc                                                                                                                                     >  2:1036108/1‑150 (MQ=255)
        cgcAATGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCa                                                                                                                                    <  1:768949/151‑1 (MQ=255)
            aTGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCg                                                                                                                                >  2:892330/1‑151 (MQ=255)
            aTGTTGCTGGCGAGTGCATCCAGTACCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCg                                                                                                                                >  2:156416/1‑151 (MQ=255)
                 gCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTg                                                                                                                            <  1:527826/150‑1 (MQ=255)
                    ggCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATaa                                                                                                                        <  1:378316/151‑1 (MQ=255)
                       gAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAggc                                                                                                                     >  1:173735/1‑151 (MQ=255)
                          tGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAggcggc                                                                                                                  >  2:89318/1‑151 (MQ=255)
                           gCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCCCCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGct                                                                                                                 <  1:221204/151‑1 (MQ=255)
                             aTCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTTGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGctct                                                                                                               <  1:449435/151‑1 (MQ=255)
                                  gTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTggcgg                                                                                                          >  2:319574/1‑151 (MQ=255)
                                   tGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTggcggc                                                                                                         >  2:998691/1‑151 (MQ=255)
                                    gCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGgcggcg                                                                                                        <  1:272197/151‑1 (MQ=255)
                                     ccGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGgcgc                                                                                                       >  1:382060/1‑151 (MQ=255)
                                      cGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACTCGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCCCTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGcgcg                                                                                                      <  2:197642/151‑1 (MQ=255)
                                      cGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGcgcg                                                                                                      <  1:26263/151‑1 (MQ=255)
                                      cGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGcgcg                                                                                                      >  2:706205/1‑151 (MQ=255)
                                           ggTGCGGGGAATCGATCTTCTCACCATTGGCATATACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTGAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATc                                                                                                 >  2:225417/1‑151 (MQ=255)
                                           ggTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATc                                                                                                 >  1:1033433/1‑151 (MQ=255)
                                              gCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATcc                                                                                                >  1:437321/1‑149 (MQ=255)
                                                    aaTCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATg                                                                                        <  2:38991/151‑1 (MQ=255)
                                                     aTCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGa                                                                                       <  2:867616/151‑1 (MQ=255)
                                                            ttctCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACt                                                                                >  1:289502/1‑151 (MQ=255)
                                                             tctcACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTg                                                                               >  2:424376/1‑151 (MQ=255)
                                                             tctcACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTg                                                                               >  2:772120/1‑151 (MQ=255)
                                                              ctcACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTgg                                                                              >  1:796604/1‑151 (MQ=255)
                                                              ctcACCATTCGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCACTCGCACCATCACTTCACCTTGTTTCAACGCATCGTATATTTCATCCGGCTGATAATTCCGCTCTTCCTGCTCTATATCCATTTCATTACTTCACt                                                                                <  1:307451/149‑1 (MQ=255)
                                                                cACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTg                                                                               <  2:932037/148‑1 (MQ=255)
                                                                   cATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGatat                                                                         >  2:735860/1‑151 (MQ=255)
                                                                   cATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGatat                                                                         <  1:453241/151‑1 (MQ=255)
                                                                         cATACACGTCGTCGCCAATGCGCTATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTc                                                                   <  2:72421/151‑1 (MQ=255)
                                                                         cATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTc                                                                   <  1:892181/151‑1 (MQ=255)
                                                                                gtcgtcGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAAc                                                            <  1:207073/151‑1 (MQ=255)
                                                                                  cgtcgCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAAcc                                                           >  2:398857/1‑150 (MQ=255)
                                                                                      gCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGc                                                      >  1:930162/1‑151 (MQ=255)
                                                                                       ccAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCt                                                     >  1:875207/1‑151 (MQ=255)
                                                                                       ccAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCt                                                     >  1:641089/1‑151 (MQ=255)
                                                                                                   aCGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTcc                                          <  2:1047407/150‑1 (MQ=255)
                                                                                                   aCGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCg                                         <  2:1033603/151‑1 (MQ=255)
                                                                                                        cAGACCACCCAGGCGCACCAGCCCTTGACCTTGTTTTAACGCTTCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCATTTCGTGCCGTGACTGGTATAGAAACTCGCCAAACCAGTGCTTGAAGTGTTTCGGCTgg                                    <  1:567650/151‑1 (MQ=255)
                                                                                                         aGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGt                                   <  2:360013/151‑1 (MQ=255)
                                                                                                           accaccCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTg                                 <  2:850009/151‑1 (MQ=255)
                                                                                                             caccCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGAt                               >  2:147474/1‑151 (MQ=255)
                                                                                                              accCAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATc                              <  2:131756/151‑1 (MQ=255)
                                                                                                                ccAGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCa                             >  1:1040953/1‑150 (MQ=255)
                                                                                                                  aGGCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAAtt                          >  1:378230/1‑151 (MQ=255)
                                                                                                                   ggCGCACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTc                         <  1:771992/151‑1 (MQ=255)
                                                                                                                     cgcACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGa                       <  2:200214/151‑1 (MQ=255)
                                                                                                                     cgcACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGa                       <  1:115645/151‑1 (MQ=255)
                                                                                                                      gcACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAg                      <  1:322694/151‑1 (MQ=255)
                                                                                                                       cACCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGc                     <  2:177967/151‑1 (MQ=255)
                                                                                                                        aCCAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGc                     >  1:382777/1‑150 (MQ=255)
                                                                                                                         ccAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGcat                   >  2:810405/1‑151 (MQ=255)
                                                                                                                          cAGCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGcatc                  <  2:34387/151‑1 (MQ=255)
                                                                                                                            gCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGGGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGcatcat                <  2:129621/151‑1 (MQ=255)
                                                                                                                            gCACTTCACCTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGcatcat                <  2:571731/151‑1 (MQ=255)
                                                                                                                                    ccTTGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCa        <  1:571685/151‑1 (MQ=255)
                                                                                                                                      ttGTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCAGt      >  2:835041/1‑151 (MQ=255)
                                                                                                                                        gTTTCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCAGttt    <  2:450227/151‑1 (MQ=255)
                                                                                                                                           tCAACGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCAGTTTAt  <  1:38519/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GCAGTCAGCGCAATGTTGCTGGCGAGTGCATCCAGTGCCGGACGGTGCGGGGAATCGATCTTCTCACCATTGGCATACACGTCGTCGCCAATGCGCAATACGCGCAGACCACCCAGGCGCACCAGCACTTCACCTTGTTTCAGCGCATCGTAGATTTCATCCGGCTGATAAGGCGGCTCTGGCGGCGCGATATCCAGTTCATGACGTGACTGGGATATAAACTCGCCAAACCATTGCTTAAAGTGTTCCGGCTGGTTGATCAATTCGAGCATCATCTCACGCAGTTTAT  >  NC_000913/1187210‑1187498

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