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% 262.8 / NA 80*228* (TAG→TAAgntXDNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/44);  total (36/44)

CGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAAT  >  NC_000913/3545433‑3545697
                                                                                                                                    |                                                                                                                                    
cGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTc                                                                                                                                >  2:1137330/1‑139 (MQ=255)
   gCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgat                                                                                                                             >  2:1395681/1‑139 (MQ=255)
         tCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       <  2:1272285/139‑1 (MQ=255)
          cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      <  2:1467329/139‑1 (MQ=255)
          cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      <  1:1395681/139‑1 (MQ=255)
                   gTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                             <  1:1330877/139‑1 (MQ=255)
                     gATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                           >  1:152945/1‑139 (MQ=255)
                      aTTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCa                                                                                                          >  2:808173/1‑139 (MQ=255)
                       ttGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCaa                                                                                                         >  1:1495159/1‑139 (MQ=255)
                            gatgatGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaa                                                                                                    >  2:770978/1‑139 (MQ=255)
                             atgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                                   <  1:1360937/139‑1 (MQ=255)
                              tgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      <  1:1168462/119‑1 (MQ=255)
                              tgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      >  2:1168462/1‑119 (MQ=255)
                                 tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGt                                                                                               <  1:1137330/139‑1 (MQ=255)
                                        accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTTGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                        <  1:169059/139‑1 (MQ=255)
                                        accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                        >  2:1228635/1‑139 (MQ=255)
                                        accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                        <  2:1495159/139‑1 (MQ=255)
                                          caccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      <  1:1228635/139‑1 (MQ=255)
                                           accGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAg                                                                                     <  2:468930/139‑1 (MQ=255)
                                            ccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAgg                                                                                    <  1:765719/139‑1 (MQ=255)
                                             cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      <  1:1230159/104‑1 (MQ=255)
                                             cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      >  2:1230159/1‑104 (MQ=255)
                                             cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   <  1:770978/139‑1 (MQ=255)
                                             cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGCTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   <  1:1129850/139‑1 (MQ=255)
                                               gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                 <  2:331851/139‑1 (MQ=255)
                                                   tACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       >  1:1387334/1‑97 (MQ=255)
                                                   tACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       <  2:1387334/97‑1 (MQ=255)
                                                   tACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTAc                                                                             <  2:152945/139‑1 (MQ=255)
                                                    aCCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                 >  2:918554/1‑86 (MQ=255)
                                                    aCCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                 <  1:918554/86‑1 (MQ=255)
                                                       gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                 >  1:1461797/1‑83 (MQ=255)
                                                       gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCt                                                                                                                                 <  2:1461797/83‑1 (MQ=255)
                                                          gCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAg                                                                                     >  1:475771/1‑124 (MQ=255)
                                                          gCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAg                                                                                     <  2:475771/124‑1 (MQ=255)
                                                          gCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                 <  2:977/128‑1 (MQ=255)
                                                          gCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                 >  1:977/1‑128 (MQ=255)
                                                                  tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                             >  2:1088879/1‑92 (MQ=255)
                                                                  tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                             <  1:1088879/92‑1 (MQ=255)
                                                                  tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   <  2:1544748/118‑1 (MQ=255)
                                                                  tGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   >  1:1544748/1‑118 (MQ=255)
                                                                   gcgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGta                                                                                                                   >  2:292898/1‑85 (MQ=255)
                                                                   gcgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGta                                                                                                                   <  1:292898/85‑1 (MQ=255)
                                                                    cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                              <  1:1330711/89‑1 (MQ=255)
                                                                    cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                              >  2:1330711/1‑89 (MQ=255)
                                                                     gcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      >  2:1572293/1‑80 (MQ=255)
                                                                     gcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                      <  1:1572293/80‑1 (MQ=255)
                                                                         ttGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgct                                                       >  1:765325/1‑139 (MQ=255)
                                                                          tGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                              >  2:717403/1‑99 (MQ=255)
                                                                          tGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                              <  1:717403/99‑1 (MQ=255)
                                                                                    cAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCa                                            <  2:1276205/139‑1 (MQ=255)
                                                                                      aTGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      <  2:723345/95‑1 (MQ=255)
                                                                                      aTGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      >  1:723345/1‑95 (MQ=255)
                                                                                               cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       <  2:1129061/53‑1 (MQ=255)
                                                                                               cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       >  1:1129061/1‑53 (MQ=255)
                                                                                               cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   <  2:1094680/89‑1 (MQ=255)
                                                                                               cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   >  1:1094680/1‑89 (MQ=255)
                                                                                                 gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      <  2:1524814/84‑1 (MQ=255)
                                                                                                 gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      >  1:1524814/1‑84 (MQ=255)
                                                                                                      ccAGGTCTGTTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGc                          <  1:808173/139‑1 (MQ=255)
                                                                                                      ccAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGc                          <  2:765325/139‑1 (MQ=255)
                                                                                                            cTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       <  1:823398/40‑1 (MQ=255)
                                                                                                            cTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                       >  2:823398/1‑40 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                             <  2:1102650/47‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                             >  1:1102650/1‑47 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                           <  2:319202/49‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAAcc                                                                                                           >  1:319202/1‑49 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                              >  1:251112/1‑62 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                              <  2:251112/62‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                        >  2:347343/1‑68 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                        <  1:347343/68‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      >  1:515868/1‑70 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTa                                                                                      <  2:515868/70‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                 <  1:1569330/75‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                 >  2:1569330/1‑75 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgctgc                                                     <  1:1608906/103‑1 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgctgc                                                     >  2:1608906/1‑103 (MQ=255)
                                                                                                               gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga                 >  2:845991/1‑139 (MQ=255)
                                                                                                                    tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   >  2:52324/1‑68 (MQ=255)
                                                                                                                    tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                   <  1:52324/68‑1 (MQ=255)
                                                                                                                              cTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAAt  >  2:185243/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAAT  >  NC_000913/3545433‑3545697

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