Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,368,966 C→T 5.5% R56C (CGT→TGT)  pspB → psp operon transcription co‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,368,9660CT 5.5% 225.6 / 2.5 72R56C (CGT→TGT) pspBpsp operon transcription co‑activator
Reads supporting (aligned to +/- strand):  ref base C (30/38);  new base T (2/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CCGTTAACCATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTATTTT  >  NC_000913/1368825‑1369114
                                                                                                                                             |                                                                                                                                                    
ccGTTAACCATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTc                                                                                                                                                <  2:215246/148‑1 (MQ=255)
 cGTTAACCATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAATGTATTCAg                                                                                                                                              >  4:169931/1‑149 (MQ=255)
 cGTTAACCATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAATGTATTCAg                                                                                                                                              >  4:169934/1‑149 (MQ=255)
       ccATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCTTTCTGGTCGCAGTGAATTGTCGCACAGTGATCAGCAGCGATTAGCGCAACTTGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTg                                                                                                                                        <  1:153056/149‑1 (MQ=255)
       ccATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTg                                                                                                                                        <  7:381966/149‑1 (MQ=255)
            tttGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAAt                                                                                                                                    <  4:163071/148‑1 (MQ=255)
                tgCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCt                                                                                                                                         <  6:174382/139‑1 (MQ=255)
                tgCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCt                                                                                                                                         >  5:174382/1‑139 (MQ=255)
                tgCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTAtt                                                                                                                               <  3:62285/149‑1 (MQ=255)
                 gCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTc                                                                                                                              >  6:36994/1‑149 (MQ=255)
                  cTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTc                                                                                                                              <  6:175584/148‑1 (MQ=255)
                   tgtttgttTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCtt                                                                                                                            >  8:223937/1‑149 (MQ=255)
                     tttgtttTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTc                                                                                                                                                >  2:261807/1‑127 (MQ=255)
                     tttgtttTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTc                                                                                                                                                <  1:261807/127‑1 (MQ=255)
                     tttgtttTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCGGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGTAATCTATTCTTGa                                                                                                                          >  3:23937/1‑149 (MQ=255)
                           ttACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAg                                                                                                                     <  6:200219/148‑1 (MQ=255)
                            tACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGa                                                                                                                    <  6:325396/148‑1 (MQ=255)
                                      gTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAATGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATcc                                                                                                              <  3:169931/144‑1 (MQ=255)
                                      gTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAATGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATcc                                                                                                              <  3:169934/144‑1 (MQ=255)
                                      gTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTa                                                                                                                                 >  2:316124/1‑125 (MQ=255)
                                      gTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTa                                                                                                                                 <  1:316124/125‑1 (MQ=255)
                                          tGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGa                                                                                                      <  7:223937/148‑1 (MQ=255)
                                            gCTGCATTACAGCAATCGTTCTGGTCGCAGTGACTTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGa                                                                                                      <  5:231764/146‑1 (MQ=255)
                                                gaTTACAGCATTCGTTCTGTTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGAGCAAGTGGCTGATGAAGCAAAAGTGATGGGCGAACGTATTCAGGCTCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCg                                                                                               <  8:448732/148‑1 (MQ=255)
                                                 aTTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                                                        >  3:399064/1‑139 (MQ=255)
                                                 aTTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                                                        <  4:399064/139‑1 (MQ=255)
                                                     cAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGCACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAAt                                                                                          >  5:445693/1‑149 (MQ=255)
                                                         aaTCGTTCTGGTCGCAGTGAATTCTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAACCAAAACGGATCCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAACGCTAATggag                                                                                      >  7:135863/1‑147 (MQ=255)
                                                          aTCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCgg                                                                                     <  5:36994/149‑1 (MQ=255)
                                                             gTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                                                        <  6:128472/127‑1 (MQ=255)
                                                             gTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                                                        >  5:128472/1‑127 (MQ=255)
                                                             gTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCg                                                                                      <  2:110897/145‑1 (MQ=255)
                                                             gTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCg                                                                                      >  1:110897/1‑145 (MQ=255)
                                                              ttCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCAt                                                                                 >  5:8136/1‑149 (MQ=255)
                                                                      gctgTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACTGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCATAACTTCCTTACTGGAGGGATCGCTAATGGCGGGCATTAATCTCa                                                                         <  7:408111/146‑1 (MQ=255)
                                                                      gCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGACCTGGAGGGATCGCTAATGGCGGGCATTAAtctc                                                                          >  3:16377/1‑148 (MQ=255)
                                                                        aGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                                                              <  6:360224/126‑1 (MQ=255)
                                                                        aGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                                                              >  5:360224/1‑126 (MQ=255)
                                                                         gTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAAt                                                                       >  8:165245/1‑148 (MQ=255)
                                                                         gTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAAt                                                                       >  8:390058/1‑148 (MQ=255)
                                                                            aaTTGTCGCAAAGTGCGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaa                                                                    <  7:387237/148‑1 (MQ=255)
                                                                            aaTTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaa                                                                    <  8:225980/148‑1 (MQ=255)
                                                                            aaTTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaa                                                                    <  3:87284/148‑1 (MQ=255)
                                                                            aaTTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaa                                                                    >  1:449381/1‑148 (MQ=255)
                                                                             aTTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaaa                                                                   <  5:178878/148‑1 (MQ=255)
                                                                                gTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAAt                                                                >  1:428352/1‑148 (MQ=255)
                                                                                 tCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTa                                                              >  8:213766/1‑149 (MQ=255)
                                                                                    cAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGATGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATgg                                                           >  6:262511/1‑149 (MQ=255)
                                                                                    cAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATgg                                                           <  4:175151/149‑1 (MQ=255)
                                                                                     aaaGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGc                                                          <  2:449381/149‑1 (MQ=255)
                                                                                       aGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGc                                                          >  1:80913/1‑147 (MQ=255)
                                                                                       aGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGc                                                          <  2:80913/147‑1 (MQ=255)
                                                                                         tGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTAt                                                             <  5:426862/142‑1 (MQ=255)
                                                                                         tGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTAt                                                             >  6:426862/1‑142 (MQ=255)
                                                                                                     aTTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGc                                                          <  2:394559/133‑1 (MQ=255)
                                                                                                     aTTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGc                                                          >  1:394559/1‑133 (MQ=255)
                                                                                                      ttAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACAACCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGc                                         <  5:262511/149‑1 (MQ=255)
                                                                                                        aGCGCAACTGGCTGTTGAAGCATATCGGATGCGCGAACGTATTCAGGCGCTGTAATATATTGTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCa                                        <  8:135863/148‑1 (MQ=255)
                                                                                                         gcgcAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATg                                      >  8:44261/1‑149 (MQ=255)
                                                                                                          cgcAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATgg                                     <  2:428352/149‑1 (MQ=255)
                                                                                                                 ggCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAATATTATGGCGTATTCCACAGCAgg                                           <  8:27549/136‑1 (MQ=255)
                                                                                                                 ggCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAgg                                           >  7:27549/1‑136 (MQ=255)
                                                                                                                        gaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTc                                                                                                                              >  6:131988/1‑46 (MQ=255)
                                                                                                                        gaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTc                                                                                                                              <  5:131988/46‑1 (MQ=255)
                                                                                                                              aaaCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCgg                  >  8:283051/1‑148 (MQ=255)
                                                                                                                               aaCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGa                >  1:188526/1‑149 (MQ=255)
                                                                                                                                   gATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGc            <  7:283051/149‑1 (MQ=255)
                                                                                                                                       cgcgAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCaa         <  4:23937/148‑1 (MQ=255)
                                                                                                                                        gcgAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAAc        >  8:269358/1‑148 (MQ=255)
                                                                                                                                         cgAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTa      <  4:16377/149‑1 (MQ=255)
                                                                                                                                          gAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGCACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTAt     >  1:30/1‑149 (MQ=255)
                                                                                                                                            aCGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTAttt   <  6:8136/149‑1 (MQ=255)
                                                                                                                                             cGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTAtttt  <  7:390058/149‑1 (MQ=255)
                                                                                                                                             cGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTAtttt  <  1:227795/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCGTTAACCATTTTTGTGCTGTTTGTTTTACCGATCTGGTTATGGCTGCATTACAGCAATCGTTCTGGTCGCAGTGAATTGTCGCAAAGTGAGCAGCAGCGATTAGCGCAACTGGCTGATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCATGGTCCGCGGCGTCTGCGCCGGGATTGCCAACTATTTT  >  NC_000913/1368825‑1369114

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: