Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,029 A→G 7.6% V37V (GTA→GTG) ‡ rsxC → SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,707,0290AG 7.6% 221.9 / 6.6 78V37V (GTA→GTG) ‡rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base G (3/3);  total (41/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 = 9.26e-01

TGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGG  >  NC_000913/1706880‑1707177
                                                                                                                                                     |                                                                                                                                                    
tGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTa                                                                                                                                                      >  2:1269864/1‑150 (MQ=255)
 gAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAc                                                                                                                                                     >  2:956112/1‑150 (MQ=255)
 gAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAc                                                                                                                                                     >  2:1161713/1‑150 (MQ=255)
  aaCACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTAcc                                                                                                                                                    <  1:1569494/150‑1 (MQ=255)
      ccATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTg                                                                                                                                                >  2:935691/1‑150 (MQ=255)
        aTTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgc                                                                                                                                              <  2:1174771/150‑1 (MQ=255)
        aTTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgc                                                                                                                                              <  2:575924/150‑1 (MQ=255)
         ttCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGgcg                                                                                                                                             >  1:964152/1‑150 (MQ=255)
           cccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCa                                                                                                                                           >  1:189183/1‑150 (MQ=255)
           cccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCa                                                                                                                                           >  1:147537/1‑150 (MQ=255)
            ccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCACGCGCAg                                                                                                                                          >  1:230891/1‑150 (MQ=255)
            ccGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCAg                                                                                                                                          <  1:494548/150‑1 (MQ=255)
             cGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGc                                                                                                                                         <  2:147537/150‑1 (MQ=255)
               tGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCAgggg                                                                                                                                       <  1:494860/150‑4 (MQ=255)
                 cGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGAATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGGGTACCTGGCGCAggggtg                                                                                                                                     <  2:766557/150‑6 (MQ=255)
                    aTCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGt                                                                                                                                  <  2:890587/150‑1 (MQ=255)
                     tCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGtt                                                                                                                                 <  1:736121/150‑1 (MQ=255)
                        ttCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTAtt                                                                                                                              <  1:337299/150‑1 (MQ=255)
                         tCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTc                                                                                                                             <  1:1201216/150‑1 (MQ=255)
                          cccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTcc                                                                                                                            <  1:956112/150‑1 (MQ=255)
                          cccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTcc                                                                                                                            <  1:729245/150‑1 (MQ=255)
                           ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCa                                                                                                                           >  1:1205910/1‑150 (MQ=255)
                           ccGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCa                                                                                                                           >  1:153007/1‑150 (MQ=255)
                                  acacCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAAc                                                                                                                    >  2:1459938/1‑150 (MQ=255)
                                     ccATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGc                                                                                                                 >  1:1453018/1‑150 (MQ=255)
                                      cATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCa                                                                                                                <  1:1319283/150‑1 (MQ=255)
                                        tGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCata                                                                                                              <  2:859401/150‑1 (MQ=255)
                                            tAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTgg                                                                                                          <  1:1161713/150‑1 (MQ=255)
                                              aGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgcg                                                                                                        <  2:964152/150‑1 (MQ=255)
                                               gTTATTCTCTGCATTCAGAAAAAATAAATTCTTGGATTTCAACGCCGGCATCCTTCCACCGGAGATGAAAACCCAGTGCAACGGTACACCCCTGCGCCAGCTACGCCTGGCGAAGCGTTTTGTTAATCCACTGAATCAGCATATTGgcgc                                                                                                       <  2:606265/150‑1 (MQ=255)
                                               gTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGgcgc                                                                                                       >  2:538217/1‑150 (MQ=255)
                                                     ctctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAgg                                                                                                 >  2:588774/1‑150 (MQ=255)
                                                       ctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                                                                               <  2:941026/150‑1 (MQ=255)
                                                       ctGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTg                                                                                               >  1:1495191/1‑150 (MQ=255)
                                                         gCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAg                                                                                             >  2:814860/1‑150 (MQ=255)
                                                            ttCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtg                                                                                          >  2:1429692/1‑150 (MQ=255)
                                                             tCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgt                                                                                         <  1:935691/150‑1 (MQ=255)
                                                              cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                                                                        >  2:565495/1‑150 (MQ=255)
                                                              cAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTtgtg                                                                                        >  2:1460784/1‑150 (MQ=255)
                                                               aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTcc                                                                                                                            >  2:722607/1‑113 (MQ=255)
                                                               aGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTcc                                                                                                                            <  1:722607/113‑1 (MQ=255)
                                                                gAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCGGCGGTTTGTTATTCCACGGAAACAGCATATTGGCGCGGAAGGTGAGGtgtggg                                                                                      >  1:195294/1‑148 (MQ=255)
                                                                   aaaataaaatCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTa                                                                                   <  1:521840/150‑1 (MQ=255)
                                                                           atCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGc                                                                           <  1:616677/150‑1 (MQ=255)
                                                                              tGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGAt                                                                        >  1:1081499/1‑150 (MQ=255)
                                                                               gggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATa                                                                       >  2:1111822/1‑150 (MQ=255)
                                                                               gggATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATa                                                                       <  1:1429314/150‑1 (MQ=255)
                                                                                  aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg                                                                    <  2:1081499/150‑1 (MQ=255)
                                                                                  aTTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAg                                                                    <  2:1453018/150‑1 (MQ=255)
                                                                                   tttCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGt                                                                   <  2:1087640/150‑1 (MQ=255)
                                                                                     tCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAt                                                                 >  2:1408100/1‑150 (MQ=255)
                                                                                      cAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTAtt                                                                >  2:1366622/1‑150 (MQ=255)
                                                                                       aaCGGCGGCATCCGTCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTg                                                               >  2:1566382/1‑150 (MQ=255)
                                                                                         cggcggCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTgcg                                                             >  2:1192686/1‑150 (MQ=255)
                                                                                               catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCa                                                       <  1:1429692/150‑1 (MQ=255)
                                                                                               catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCa                                                       <  2:189183/150‑1 (MQ=255)
                                                                                               catccatcCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGAGATAAAGTATTGCGCGGCCa                                                       >  2:737535/1‑150 (MQ=255)
                                                                                                        aCCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTAc                                              <  1:1138176/150‑1 (MQ=255)
                                                                                                         ccGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTAcc                                             <  1:1566382/150‑1 (MQ=255)
                                                                                                              gaTGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTgg                                        <  1:1105165/150‑1 (MQ=255)
                                                                                                                tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTc                                      >  1:1251788/1‑150 (MQ=255)
                                                                                                                tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTc                                      <  1:459859/150‑1 (MQ=255)
                                                                                                                tGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTc                                      >  2:55806/1‑150 (MQ=255)
                                                                                                                 gAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTTCGCGGCCAGCCGCTTACCCGTGGTcg                                     >  1:635827/1‑150 (MQ=255)
                                                                                                                      cccAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCa                                <  1:538217/150‑1 (MQ=255)
                                                                                                                          gTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAt                            >  1:705918/1‑150 (MQ=255)
                                                                                                                           tCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGTGTTAGCGTCGGTGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtg                           >  2:693392/1‑150 (MQ=255)
                                                                                                                           gccAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCGGAAGGGGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtg                           <  1:28306/149‑1 (MQ=255)
                                                                                                                            ccAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtgc                          >  1:1343733/1‑150 (MQ=255)
                                                                                                                             cAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAAtgct                         >  1:1220160/1‑150 (MQ=255)
                                                                                                                                cGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGcc                      <  1:1459938/150‑1 (MQ=255)
                                                                                                                                       ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAc               >  2:869298/1‑150 (MQ=255)
                                                                                                                                       ccccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAc               >  2:1261338/1‑150 (MQ=255)
                                                                                                                                       acgcTGCGCCAGGGGTACATGGCGCAGCGTTTTGTTATTCCACTGAAAAAGCATATTGGAGCTGAACGAGAGTTGTGCCTTAGCGTCGGCGATAAAGTATTGCGACGCCAGCCGCTTACCCGTGGTCGCGGAAAAATGATGCATGTTAAc               >  2:296006/4‑150 (MQ=255)
                                                                                                                                        cgcTGCGCCAGGGGTACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcg              >  1:1304730/3‑150 (MQ=255)
                                                                                                                                        cccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcg              <  1:588774/150‑1 (MQ=255)
                                                                                                                                         gcTGCGCCAGGGGTACCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgc             >  1:13459/2‑150 (MQ=255)
                                                                                                                                         ccTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCAcgc             >  1:5014/1‑150 (MQ=255)
                                                                                                                                                    tACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCgg  >  1:105302/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAACACCATTCCCGTGCGTATCATTCCCGTGGAACACCATGCTTAAGTTATTCTCTGCATTCAGAAAAAATAAAATCTGGGATTTCAACGGCGGCATCCATCCACCGGAGATGAAAACCCAGTCCAACGGTACACCCCTGCGCCAGGTACCCCTGGCGCAGCGTTTTGTTATTCCACTGAAACAGCATATTGGCGCTGAAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGG  >  NC_000913/1706880‑1707177

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