Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,639 A→C S105A (TCA→GCA)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,6390AC96.7% 333.8 / 1.1 90S105A (TCA→GCA) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (37/50);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-02

GCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCC  >  NC_000913/4485492‑4485784
                                                                                                                                                   |                                                                                                                                                 
gctgacgacGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGcc                                                                                                                                                <  1:354372/148‑1 (MQ=255)
gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGcc                                                                                                                                                >  1:427929/1‑151 (MQ=255)
   gacgacGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGt                                                                                                                                             <  2:974372/151‑1 (MQ=255)
    acgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGt                                                                                                                                             >  2:1143020/1‑150 (MQ=255)
       acgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCa                                                                                                                                         <  2:183846/151‑1 (MQ=255)
       acgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCa                                                                                                                                         <  1:454617/151‑1 (MQ=255)
         gcgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGAACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCAtt                                                                                                                                       <  2:875674/151‑1 (MQ=255)
            gTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCTTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAg                                                                                                                                    <  2:1183175/151‑1 (MQ=255)
            gTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAg                                                                                                                                    <  2:813785/151‑1 (MQ=255)
                gCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCg                                                                                                                                  >  1:78542/1‑149 (MQ=255)
                gCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGt                                                                                                                                 >  1:759583/1‑150 (MQ=255)
                     tGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTaa                                                                                                                            <  2:529205/150‑1 (MQ=255)
                         cGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGtt                                                                                                                       <  2:82733/151‑1 (MQ=255)
                         cGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGtt                                                                                                                       <  1:874514/151‑1 (MQ=255)
                         cGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGt                                                                                                                        <  1:262494/150‑1 (MQ=255)
                         cGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGt                                                                                                                        <  2:1218656/150‑1 (MQ=255)
                            tttCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGtttt                                                                                                                     >  2:288135/1‑150 (MQ=255)
                             ttCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTc                                                                                                                    <  1:1234291/150‑1 (MQ=255)
                               cAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTg                                                                                                                  <  1:912148/150‑1 (MQ=255)
                               cAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTg                                                                                                                  >  1:1288963/1‑150 (MQ=255)
                               cAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGACTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCt                                                                                                                   >  2:1257838/1‑149 (MQ=255)
                                aGCTGATCGAAACTGTAGAGCGTATCTTTTGCCGTCTCGTCAGCCTGCCGCGCTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGaa                                                                                                                <  2:304464/151‑1 (MQ=255)
                                  cggATCGAAACTGAAGAGCGACCCTTTCGCCGTCTCGACAGCTTGACGCACTTTCAAGAAGTTGTTACGGCCTTTAACGTGCAGCTCAGACAGAAAGCAGACCGCCCCTATCGCTTCATTATCATTCAGCGTATTAATGGGTTTCt                                                                                                                   >  1:749959/3‑146 (MQ=255)
                                    gATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAAc                                                                                                            <  1:1026182/151‑1 (MQ=255)
                                    gATCGAAACTGTAGAGCGTCTCTTTTGCCGGCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAAc                                                                                                            <  2:725375/151‑1 (MQ=255)
                                     aTCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATa                                                                                                              <  1:218157/148‑1 (MQ=255)
                                      tCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCt                                                                                                          >  1:903949/1‑151 (MQ=255)
                                       cGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTg                                                                                                         <  1:1032066/151‑1 (MQ=255)
                                        gAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGc                                                                                                        <  1:1189556/151‑1 (MQ=255)
                                        gAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGc                                                                                                        >  2:790980/1‑151 (MQ=255)
                                           aCTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTg                                                                                                     <  2:514221/151‑1 (MQ=255)
                                             tGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGt                                                                                                    >  1:459219/1‑150 (MQ=255)
                                              gTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTAc                                                                                                  <  2:383455/151‑1 (MQ=255)
                                                    cGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGGGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGc                                                                                            >  2:270272/1‑151 (MQ=255)
                                                      tctcTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTc                                                                                          >  1:1172570/1‑151 (MQ=255)
                                                       ctctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTc                                                                                          >  2:591529/1‑150 (MQ=255)
                                                        tctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCAt                                                                                        >  2:334476/1‑151 (MQ=255)
                                                         ctTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTc                                                                                          <  2:52493/148‑1 (MQ=255)
                                                           tttGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCAt                                                                                        <  1:1199601/148‑1 (MQ=255)
                                                                      gACAGACTGACGCACTTTCAAATAGTTGTGACGGTCTTTAACGTGCAGCTCAGTGAGAAAGCAGACCGCTTCCATTGCGCCAGAATCATTCAGCGAATAAATGGTTATCTGAATAACCTGCTTGTACTGACGCTCAACCAGCTCACGtt                                                                            <  1:550092/149‑1 (MQ=255)
                                                                       aCAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCtt                                                                         <  1:1185756/151‑1 (MQ=255)
                                                                       aCAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCtt                                                                         >  2:882272/1‑151 (MQ=255)
                                                                        cAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCTCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCttt                                                                        <  1:187478/151‑1 (MQ=255)
                                                                          gCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGc                                                                      <  2:461690/151‑1 (MQ=255)
                                                                           ccTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGcc                                                                     >  1:261007/1‑151 (MQ=255)
                                                                           ccTCACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGcc                                                                     >  1:1058335/1‑151 (MQ=255)
                                                                             tGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGc                                                                   >  1:530834/1‑151 (MQ=255)
                                                                              gACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTGCTGACGCTCATCCAGCTCACGTTCTTTGCCGc                                                                   <  2:181816/150‑1 (MQ=255)
                                                                              gACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCg                                                                    <  1:766002/149‑1 (MQ=255)
                                                                              gACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGc                                                                   <  2:505604/150‑1 (MQ=255)
                                                                                  cACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAAcc                                                               >  1:459431/1‑150 (MQ=255)
                                                                                   aCTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCaa                                                             <  2:373454/151‑1 (MQ=255)
                                                                                      ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                           >  2:976026/1‑150 (MQ=255)
                                                                                      ttCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                           <  1:676066/150‑1 (MQ=255)
                                                                                       tCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACTGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATa                                                           >  1:278872/1‑149 (MQ=255)
                                                                                         cAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGga                                                       <  1:1261720/151‑1 (MQ=255)
                                                                                           gtagtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAg                                                         >  1:547044/1‑147 (MQ=255)
                                                                                             agtTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAAt                                                   <  1:484791/151‑1 (MQ=255)
                                                                                                tgttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTc                                                 <  2:868138/150‑1 (MQ=255)
                                                                                                tgttACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCg                                                >  2:1181405/1‑151 (MQ=255)
                                                                                                     cGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                           >  1:910451/1‑151 (MQ=255)
                                                                                                      ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                           >  2:1053097/1‑150 (MQ=255)
                                                                                                      ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                           >  1:321474/1‑150 (MQ=255)
                                                                                                      ggCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCgg                                           >  1:31767/1‑150 (MQ=255)
                                                                                                        ccTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGaaa                                        <  2:663598/151‑1 (MQ=255)
                                                                                                          tttAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGt                                      >  2:914454/1‑151 (MQ=255)
                                                                                                            tAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTAc                                    >  1:151681/1‑151 (MQ=255)
                                                                                                             aaCGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTAc                                    <  2:770656/150‑1 (MQ=255)
                                                                                                              aCGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACAt                                  <  2:909309/151‑1 (MQ=255)
                                                                                                                gTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTc                                <  1:237833/151‑1 (MQ=255)
                                                                                                                gTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTAACATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTc                                >  1:1211652/1‑151 (MQ=255)
                                                                                                                 tGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCg                               >  1:819624/1‑151 (MQ=255)
                                                                                                                    aGCTCAGTCAGAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGa                             >  1:902443/1‑150 (MQ=255)
                                                                                                                        cagtcagAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGTGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATg                        <  2:173724/151‑1 (MQ=255)
                                                                                                                              gAAAGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgc                    >  2:768069/1‑149 (MQ=255)
                                                                                                                               aaaGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcagc                 >  2:588075/1‑151 (MQ=255)
                                                                                                                               aaaGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcagc                 <  1:405883/151‑1 (MQ=255)
                                                                                                                               aaaGCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTgcagc                 <  1:717926/151‑1 (MQ=255)
                                                                                                                                  gCAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGcaac              <  2:594526/151‑1 (MQ=255)
                                                                                                                                   cAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGcaaca             <  1:1142957/151‑1 (MQ=255)
                                                                                                                                   cAGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGcaaca             <  1:100898/151‑1 (MQ=255)
                                                                                                                                    aGACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGcaaca             >  2:846731/1‑150 (MQ=255)
                                                                                                                                     gACCGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAAt           >  2:1147205/1‑151 (MQ=255)
                                                                                                                                        cGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTc        >  1:287222/1‑151 (MQ=255)
                                                                                                                                        cGCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTc        <  2:451552/151‑1 (MQ=255)
                                                                                                                                         gCTTCCATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCt       <  1:347301/151‑1 (MQ=255)
                                                                                                                                         gCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCt       <  2:903985/151‑1 (MQ=255)
                                                                                                                                             ccATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTg      <  2:1084462/148‑1 (MQ=255)
                                                                                                                                             ccATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGccc   <  2:920659/151‑1 (MQ=255)
                                                                                                                                              cATTGCGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGcccc  >  1:869126/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACGGCCTTTAACGTGCAGCTCAGTCAGAAAGCAGACCGCTTCCATTGAGCCAGTATCATTCAGCGTATTAATGGTTTTCTGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCC  >  NC_000913/4485492‑4485784

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