Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 678,688 A→C 5.9% N425T (AAT→ACT)  djlB → putative HscC co‑chaperone, uncharacterized J domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913678,6880AC 5.9% 229.7 / 6.6 68N425T (AAT→ACT) djlBputative HscC co‑chaperone, uncharacterized J domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TTCACTGGTCGAGATGGTTTTGCACAAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGGAAAATACTCCATCTGTTCGGGTGACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTAATCTTTCCGCATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTA  >  NC_000913/678464‑678951
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                       
ttCACTGGTCGAGATGGTTTTGCACAAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATg                                                                                                                                                                                                                        >  2:507022/1‑274 (MQ=255)
    cTGGTCGAGATGGTTTTGCACAAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCg                                                                                                                                                                                                                    <  4:507021/274‑1 (MQ=255)
                  tttGCACAAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGAc                                                                                                                                                                                                               >  4:478677/1‑265 (MQ=255)
                     tcccaAAGCCCGTAAGCCGTGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCCAGTTGCCAAACTATTAGTGCAGGCCGGCCAACCTACATAATAATCATTAATGCCGATTGTCAGCCGCCTTTATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATGCAGGCGGCGCATTTAACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCa                                                                                                                                                                                                   <  4:58121/271‑1 (MQ=255)
                      cacaAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAtg                                                                                                                                                                                                           <  2:478678/265‑1 (MQ=255)
                       acaAAGGACGTAAGCCGCGGATCTAACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAgg                                                                                                                                                                                                 <  2:476694/274‑1 (MQ=255)
                                  aaGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAAGGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGcat                                                                                                                                                                                     <  4:248439/275‑1 (MQ=255)
                                          gATCTCACTCACGCGTGATCCAGAACATACCTTATGGCCATTATAACTAGTTGCCCAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATGGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGAAAATTGTTGAACCGGAGCCGCAACCGCATCCAATGCCCGGTGACGATGGGGGGCCAGGCTGGCTGGGCATCATTaaaa                                                                                                                                                                             >  4:285695/1‑275 (MQ=255)
                                              tcactcacGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTggg                                                                                                                                                                                        <  2:285696/260‑1 (MQ=255)
                                                            tACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTa                                                                                                                                                           >  4:234798/1‑275 (MQ=255)
                                                                   cATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGc                                                                                                                                                    <  2:234799/275‑1 (MQ=255)
                                                                      aCCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGc                                                                                                                                                                                                     >  2:338343/1‑223 (MQ=255)
                                                                          tATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAgg                                                                                                                                                                                                 <  4:338342/223‑1 (MQ=255)
                                                                           aTGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGc                                                                                                                                                                                                     >  2:67437/1‑218 (MQ=255)
                                                                               ccATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAgg                                                                                                                                                                                                 <  4:67436/218‑1 (MQ=255)
                                                                                cATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATa                                                                                                                                          <  2:224107/272‑1 (MQ=255)
                                                                                     tatCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCg                                                                                                                                                                                                                                                                 <  7:170275/148‑1 (MQ=255)
                                                                                        cTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGACTTTTACCgaga                                                                                                                                                                                                                                                              >  5:450034/1‑148 (MQ=255)
                                                                                        cTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGACTTTTACCgaga                                                                                                                                                                                                                                                              >  5:249236/1‑148 (MQ=255)
                                                                                            ttGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAAttt                                                                                                                                                                                                                                                                      >  1:136869/1‑136 (MQ=255)
                                                                                                cAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTAcc                                                                                                                                                                                                                                                                  <  3:136869/136‑1 (MQ=255)
                                                                                                  aaCAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAAtttt                                                                                                                                                                                                                                                                     >  7:37006/1‑131 (MQ=255)
                                                                                                  aaCAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAAtttt                                                                                                                                                                                                                                                                     <  8:37006/131‑1 (MQ=255)
                                                                                                  aaCAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTg                                                                                                                                                                                                                                                     <  8:861026/147‑1 (MQ=255)
                                                                                                  aaCAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTg                                                                                                                                                                                                                                                     >  7:861026/1‑147 (MQ=255)
                                                                                                      aTTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGaa                                                                                                                                                                                                                                                >  8:32226/1‑148 (MQ=255)
                                                                                                         aGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGACTTTTACCGAGAAGCAAACTGTTGAACCg                                                                                                                                                                                                                                             <  6:450034/148‑1 (MQ=255)
                                                                                                         aGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGACTTTTACCGAGAAGCAAACTGTTGAACCg                                                                                                                                                                                                                                             <  6:249236/148‑1 (MQ=255)
                                                                                                                gCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGt                                                                                                                                                                                                                                                    <  6:501304/134‑1 (MQ=255)
                                                                                                                gCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGt                                                                                                                                                                                                                                                    >  5:501304/1‑134 (MQ=255)
                                                                                                                gCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGc                                                                                                                                                                                                     >  4:90248/1‑181 (MQ=255)
                                                                                                                    gCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAgg                                                                                                                                                                                                 <  2:90249/181‑1 (MQ=255)
                                                                                                                        aCCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGaa                                                                                                                                                                                                                                                <  7:549962/130‑1 (MQ=255)
                                                                                                                        aCCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGaa                                                                                                                                                                                                                                                >  8:549962/1‑130 (MQ=255)
                                                                                                                            aCAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCaa                                                                                                                                                                                                                          >  8:534719/1‑148 (MQ=255)
                                                                                                                               gaagaaTCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGc                                                                                                                                                                                                                       <  7:32226/148‑1 (MQ=255)
                                                                                                                                     tCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTg                                                                                                                                                                                                                 <  7:534719/148‑1 (MQ=255)
                                                                                                                                            tGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAACCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGTAAAATACTACATCTGTTCGGGTGACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTaa                                                                           >  2:414989/1‑275 (MQ=255)
                                                                                                                                                            ccTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTg                                                                                                                                                                                          <  8:215937/148‑1 (MQ=255)
                                                                                                                                                            ccTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTg                                                                                                                                                                                          >  7:215937/1‑148 (MQ=255)
                                                                                                                                                            ccTTGATGCAGAAAATCGTTGTCCACTGGAAGCAGTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGAAAACTGTTGAACAGAGGCCGCAACCACAGTGAATGCCCGTTGACGACGGTGGGGCNGGCAGTCAGGGCCGCATTGAAATTATTTACTATATTTTTTTATATGCTGGTTTGATTGGTAAATTACTCCGTCGGTTGGGNGGCCAGTCTGCGTTCAGGTGGGGGGNACACGACCACAAATCTTTCCGGATATTTg                                                           >  3:89215/1‑275 (MQ=255)
                                                                                                                                                                   gCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAtggt                                                                                                                                                                                                         >  6:252351/1‑126 (MQ=255)
                                                                                                                                                                   gCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAtggt                                                                                                                                                                                                         <  5:252351/126‑1 (MQ=255)
                                                                                                                                                                    cAGATAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAAtttt                                                                                                                                                                                                                                                                     >  7:473917/1‑65 (MQ=255)
                                                                                                                                                                    cAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAAGCGGCGAAtttt                                                                                                                                                                                                                                                                     <  8:473917/65‑1 (MQ=255)
                                                                                                                                                                     aGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGc                                                                                                                                                                                                                                         <  5:591427/92‑1 (MQ=255)
                                                                                                                                                                     aGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGc                                                                                                                                                                                                                                         >  6:591427/1‑92 (MQ=255)
                                                                                                                                                                      gAAAATCGTTGTCCTCTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGcatca                                                                                                                                                                                   <  8:613146/145‑1 (MQ=255)
                                                                                                                                                                      gAAAATCGTTGTCCTCTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGcatca                                                                                                                                                                                   >  7:613146/1‑145 (MQ=255)
                                                                                                                                                                       aaaaTCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCt                                                                                                                                                                                           >  6:471108/1‑136 (MQ=255)
                                                                                                                                                                       aaaaTCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCt                                                                                                                                                                                           <  5:471108/136‑1 (MQ=255)
                                                                                                                                                                              ttGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGca                                                                                                                                                                                      >  1:78618/1‑134 (MQ=255)
                                                                                                                                                                                  ccACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGcatcat                                                                                                                                                                                  <  3:78618/134‑1 (MQ=255)
                                                                                                                                                                                           gCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAtggt                                                                                                                                                                                                         >  7:576173/1‑102 (MQ=255)
                                                                                                                                                                                           gCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAtggt                                                                                                                                                                                                         <  8:576173/102‑1 (MQ=255)
                                                                                                                                                                                               taattaTTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAGACTGTTGAACCGGAGCCGCAACCGCAGCCAATGcc                                                                                                                                                                                                                      >  8:763951/1‑85 (MQ=255)
                                                                                                                                                                                               taattaTTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGcc                                                                                                                                                                                                                      <  7:763951/85‑1 (MQ=255)
                                                                                                                                                                                                       cgtcgATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGGAAAATACTCCATCTGTTCGGGTGACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTAATCTTTCCGCATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCAt                >  2:508299/1‑275 (MQ=255)
                                                                                                                                                                                                           gATTATTAATTCAGGGGGCCAATTTTACCGGGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGGAAAATACTCCATCTGTTCGGGTGACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTAATCTTTCCGCATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATa               <  4:508298/272‑1 (MQ=255)
                                                                                                                                                                                                            attattAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAattatt                                                                                                                                                                        >  8:906077/1‑118 (MQ=255)
                                                                                                                                                                                                            attattAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAattatt                                                                                                                                                                        <  7:906077/118‑1 (MQ=255)
                                                                                                                                                                                                                ttaATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTAttt                                                                                                                                                                       <  6:833642/115‑1 (MQ=255)
                                                                                                                                                                                                                ttaATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTAttt                                                                                                                                                                       >  5:833642/1‑115 (MQ=255)
                                                                                                                                                                                                                 taATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAt                                                                                                                                                                                                            >  7:558393/1‑77 (MQ=255)
                                                                                                                                                                                                                 taATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGAt                                                                                                                                                                                                            <  8:558393/77‑1 (MQ=255)
                                                                                                                                                                                                                     tCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTg                                                                                                                                                                                          >  8:237578/1‑91 (MQ=255)
                                                                                                                                                                                                                     tCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTg                                                                                                                                                                                          <  7:237578/91‑1 (MQ=255)
                                                                                                                                                                                                                         gcggcgAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGTCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGGAAAATACTCCTTCTGTTCGGGTTACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTAATCTTTCCGCATTTTTTTAATGTTGTTTTTCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTCTACATTa  >  1:42601/1‑271 (MQ=255)
                                                                                                                                                                                                                               aaTTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAatt                                                                                                                                                                           <  6:904415/96‑1 (MQ=255)
                                                                                                                                                                                                                               aaTTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAatt                                                                                                                                                                           >  5:904415/1‑96 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                       
TTCACTGGTCGAGATGGTTTTGCACAAAGGACGTAAGCCGCGGATCTCACTCACGCGTGATACAGAACATACCTTATGGCCATTATATCTAGTTGCCAAACAATTAGTGCAGGCCTGCCAACCTACAGAAGAATCATTAATGCCGATTGTGAGCCGCCTTGATGCAGAAAATCGTTGTCCACTGGAAGCATTAATTATTCGTCGATTATTAATTCAGGCGGCGAATTTTACCGAGAAGCAAACTGTTGAACCGGAGCCGCAACCGCAGCCAATGCCCGTTGACGATGGTGGGCCAGGCTGTCTGGGCATCATTAAAATTATTTTCTATATTTTTATCTTTGCTGGTTTGATAGGGAAAATACTCCATCTGTTCGGGTGACAGTCGACGTTAATGTGCGGGAAGCCCGGCAAGTAATCTTTCCGCATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTA  >  NC_000913/678464‑678951

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