Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,545,565 G→A *228* (TAG→TAA gntX → DNA catabolic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,545,5650GA100.0% 291.1 / NA 89*228* (TAG→TAAgntXDNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (44/45);  total (44/45)

CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACA  >  NC_000913/3545431‑3545692
                                                                                                                                      |                                                                                                                               
cccGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                             >  1:835309/1‑139 (MQ=255)
cccGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                             >  2:753055/1‑139 (MQ=255)
    tGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCga                                                                                                                         >  1:987578/1‑139 (MQ=255)
       aaGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                      <  2:944909/139‑1 (MQ=255)
           tCGCCATATGGTGATTGTGGATGATGTCGTTACGACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  <  2:190729/139‑1 (MQ=255)
            cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 <  1:753055/139‑1 (MQ=255)
                       gATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                      <  1:700426/139‑1 (MQ=255)
                         ttGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgat                                                                                                                        >  2:289572/1‑119 (MQ=255)
                         ttGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgat                                                                                                                        <  1:289572/119‑1 (MQ=255)
                         ttGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCaa                                                                                                    >  1:369225/1‑139 (MQ=255)
                          tgtgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  <  1:80814/124‑1 (MQ=255)
                          tgtgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  >  2:80814/1‑124 (MQ=255)
                              gatgatGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 <  2:363992/121‑1 (MQ=255)
                              gatgatGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 >  1:363992/1‑121 (MQ=255)
                                    gTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         >  1:1141304/1‑139 (MQ=255)
                                          accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   <  2:835309/139‑1 (MQ=255)
                                             accGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                              >  2:647470/1‑125 (MQ=255)
                                             accGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                              <  1:647470/125‑1 (MQ=255)
                                              ccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGTTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAgg                                                                               <  2:867823/139‑1 (MQ=255)
                                              ccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGCCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                         <  1:540942/113‑1 (MQ=255)
                                              ccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGCCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                         >  2:540942/1‑113 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 <  2:574171/104‑1 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 >  1:574171/1‑104 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                              <  2:446306/139‑1 (MQ=255)
                                                 gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         >  2:579166/1‑126 (MQ=255)
                                                 gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         <  1:579166/126‑1 (MQ=255)
                                                    gTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTCGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   >  1:283061/1‑129 (MQ=255)
                                                    gTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTCGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   <  2:283061/129‑1 (MQ=255)
                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  <  2:1328665/93‑1 (MQ=255)
                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  >  1:1328665/1‑93 (MQ=255)
                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                              >  1:1137036/1‑129 (MQ=255)
                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                              <  2:1137036/129‑1 (MQ=255)
                                                              agagATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                             <  2:133025/77‑1 (MQ=255)
                                                              agagATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                             >  1:133025/1‑77 (MQ=255)
                                                                  aTTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   <  2:878187/115‑1 (MQ=255)
                                                                  aTTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   >  1:878187/1‑115 (MQ=255)
                                                                    tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                      <  1:1246057/78‑1 (MQ=255)
                                                                    tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                      >  2:1246057/1‑78 (MQ=255)
                                                                      cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCg                                                       >  2:849523/1‑139 (MQ=255)
                                                                           ttGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                            <  2:719558/113‑1 (MQ=255)
                                                                           ttGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                            >  1:719558/1‑113 (MQ=255)
                                                                            tGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgctgc                                                >  2:955760/1‑139 (MQ=255)
                                                                             gCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         <  2:347009/98‑1 (MQ=255)
                                                                             gCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         >  1:347009/1‑98 (MQ=255)
                                                                              cTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGca                                               <  2:369225/139‑1 (MQ=255)
                                                                                   aCGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAgc                                          <  2:1006026/139‑1 (MQ=255)
                                                                                   aCGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAgcg                                         <  1:955760/139‑1 (MQ=255)
                                                                                     gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                     <  1:1240962/62‑1 (MQ=255)
                                                                                     gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                     >  2:1240962/1‑62 (MQ=255)
                                                                                     gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                      <  1:498940/77‑1 (MQ=255)
                                                                                     gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                      >  2:498940/1‑77 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                            <  1:148959/99‑1 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                            >  2:148959/1‑99 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACtt                                    >  2:518257/1‑139 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACtt                                    >  1:459381/1‑139 (MQ=255)
                                                                                                 cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                      <  2:248683/65‑1 (MQ=255)
                                                                                                 cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                      >  1:248683/1‑65 (MQ=255)
                                                                                                 cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAAc                            >  1:884391/1‑139 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                              <  1:417194/71‑1 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                              >  2:417194/1‑71 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         >  2:300329/1‑76 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         <  1:300329/76‑1 (MQ=255)
                                                                                                      gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  <  2:374194/48‑1 (MQ=255)
                                                                                                      gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  >  1:374194/1‑48 (MQ=255)
                                                                                                        ccAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGc                     <  1:849523/139‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  <  1:449211/37‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                  >  2:449211/1‑37 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 >  1:1213468/1‑38 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                 <  2:1213468/38‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         >  1:484285/1‑62 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                         <  2:484285/62‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                 >  1:1210210/1‑70 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                 <  2:1210210/70‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTc                                                         >  1:588864/1‑94 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTc                                                         <  2:588864/94‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATc                >  2:838816/1‑135 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATc                <  1:838816/135‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga            >  2:1170300/1‑139 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGAAGATGGGCGTATTATAACCAACTaaaa                                                                                              >  2:528858/1‑57 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGAAGATGGGCGTATTATAACCAACTaaaa                                                                                              <  1:528858/57‑1 (MQ=255)
                                                                                                                    ccTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                              >  2:328790/1‑70 (MQ=255)
                                                                                                                    ccTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                              <  1:328790/70‑1 (MQ=255)
                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   <  2:1141324/63‑1 (MQ=255)
                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                   >  1:1141324/1‑63 (MQ=255)
                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGaagaag       >  1:390076/1‑139 (MQ=255)
                                                                                                                      tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGaagaag       <  1:518257/139‑1 (MQ=255)
                                                                                                                           cGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGaca  <  2:459381/139‑1 (MQ=255)
                                                                                                                               ccTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                              <  1:335117/75‑1 (MQ=255)
                                                                                                                               ccTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                              >  2:335117/1‑75 (MQ=255)
                                                                                                                                      |                                                                                                                               
CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACA  >  NC_000913/3545431‑3545692

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