Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,992,898 C→T S582S (AGC→AGT cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,8980CT98.5% 202.4 / ‑4.3 66S582S (AGC→AGTcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/35);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

CACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAA  >  NC_000913/3992762‑3993031
                                                                                                                                        |                                                                                                                                     
caccgacaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGttt                                                                                                                                     >  1:915612/1‑139 (MQ=255)
caccgacaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGttt                                                                                                                                     >  1:1570768/1‑139 (MQ=255)
   cgacaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTgg                                                                                                                                  >  1:103691/1‑139 (MQ=255)
     acaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCg                                                                                                                                >  2:658553/1‑139 (MQ=255)
     acaccgaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCg                                                                                                                                <  1:1704308/139‑1 (MQ=255)
          gaAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAAt                                                                                                                           >  1:39161/1‑139 (MQ=255)
           aaGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATc                                                                                                                          >  2:1737564/1‑139 (MQ=255)
                 ctctACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAAtt                                                                                                                    >  2:898080/1‑139 (MQ=255)
                   ctACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGc                                                                                                                  <  1:161422/139‑1 (MQ=255)
                     aCAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCt                                                                                                                <  2:103691/139‑1 (MQ=255)
                           cGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTggtagg                                                                                                          <  1:658553/139‑1 (MQ=255)
                           cGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTggtagg                                                                                                          <  1:1737564/139‑1 (MQ=255)
                              gtgAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGgtaggtag                                                                                                       <  2:371278/139‑1 (MQ=255)
                              gtgAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGgtaggtag                                                                                                       >  1:189296/1‑139 (MQ=255)
                                   aTCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTg                                                                                                  <  2:55735/139‑1 (MQ=255)
                                   aTCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTg                                                                                                  <  1:524914/139‑1 (MQ=255)
                                      cGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTggtagg                                                                                                          <  2:1682195/128‑1 (MQ=255)
                                      cGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTggtagg                                                                                                          >  1:1682195/1‑128 (MQ=255)
                                      cGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGAcc                                                                                               <  1:1478509/139‑1 (MQ=255)
                                      cGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGAcc                                                                                               <  1:117480/139‑1 (MQ=255)
                                        ccATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACctg                                                                                             >  2:1467641/1‑139 (MQ=255)
                                                       cGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCg                                                                              >  2:1804361/1‑139 (MQ=255)
                                                             ccTGGAATATGACCCGGCAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGGTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAAc                                                                        <  1:1407140/139‑1 (MQ=255)
                                                                ggAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTc                                                                               <  1:659780/129‑1 (MQ=255)
                                                                ggAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTc                                                                               >  2:659780/1‑129 (MQ=255)
                                                                ggAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      >  2:1192514/1‑138 (MQ=255)
                                                                ggAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      <  1:1192514/138‑1 (MQ=255)
                                                                 gAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      >  2:1107740/1‑137 (MQ=255)
                                                                 gAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      <  1:1107740/137‑1 (MQ=255)
                                                                 gAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAg                                                                    >  2:1665398/1‑139 (MQ=255)
                                                                    tatGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGc                                                                 <  2:1429309/139‑1 (MQ=255)
                                                                      tGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGt                                                               >  2:415130/1‑139 (MQ=255)
                                                                       gACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTa                                                              <  1:1467641/139‑1 (MQ=255)
                                                                         cccGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACg                                                            <  1:1804361/139‑1 (MQ=255)
                                                                         cccGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACg                                                            <  1:1665398/139‑1 (MQ=255)
                                                                            gACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      >  1:786534/1‑126 (MQ=255)
                                                                            gACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGa                                                                      <  2:786534/126‑1 (MQ=255)
                                                                            gACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTg                                                         >  2:542569/1‑139 (MQ=255)
                                                                             aCAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGc                                                        >  2:1066067/1‑139 (MQ=255)
                                                                                 cggcgTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACtt                                                    <  2:411275/139‑1 (MQ=255)
                                                                                         cGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTg                                                                                                               <  1:677711/72‑1 (MQ=255)
                                                                                         cGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTg                                                                                                               >  2:677711/1‑72 (MQ=255)
                                                                                           cAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAg                                          <  1:1066067/139‑1 (MQ=255)
                                                                                           cAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAg                                          <  2:915612/139‑1 (MQ=255)
                                                                                                  gtggtgCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCga                                   >  2:1055269/1‑139 (MQ=255)
                                                                                                   tggtgCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCgat                                  >  2:1366378/1‑139 (MQ=255)
                                                                                                     gtgCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCg                                            <  1:268188/127‑1 (MQ=255)
                                                                                                     gtgCATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCg                                            >  2:268188/1‑127 (MQ=255)
                                                                                                        cATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCga                                   <  2:1430244/133‑1 (MQ=255)
                                                                                                        cATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCga                                   >  1:1430244/1‑133 (MQ=255)
                                                                                                        cATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCg                             <  1:1055269/139‑1 (MQ=255)
                                                                                                        cATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCg                             >  2:792913/1‑139 (MQ=255)
                                                                                                        cATTTCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCg                             <  2:39161/139‑1 (MQ=255)
                                                                                                          tttCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGc                                                                 >  2:151462/1‑101 (MQ=255)
                                                                                                          tttCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGc                                                                 <  1:151462/101‑1 (MQ=255)
                                                                                                          tttCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGaa                           <  1:1046680/139‑1 (MQ=255)
                                                                                                          tttCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGaa                           <  2:1276426/139‑1 (MQ=255)
                                                                                                            tCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGaa                                                                                                                            <  2:500207/40‑1 (MQ=255)
                                                                                                            tCGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGaa                                                                                                                            >  1:500207/1‑40 (MQ=255)
                                                                                                             cGATTTCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGcc                        <  1:898080/139‑1 (MQ=255)
                                                                                                                tttCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGgtaggtag                                                                                                       >  2:128735/1‑57 (MQ=255)
                                                                                                                tttCCGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGgtaggtag                                                                                                       <  1:128735/57‑1 (MQ=255)
                                                                                                                    cGTAAGCTGGATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGaaaa                 <  2:189296/139‑1 (MQ=255)
                                                                                                                             gATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAg                                          <  2:152193/105‑1 (MQ=255)
                                                                                                                             gATGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAg                                          >  1:152193/1‑105 (MQ=255)
                                                                                                                              aTGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCgg       <  1:1366378/139‑1 (MQ=255)
                                                                                                                              aTGTCTTCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCgg       >  2:959307/1‑139 (MQ=255)
                                                                                                                                   ttCAGTTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCaaaa  >  1:333379/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGGTAGCGTTGACCTGCTGTACCGCAACTCGTGGAACGAAGTGCGTACGCTGCACTTCAACGGCGAGCAATCGATGATCGAAGCCCTGAAAACTATTCTCGGCAAAA  >  NC_000913/3992762‑3993031

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