Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 554,198 A→G R80R (CGT→CGC)  ppiB ← peptidyl‑prolyl cis‑trans isomerase B (rotamase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913554,1980AG99.2% 485.1 / ‑4.5 133R80R (CGT→CGC) ppiBpeptidyl‑prolyl cis‑trans isomerase B (rotamase B)
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (58/74);  total (59/74)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

Only 100 of 133 total aligned reads displayed.

GCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCACGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAACCTTCGCG  >  NC_000913/554051‑554344
                                                                                                                                                   |                                                                                                                                                  
gCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCgg                                                                                                                                                 <  1:153389/151‑1 (MQ=255)
  cGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGt                                                                                                                                               >  3:165143/1‑151 (MQ=255)
  cGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGt                                                                                                                                               >  1:360783/1‑151 (MQ=255)
  cGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGt                                                                                                                                               <  1:165129/151‑1 (MQ=255)
   gTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTa                                                                                                                                              <  4:1026341/151‑1 (MQ=255)
    tCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTAt                                                                                                                                             >  2:479902/1‑151 (MQ=255)
     cAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTAtt                                                                                                                                            >  2:1567631/1‑151 (MQ=255)
      aaCCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCGGAGTACGTGCCATTGCCAGCGTACCGCGGGTAttt                                                                                                                                           <  1:446625/151‑1 (MQ=255)
        ccACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTAttttt                                                                                                                                         <  4:1298764/151‑1 (MQ=255)
        ccACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTAttttt                                                                                                                                         >  1:329237/1‑151 (MQ=255)
        ccACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTAttttt                                                                                                                                         <  4:792468/151‑1 (MQ=255)
         cACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTc                                                                                                                                        <  1:58188/151‑1 (MQ=255)
          aCTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCa                                                                                                                                       <  4:1369905/151‑1 (MQ=255)
          aCTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCa                                                                                                                                       <  1:467998/151‑1 (MQ=255)
               aGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCg                                                                                                                                  <  2:402549/151‑1 (MQ=255)
                 cAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgtt                                                                                                                                <  4:361432/151‑1 (MQ=255)
                 cAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgtt                                                                                                                                <  2:480490/151‑1 (MQ=255)
                 cAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgtt                                                                                                                                <  1:227498/151‑1 (MQ=255)
                  aaaCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgttg                                                                                                                               <  4:24803/151‑1 (MQ=255)
                  aaaCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgttg                                                                                                                               <  1:335962/151‑1 (MQ=255)
                  aaaCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgttg                                                                                                                               <  2:854228/151‑1 (MQ=255)
                   aaCACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCgttgt                                                                                                                              >  4:549627/1‑151 (MQ=255)
                        gCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCt                                                                                                                         >  4:901848/1‑151 (MQ=255)
                           gTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCg                                                                                                                      >  3:4292/1‑151 (MQ=255)
                            tAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGt                                                                                                                     >  1:145111/1‑151 (MQ=255)
                             aGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGtt                                                                                                                    <  1:40609/151‑1 (MQ=255)
                                cccAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGttttt                                                                                                                 >  3:104789/1‑151 (MQ=255)
                                 ccAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTg                                                                                                                <  4:1095256/151‑1 (MQ=255)
                                   aaCCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGAt                                                                                                              >  3:332721/1‑151 (MQ=255)
                                    aCCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATc                                                                                                             >  1:281909/1‑151 (MQ=255)
                                      cTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCgg                                                                                                           >  3:142334/1‑151 (MQ=255)
                                          cAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCt                                                                                                       >  2:735697/1‑151 (MQ=255)
                                           aGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCtt                                                                                                      <  4:477727/151‑1 (MQ=255)
                                            ggCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCttt                                                                                                     <  4:367474/151‑1 (MQ=255)
                                            ggCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCttt                                                                                                     <  2:516804/151‑1 (MQ=255)
                                              cTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGGTGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTtg                                                                                                    >  1:339605/1‑150 (MQ=255)
                                                 tCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTtggtg                                                                                                 <  4:1300142/150‑1 (MQ=255)
                                                 tCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTtggtg                                                                                                 <  3:174516/150‑1 (MQ=255)
                                                   gCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCt                                                                                              >  2:1622920/1‑151 (MQ=255)
                                                     cAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCttt                                                                                            <  4:346309/151‑1 (MQ=255)
                                                     cAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCttt                                                                                            >  2:860201/1‑151 (MQ=255)
                                                       gagaAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCttttt                                                                                          <  1:474566/151‑1 (MQ=255)
                                                         gaAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGt                                                                                        >  2:1565869/1‑151 (MQ=255)
                                                         gaAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGt                                                                                        <  4:1184025/151‑1 (MQ=255)
                                                           aGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGttt                                                                                      >  4:241455/1‑151 (MQ=255)
                                                           aGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGttt                                                                                      <  3:196375/151‑1 (MQ=255)
                                                            gTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTc                                                                                     <  1:388690/151‑1 (MQ=255)
                                                              tCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCAt                                                                                   >  4:499481/1‑151 (MQ=255)
                                                               cAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCAt                                                                                   <  2:1027313/150‑1 (MQ=255)
                                                               cAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATg                                                                                  <  2:1027277/151‑1 (MQ=255)
                                                                aGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGc                                                                                 >  4:14141/1‑151 (MQ=255)
                                                                  gAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGccc                                                                               >  4:862686/1‑151 (MQ=255)
                                                                        gTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCa                                                                         <  2:498781/151‑1 (MQ=255)
                                                                          tATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCaa                                                                        >  4:970389/1‑150 (MQ=255)
                                                                           aTCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCAGTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCaaa                                                                       <  4:57785/150‑1 (MQ=255)
                                                                            tCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAc                                                                     >  2:81448/1‑151 (MQ=255)
                                                                            tCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAc                                                                     >  3:215309/1‑151 (MQ=255)
                                                                             cAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAcc                                                                    >  4:1436723/1‑151 (MQ=255)
                                                                             cAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAcc                                                                    <  4:851655/151‑1 (MQ=255)
                                                                             cAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAcc                                                                    <  4:732516/151‑1 (MQ=255)
                                                                             cAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAcc                                                                    <  3:360800/151‑1 (MQ=255)
                                                                                 cACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAAccgccg                                                                >  2:767969/1‑151 (MQ=255)
                                                                                  aCGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACcgccgc                                                               >  1:37953/1‑151 (MQ=255)
                                                                                      tgatgaAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCt                                                            <  2:1282844/150‑1 (MQ=255)
                                                                                           aagaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATc                                                       >  4:267137/1‑150 (MQ=255)
                                                                                           aagaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATc                                                       >  1:30649/1‑150 (MQ=255)
                                                                                            agaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCAt                                                     <  2:406590/151‑1 (MQ=255)
                                                                                            agaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCACGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCa                                                      >  2:360279/1‑150 (MQ=255)
                                                                                             gaaCTGTGCAGTTGCATAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATa                                                    <  4:1511151/151‑1 (MQ=255)
                                                                                             gaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCAt                                                     >  3:273365/1‑150 (MQ=255)
                                                                                             gaaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATa                                                    <  3:381563/151‑1 (MQ=255)
                                                                                              aaCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATa                                                    >  4:1430670/1‑150 (MQ=255)
                                                                                               aCTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATaa                                                   >  4:1103070/1‑150 (MQ=255)
                                                                                                cTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATaa                                                   >  2:136536/1‑149 (MQ=255)
                                                                                                  gtgCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAg                                                 <  3:145125/149‑1 (MQ=255)
                                                                                                    gCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCg                                              <  2:772921/150‑1 (MQ=255)
                                                                                                    gCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGt                                             <  2:1117183/151‑1 (MQ=255)
                                                                                                          gCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATaa                                       <  4:540532/151‑1 (MQ=255)
                                                                                                          gCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATa                                        >  4:1044829/1‑150 (MQ=255)
                                                                                                          gCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATa                                        <  1:271278/150‑1 (MQ=255)
                                                                                                           cAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATaa                                       >  2:1293908/1‑150 (MQ=255)
                                                                                                            tgagTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGAATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAac                                      >  2:377988/2‑150 (MQ=255)
                                                                                                            agagTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAaca                                     <  2:1286471/151‑1 (MQ=255)
                                                                                                            agagTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAac                                      >  2:775539/1‑150 (MQ=255)
                                                                                                              agTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACg                                   >  2:1450734/1‑151 (MQ=255)
                                                                                                                tGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGt                                 >  4:481732/1‑151 (MQ=255)
                                                                                                                  cGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTgg                               <  2:477872/151‑1 (MQ=255)
                                                                                                                  cGGAGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTgg                               <  4:1569577/151‑1 (MQ=255)
                                                                                                                     aGCCTGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGaaa                            <  4:132181/151‑1 (MQ=255)
                                                                                                                         tGAGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATgg                        >  4:1179446/1‑151 (MQ=255)
                                                                                                                           aGTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGgtg                      >  3:376815/1‑151 (MQ=255)
                                                                                                                            gTACGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGgtgt                     >  2:35980/1‑151 (MQ=255)
                                                                                                                              aCGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGtgttg                   <  4:1440971/151‑1 (MQ=255)
                                                                                                                              aCGTGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGtgttg                   <  4:1526093/151‑1 (MQ=255)
                                                                                                                                 tGCCATTGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGtgttgtt                 <  1:372225/150‑1 (MQ=255)
                                                                                                                                       tGCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTaaaa           >  2:805634/1‑150 (MQ=255)
                                                                                                                                        gCCAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAAc         >  2:1453306/1‑151 (MQ=255)
                                                                                                                                         ccAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAAcc        >  4:1475192/1‑151 (MQ=255)
                                                                                                                                          cAGCGTACCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAACCt       >  2:34638/1‑151 (MQ=255)
                                                                                                                                                aCCGCGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAACCTTcgca  <  2:900677/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGTCAACCACTTCAGCAAACACGCAGTAGCCCCAACCTTGCAGGCTTTCGCCAGAGAAGTTCAGGAAGTCGTTATCAACCACGTTGATGAAGAACTGTGCAGTTGCAGAGTGCGGAGCCTGAGTACGTGCCATTGCCAGCGTACCACGGGTATTTTTCAGGCCGTTGTTGGCTTCGTTTTTGATCGGTTCTTTGGTGGCTTTTTGTTTCATGCCCGGTTCAAAACCGCCGCCCTGAATCATAAAGCCGTTGATAACACGGTGGAAAATGGTGTTGTTGTAAAAACCTTCGCG  >  NC_000913/554051‑554344

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