Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,705 C→T 8.4% intergenic (+42/+18) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,7050CT 8.4% 282.7 / 9.0 96intergenic (+42/+18)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base C (44/42);  new base T (4/4);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTC  >  NC_000913/681561‑681848
                                                                                                                                                |                                                                                                                                               
cgaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGAAAGGGTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTATTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATc                                                                                                                                            <  1:1586225/149‑1 (MQ=255)
 gaGGGTGACAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATcc                                                                                                                                           <  2:475665/150‑1 (MQ=255)
  ccggTGAAAAACCACAGAGCATTACCAGCATAATCAAACTCAGTTTAACGGTGCTCGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCg                                                                                                                                          <  1:681685/148‑1 (MQ=255)
   gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                                                                         <  2:983543/150‑1 (MQ=255)
   gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                                                                         <  1:793356/150‑1 (MQ=255)
   gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                                                                         >  2:339548/1‑150 (MQ=255)
   gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                                                                         >  1:1698763/1‑150 (MQ=255)
   gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                                                                         >  2:707778/1‑150 (MQ=255)
    ggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGc                                                                                                                                        <  1:1593271/150‑1 (MQ=255)
        aaaaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGcc                                                                                                                                    >  1:540349/1‑150 (MQ=255)
         aaaaCCACAGAGCAATACCACCATAATCAAAATCAGTTTAAGGGTGCTGGTGATAGGATATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCatcgg                                                                                                                                   <  2:1103609/150‑5 (MQ=38)
          aaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgt                                                                                                                                  >  2:1723961/1‑150 (MQ=255)
           aaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                                                                                 <  2:1132991/150‑1 (MQ=255)
           aaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                                                                                 <  2:1486223/150‑1 (MQ=255)
           aaCCACAGAGCATAACCAGCATTATCAAAATCTGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAACTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCGGGCAGCCgtg                                                                                                                                 <  2:495710/150‑1 (MQ=255)
           aaCCACACTGCATTACCAGCATTATCAAAATCAGTTTAACGGAGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                                                                                 <  2:844162/150‑1 (MQ=255)
            aCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGCGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                                                                                <  1:143630/150‑1 (MQ=255)
             caacaTAGCACTACCAGCAATAGCAAAAACAGTTTAACGCTGCTGGTGATAGGTTATGCTCTCGTCAAAATCGCGATGTTGGTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCACCCGGCAGCCGTGtt                                                                                                                               <  1:852369/148‑1 (MQ=255)
               acaGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                                                                             >  2:493790/1‑150 (MQ=255)
                  gagCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCTCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCAGTGTTTATgg                                                                                                                          >  2:874428/1‑150 (MQ=255)
                  gagCAATACCAGCATTATCAAAATCAGTTTAACGGTGCGGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCAGCAGCCGTGTTTATgg                                                                                                                          <  2:1246797/150‑1 (MQ=255)
                   cgCATTACCAGCATTATCAAAATACGTTTAACGGTGCTGGTGATAGGTAATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                                                                         <  2:1702304/149‑1 (MQ=255)
                       ttACCAGCATTATCAAAATCATTTTAACGGTGCTGGTGATAGGTTATGCTCTCGTCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtct                                                                                                                     <  2:1630570/150‑1 (MQ=255)
                        tACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtctc                                                                                                                    >  2:462374/1‑150 (MQ=255)
                         aCCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                                                                                   <  2:1526356/150‑1 (MQ=255)
                           cAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                                 >  2:258811/1‑150 (MQ=255)
                            aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                <  1:976996/150‑1 (MQ=255)
                            aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                <  2:645443/150‑1 (MQ=255)
                            aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                <  2:168470/150‑1 (MQ=255)
                            aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTAATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                <  1:1126250/150‑1 (MQ=255)
                                ctATCAAACTAAGATTAACGGTGCGGCGGATAGGTTATTCTCTCGGCAAAATCGCGATGAGGTTTAGATGACGATGTGAGTCTTGCGCCGATGCCGGACGCAGCGTTTACGCCTCATCCGGCAGCCGTGTTTATCCGAGCTCAATGGCTa                                                                                                            <  2:100635/149‑1 (MQ=255)
                                 tATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaa                                                                                                           >  2:750531/1‑150 (MQ=255)
                                 tATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaa                                                                                                           >  2:744158/1‑150 (MQ=255)
                                  aTCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaaa                                                                                                          >  1:1121848/1‑150 (MQ=255)
                                   acAAAATAAGTTTAACGCTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGAAGTTGTTAAGCTTAACATCTGTTTTTTGCGCCGATGCCGGCTTCAGCGTTTACTCAGCATCCGTCTGACGTGATAATGGCATCTCAAAGCCAAaat                                                                                                         >  1:1520696/2‑150 (MQ=255)
                                    cAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaata                                                                                                        >  1:1704654/1‑150 (MQ=255)
                                        aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                                                                                    <  2:407001/150‑1 (MQ=255)
                                           aGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAAGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctc                                                                                                 >  1:540032/1‑150 (MQ=255)
                                            gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                                                                                >  1:1470002/1‑150 (MQ=255)
                                            gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTAGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCCGCCGTGTTTATGGGATCGCAAAGGCTAAAAAAtcgctcg                                                                                                >  1:850328/1‑150 (MQ=255)
                                                aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccc                                                                                            >  2:1576816/1‑149 (MQ=255)
                                                 aCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                                                                           >  1:110414/1‑150 (MQ=255)
                                                          gTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTATATAATCGCTCGCCCGCCGCCGAACt                                                                                  >  1:1642118/1‑150 (MQ=255)
                                                           tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGAAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                                                                                 >  2:1246655/1‑150 (MQ=255)
                                                            gATAGGTTATGCTCTCGGCTAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTc                                                                                <  1:1730552/150‑1 (MQ=255)
                                                            gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGAGTTTACGCCGCATCCGGCAGCAGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCAGCCGCCGCACtta                                                                                >  1:1440464/1‑149 (MQ=255)
                                                             aTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                                                                               <  2:471744/150‑1 (MQ=255)
                                                             aTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                                                                               <  2:427813/150‑1 (MQ=255)
                                                               aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                             <  1:808768/150‑1 (MQ=255)
                                                               aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                             >  2:1624835/1‑150 (MQ=255)
                                                                ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                            <  2:110414/150‑1 (MQ=255)
                                                                ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                            >  1:572350/1‑150 (MQ=255)
                                                                ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                            >  1:753263/1‑150 (MQ=255)
                                                                ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATAGCTCGCCCGCCGCCGAAGTTCATcc                                                                            >  2:136839/1‑150 (MQ=255)
                                                                   tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                             <  2:1180312/146‑1 (MQ=255)
                                                                   tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                             >  1:1180312/1‑146 (MQ=255)
                                                                     tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgcg                                                                       >  2:13190/1‑150 (MQ=255)
                                                                      gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGTTCTCAATGGCTAAATAATCGCTCGCCCGCCCCCGCACTTCATCCACGCgg                                                                      >  2:1039710/1‑150 (MQ=255)
                                                                        tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGc                                                                    <  1:258811/150‑1 (MQ=255)
                                                                           cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat                                                                 >  2:1652091/1‑150 (MQ=255)
                                                                              cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATc                                                              >  1:1548247/1‑150 (MQ=255)
                                                                                 aaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt                                                           <  2:469081/150‑1 (MQ=255)
                                                                                   tCGCGTTGTTGTTTAGCTGACGATGTGTTTTTTGCTCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGGAGCCTTGTTTATGGGTTCTCAATGGCTAAATAATGGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                         <  1:120199/150‑1 (MQ=255)
                                                                                   tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                         <  1:959651/150‑1 (MQ=255)
                                                                                   tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                         <  1:480234/150‑1 (MQ=255)
                                                                                   tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                         <  2:208061/150‑1 (MQ=255)
                                                                                    cgcaATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGa                                                        >  2:1649696/1‑150 (MQ=255)
                                                                                     gcgATGTTGTTTAGCTGACGCTGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGAc                                                       <  1:837374/150‑1 (MQ=255)
                                                                                            tgttTAGCTGCGGATGTGTTTTTTGCGCCGCTGCCGTATGCTGCGTTTACGCGGCATCCGGCAGCCGTGTTTATGGTATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcg                                                <  2:1520833/150‑1 (MQ=255)
                                                                                                tAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCCTCCGGCCGCCGGGTTTATGGGATCGCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCCCGCGCGCCATATCGTTTGCCGGCAAcgccgc                                            >  2:58312/1‑150 (MQ=255)
                                                                                                 aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCa                                           >  2:1197154/1‑150 (MQ=255)
                                                                                                   cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATc                                         >  2:1580176/1‑150 (MQ=255)
                                                                                                     gACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCaa                                       >  2:704897/1‑150 (MQ=255)
                                                                                                     gACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCaa                                       >  2:531054/1‑150 (MQ=255)
                                                                                                     gACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTAATGGGATCACAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCCCGCGGGCCATATCGTTTGACTGCAACGCCGCATCaa                                       >  2:1539901/1‑150 (MQ=255)
                                                                                                         aTGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATc                                   <  1:13190/150‑1 (MQ=255)
                                                                                                             gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                               <  2:1631565/150‑1 (MQ=255)
                                                                                                             gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGACGCCGAACTACATCCCCGCGGGCCATATCGTTTGACTGCACCGGCCCATCAAAATCGGTg                               >  2:323841/1‑150 (MQ=255)
                                                                                                              ttttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGa                              <  2:1698763/150‑1 (MQ=255)
                                                                                                                        ggctGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGcc                    >  1:1058402/4‑150 (MQ=255)
                                                                                                                         gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGCCTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCg                   >  2:1282540/1‑150 (MQ=255)
                                                                                                                          aTGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCgc                  >  1:542993/1‑149 (MQ=255)
                                                                                                                            gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtc                <  2:618870/150‑1 (MQ=255)
                                                                                                                             ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtct               >  1:1720001/1‑150 (MQ=255)
                                                                                                                             ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtct               <  1:1723961/150‑1 (MQ=255)
                                                                                                                               ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct             >  2:523584/1‑150 (MQ=255)
                                                                                                                               ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct             >  2:3263/1‑150 (MQ=255)
                                                                                                                               ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct             >  1:698530/1‑150 (MQ=255)
                                                                                                                               ggATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct             <  2:572350/150‑1 (MQ=255)
                                                                                                                                   gCAGCGTTTACGCCGCATCCGTCAGCCGGGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCCCGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTtctc         <  2:237777/150‑1 (MQ=255)
                                                                                                                                    cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcg        >  1:402311/1‑150 (MQ=255)
                                                                                                                                    cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcg        >  2:183072/1‑150 (MQ=255)
                                                                                                                                      gCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgct      <  1:1651710/150‑1 (MQ=255)
                                                                                                                                      gCGTTTACGCCGCAGCGGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgct      <  1:1334022/150‑1 (MQ=255)
                                                                                                                                        gTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgctcg    <  2:1345705/150‑1 (MQ=255)
                                                                                                                                          ttACGCCGCATCCGGCAGCCGGGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTGCCTCCCCGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGGGAACACCCGGCCAAGCTCTTCTCGCTCGTc  >  2:881777/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTC  >  NC_000913/681561‑681848

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