Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,618,075 C→T *234* (TAG→TAA hda ← ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,0750CT97.3% 241.4 / ‑1.2 75coding (702/702 nt)hdaATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/35);  minor base . (1/1);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCCC  >  NC_000913/2617947‑2618209
                                                                                                                                |                                                                                                                                      
cATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCa                                                                                                                              >  2:1056140/1‑139 (MQ=255)
 aTCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTAc                                                                                                                                     >  2:370561/1‑131 (MQ=255)
 aTCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTAc                                                                                                                                     <  1:370561/131‑1 (MQ=255)
  tCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGa                                                                                                                            >  2:1111166/1‑139 (MQ=255)
          tGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttcttt                                                                                                                    >  1:1641164/1‑139 (MQ=255)
                  aTTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACg                                                                                                            >  2:392830/1‑139 (MQ=255)
                    tGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                          >  2:951954/1‑139 (MQ=255)
                    tGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                          >  1:1575006/1‑139 (MQ=255)
                    tGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                          <  1:1233970/139‑1 (MQ=255)
                      ccGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAAt                                                                                                        >  1:1203705/1‑139 (MQ=255)
                        ggATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAGACGGAATgg                                                                                                      <  1:1111166/139‑1 (MQ=255)
                        ggATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                      <  1:392830/139‑1 (MQ=255)
                          aTGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTc                                                                                                    <  2:1203705/139‑1 (MQ=255)
                              gACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCt                                                                                                >  2:263999/1‑139 (MQ=255)
                                       cACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTgcg                                                                                       <  2:839795/139‑1 (MQ=255)
                                               aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGa                                                                                                                            >  2:1298809/1‑94 (MQ=255)
                                               aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGa                                                                                                                            <  1:1298809/94‑1 (MQ=255)
                                               aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATc                                                                               >  1:1410052/1‑139 (MQ=255)
                                                   ggCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                           <  2:1410052/139‑1 (MQ=255)
                                                   ggCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                           <  2:1575006/139‑1 (MQ=255)
                                                    gCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                          >  2:451376/1‑139 (MQ=255)
                                                    gCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                          >  1:53565/1‑139 (MQ=255)
                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  1:951954/139‑1 (MQ=255)
                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  2:1341491/139‑1 (MQ=255)
                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  1:169577/139‑1 (MQ=255)
                                                     ccTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         >  1:1341491/1‑139 (MQ=255)
                                                          gAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATc                                                                    >  2:841844/1‑139 (MQ=255)
                                                               cTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACt                                                               >  1:1586134/1‑139 (MQ=255)
                                                                  tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  1:646313/126‑1 (MQ=255)
                                                                  tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         >  2:646313/1‑126 (MQ=255)
                                                                  tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                                                            >  1:1261234/1‑139 (MQ=255)
                                                                   ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATc                                                           >  1:1656892/1‑139 (MQ=255)
                                                                   ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATc                                                           >  1:149074/1‑139 (MQ=255)
                                                                    gATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATcc                                                          >  2:63389/1‑139 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAATGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  1:385661/139‑1 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCtt                                                                                               >  2:1565967/1‑101 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCtt                                                                                               <  1:1565967/101‑1 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  2:448298/139‑1 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  1:1056140/139‑1 (MQ=255)
                                                                     aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  1:841844/139‑1 (MQ=255)
                                                                        aaCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACg                                                      >  2:1457781/1‑139 (MQ=255)
                                                                            tCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCAt                                                  >  1:644764/1‑139 (MQ=255)
                                                                              gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  1:938096/114‑1 (MQ=255)
                                                                              gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         >  2:938096/1‑114 (MQ=255)
                                                                              gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAAc                                                                >  1:1321713/1‑123 (MQ=255)
                                                                              gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAAc                                                                <  2:1321713/123‑1 (MQ=255)
                                                                               cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACACTAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                                                  <  1:879844/120‑1 (MQ=255)
                                                                               cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                                                  >  2:879844/1‑120 (MQ=255)
                                                                                        gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTgcg                                      <  2:644764/139‑1 (MQ=255)
                                                                                             ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAttt                                 <  2:1656892/139‑1 (MQ=255)
                                                                                             ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAttt                                 <  1:263999/139‑1 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                               <  1:801778/60‑1 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACaa                                                                                                               >  2:801778/1‑60 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                          <  1:1243037/65‑1 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGa                                                                                                          >  2:1243037/1‑65 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  2:160087/114‑1 (MQ=255)
                                                                                              gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         >  1:160087/1‑114 (MQ=255)
                                                                                                  tCGCGCCGCATCCGACAATAAACACCTTATTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         >  1:778994/1‑109 (MQ=255)
                                                                                                  tCGCGCCGCATCCGACAATAAACACCTTATTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                         <  2:778994/109‑1 (MQ=255)
                                                                                                       ccgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAg                       <  2:1261234/139‑1 (MQ=255)
                                                                                                        cgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGc                                           >  1:736940/1‑118 (MQ=255)
                                                                                                        cgcATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGc                                           <  2:736940/118‑1 (MQ=255)
                                                                                                            tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                  <  2:296469/139‑1 (MQ=255)
                                                                                                                    tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                          <  1:391023/75‑1 (MQ=255)
                                                                                                                    tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                          >  2:391023/1‑75 (MQ=255)
                                                                                                                           cTTATTTGCAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                  >  1:1191716/1‑124 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                           <  1:1440439/67‑1 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                           >  2:1440439/1‑67 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa        >  2:625815/1‑134 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGaa        <  1:625815/134‑1 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGcc   >  2:1057666/1‑139 (MQ=255)
                                                                                                                           cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACCATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt                  <  2:1191716/124‑1 (MQ=255)
                                                                                                                            ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                         <  1:1021842/68‑1 (MQ=255)
                                                                                                                            ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGccc  >  2:1206168/1‑139 (MQ=255)
                                                                                                                            ttATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTACGCGGTAATCGACGCa                                                                         >  2:1021842/1‑68 (MQ=255)
                                                                                                                                |                                                                                                                                      
CATCGCTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCCC  >  NC_000913/2617947‑2618209

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