Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,406,145 T→C I56T (ATC→ACC)  smrA → DNA endonuclease SmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,406,1450TC98.6% 251.9 / ‑6.2 72I56T (ATC→ACC) smrADNA endonuclease SmrA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/38);  minor base G (1/0);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATG  >  NC_000913/1406008‑1406281
                                                                                                                                         |                                                                                                                                        
tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCAccc                                                                                                                                        >  2:1544205/1‑140 (MQ=255)
  aCGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCa                                                                                                                                      <  1:1109092/140‑1 (MQ=255)
    gCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACt                                                                                                                                    <  2:312051/140‑1 (MQ=255)
                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTgg                                                                                                                      >  2:213284/1‑140 (MQ=255)
                    aGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAg                                                                                                                    >  1:521590/1‑140 (MQ=255)
                          tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTcgg                                                                                                              >  1:1363537/1‑140 (MQ=255)
                           gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTcggc                                                                                                             >  1:575447/1‑140 (MQ=255)
                           gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTcggc                                                                                                             >  2:1258957/1‑140 (MQ=255)
                               cGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGa                                                                                                         <  1:213284/140‑1 (MQ=255)
                                      ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTg                                                                                                  <  1:1258957/140‑1 (MQ=255)
                                      ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTg                                                                                                  >  2:1637816/1‑140 (MQ=255)
                                        gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCa                                                                                                <  2:968968/140‑1 (MQ=255)
                                           gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACa                                                                                             >  1:84583/1‑140 (MQ=255)
                                             ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg                                                                                           <  1:1544205/140‑1 (MQ=255)
                                              cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg                                                                                          <  2:1220193/140‑1 (MQ=255)
                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                                                        >  1:1170266/1‑140 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg                                                                                         <  1:916605/138‑1 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg                                                                                         >  2:916605/1‑138 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                                                       <  2:353281/140‑1 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                                                       >  2:1191140/1‑140 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                                                       <  2:661009/140‑1 (MQ=255)
                                                 tGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                                                       <  2:708921/140‑1 (MQ=255)
                                                   gCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGc                                                                                     <  1:953142/140‑1 (MQ=255)
                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTGAGTCAGCCGCTGGAGTTTAGGCGGGAAGGGTTGCAACATGGGGTGCTgg                                                                                  <  1:1183403/140‑1 (MQ=255)
                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTgg                                                                                  >  1:1149418/1‑140 (MQ=255)
                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                <  1:1438082/140‑1 (MQ=255)
                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                >  1:1570584/1‑140 (MQ=255)
                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                <  1:1620058/140‑1 (MQ=255)
                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                >  2:1201669/1‑140 (MQ=255)
                                                           cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg                                                                                           <  2:406855/126‑1 (MQ=255)
                                                           cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg                                                                                           >  1:406855/1‑126 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGGAGGGTTTCAACATGGGGGGCTGGATAAGCTgcg                                                                        >  2:1659548/1‑140 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  2:1363537/140‑1 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  2:274023/140‑1 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  2:1412684/140‑1 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  1:1242174/140‑1 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  2:1170266/140‑1 (MQ=255)
                                                                aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                        <  2:84583/140‑1 (MQ=255)
                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAg                                                                             >  1:185213/1‑134 (MQ=255)
                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAg                                                                             <  2:185213/134‑1 (MQ=255)
                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                                                       <  1:1630125/140‑1 (MQ=255)
                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                                                       <  1:29378/140‑1 (MQ=255)
                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                                                       <  2:521590/140‑1 (MQ=255)
                                                                  ccAACGTGCGCCGCAGCGTATTGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCa                                                                      <  2:621762/140‑1 (MQ=255)
                                                                      cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGGTGCAACATGGGGTGCTGGATAAGCTGCGCAGTgg                                                                  >  2:159611/1‑140 (MQ=255)
                                                                            ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                                                                 <  1:1340098/135‑1 (MQ=255)
                                                                            ccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                                                                 >  2:1340098/1‑135 (MQ=255)
                                                                                  cGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCa                                                      <  1:1191140/140‑1 (MQ=255)
                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAAc                                                    >  1:1523942/1‑140 (MQ=255)
                                                                                        gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGc                                                <  1:933555/140‑1 (MQ=255)
                                                                                            cGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGc                                                <  2:662677/136‑1 (MQ=255)
                                                                                            cGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGc                                                >  1:662677/1‑136 (MQ=255)
                                                                                            cGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGc                                            >  2:420305/1‑140 (MQ=255)
                                                                                              cTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCAGCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt                                          >  2:1478137/1‑140 (MQ=255)
                                                                                               tGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTg                                         >  2:1351610/1‑140 (MQ=255)
                                                                                                gCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGGTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGa                                        >  2:1592224/1‑140 (MQ=255)
                                                                                                 cAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGaa                                       <  1:770724/140‑1 (MQ=255)
                                                                                                     ttGATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCtt                                   <  2:1437842/140‑1 (MQ=255)
                                                                                                       gATAATTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCtttt                                 >  1:1150229/1‑140 (MQ=255)
                                                                                                           aTTTCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTgcgc                             >  2:1424813/1‑140 (MQ=255)
                                                                                                            tttCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGcc                            >  2:507910/1‑140 (MQ=255)
                                                                                                              tCCTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAg                          >  2:438050/1‑140 (MQ=255)
                                                                                                               ccTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGc                         <  1:1201669/140‑1 (MQ=255)
                                                                                                                cTCACCACCGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGcc                        >  2:513118/1‑140 (MQ=255)
                                                                                                                   accaccGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAATCTGCGCCGGGGTAAAAATCCCCAACAGGAGAGCCTGAATCTTTTGGGCCAGCCGGt                     >  2:1586323/1‑140 (MQ=255)
                                                                                                                   accaccGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCa                           <  2:348097/134‑1 (MQ=255)
                                                                                                                   accaccGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCa                           >  1:348097/1‑134 (MQ=255)
                                                                                                                   accaccGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCCTTTGCGCCAGCCGGt                     >  2:1187985/1‑140 (MQ=255)
                                                                                                                       ccGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGaa                 <  2:1570584/140‑1 (MQ=255)
                                                                                                                        cGGATTTCTCGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAag                >  1:1506992/1‑140 (MQ=255)
                                                                                                                               ctcGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCg         <  1:507910/140‑1 (MQ=255)
                                                                                                                                 cGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCa       >  2:599022/1‑140 (MQ=255)
                                                                                                                                 cGACATCACCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCa       <  1:1659548/140‑1 (MQ=255)
                                                                                                                                    catcacCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCaaaa    <  1:420305/140‑1 (MQ=255)
                                                                                                                                    catcacCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCaaaa    >  1:1360131/1‑140 (MQ=255)
                                                                                                                                     atcaccCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCaaaa    >  2:198758/1‑139 (MQ=255)
                                                                                                                                     atcaccCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCaaaa    <  1:198758/139‑1 (MQ=255)
                                                                                                                                      tcaccCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATg  <  2:575447/140‑1 (MQ=255)
                                                                                                                                      tcaccCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATg  <  2:1360131/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATG  >  NC_000913/1406008‑1406281

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