Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,847 C→A 5.0% A252S (GCA→TCA)  yfjH ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,8470CA 5.0% 285.9 / 2.2 89A252S (GCA→TCA) yfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base A (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTA  >  NC_000913/2757699‑2757995
                                                                                                                                                    |                                                                                                                                                    
tCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGc                                                                                                                                                      >  7:159859/1‑149 (MQ=255)
             aCCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCt                                                                                                                                                     <  7:219462/137‑1 (MQ=255)
             aCCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCt                                                                                                                                                     >  8:219462/1‑137 (MQ=255)
             aCCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTAttt                                                                                                                                         <  7:233899/149‑1 (MQ=255)
               cATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTAttttt                                                                                                                                       >  8:95463/1‑149 (MQ=255)
                aTTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTc                                                                                                                                      >  1:220709/1‑149 (MQ=255)
                    aTTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCAc                                                                                                                                   <  2:84785/148‑1 (MQ=255)
                     tttAAAATAGTATCCACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAAGACCATAGAGGTTCGCTAAGTTCATTAACTGCTACAACCTTATTTTTcgcccc                                                                                                                                 >  8:88587/1‑144 (MQ=255)
                       tAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCca                                                                                                                                <  2:225904/148‑1 (MQ=255)
                        aaaaTATTTTAAACAAAAGCATCTTTACTCATTGTTTATACACGATATTTTTTATTTTTATTA‑TTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCcatc                                                                                                                              <  5:102216/148‑1 (MQ=255)
                        aaaaTAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCcat                                                                                                                               <  3:220201/148‑1 (MQ=255)
                            tAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTc                                                                                                                          >  4:197754/1‑149 (MQ=255)
                               tATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCAGTAACTGCTACAACCTTATTTTTCCACCCATCATTCgg                                                                                                                        >  6:59707/1‑148 (MQ=255)
                                  cAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTTTCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGaa                                                                                                                      >  2:32723/1‑147 (MQ=255)
                                   aaCCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATTATTCGGaa                                                                                                                      >  4:122734/1‑146 (MQ=255)
                                   aaCCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGaa                                                                                                                      >  5:167550/1‑146 (MQ=255)
                                    aCCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACaga                                                                                                                  <  6:48265/149‑1 (MQ=255)
                                    aCCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACa                                                                                                                    >  5:122266/1‑147 (MQ=255)
                                    aCCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACa                                                                                                                    >  1:242434/1‑147 (MQ=255)
                                     ccAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACagag                                                                                                                 >  3:231612/1‑149 (MQ=255)
                                      cAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGt                                                                                                                <  2:220709/149‑1 (MQ=255)
                                        aaGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCa                                                                                                              >  2:159100/1‑149 (MQ=255)
                                         aGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGTGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCat                                                                                                             <  8:249636/149‑1 (MQ=255)
                                          gCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCatt                                                                                                            <  2:54568/149‑1 (MQ=255)
                                            aTCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCattaa                                                                                                          >  1:208734/1‑149 (MQ=255)
                                              cTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCattaat                                                                                                         >  6:88422/1‑148 (MQ=255)
                                                 gACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCattaatta                                                                                                       <  5:163126/147‑1 (MQ=255)
                                                 gACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCattaat                                                                                                         <  8:159859/145‑1 (MQ=255)
                                                 gACTCATTGTTGCTACACGATATTTTTTATTTTGATTAATTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCattaatta                                                                                                       <  5:163123/147‑1 (MQ=255)
                                                    tCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTc                                                                                                   >  6:139904/1‑148 (MQ=255)
                                                      aTTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCTATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTcaac                                                                                                <  2:185313/149‑1 (MQ=255)
                                                         gttgCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTcaaccaa                                                                                             >  1:69536/1‑149 (MQ=255)
                                                          ttgCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCAATCATTCGGAACAGAGTCATTAATTAGCTccaccaac                                                                                            >  7:57613/1‑149 (MQ=255)
                                                          ttgCTACACGATATTTTTTATTTTGATTAATTATCTCAAACTCTGCCATTTGAATAAAATCCCATAGATGTTCGATAAGTTCATTAACTGCTACAACCTGATTTTTCCAAACAACATTCGAAACATAGTCATTAATAAGATcacccacc                                                                                            >  2:286/1‑147 (MQ=255)
                                                              tACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGCCATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACtt                                                                                          >  6:22151/1‑147 (MQ=255)
                                                               acacGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAg                                                                                       >  3:9007/1‑149 (MQ=255)
                                                                cacGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCa                                                                                        <  5:96482/147‑1 (MQ=255)
                                                                cacGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCa                                                                                        <  1:159100/147‑1 (MQ=255)
                                                                cacGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCAATATCGTTTCGATAAGTTCATTAACTGATACAAACTTATTTTTCAACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTAAGt                                                                                      >  3:183977/1‑149 (MQ=255)
                                                                cacGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCACATAGCGGTTCGAGAAGTTCATTAACTGATACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGt                                                                                      >  8:110064/1‑149 (MQ=255)
                                                                cacGATATTTTTTATTTTCATTATTTATCTCAAACTCTCACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGATACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGt                                                                                      >  8:106665/1‑149 (MQ=255)
                                                                 acGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGtt                                                                                     <  7:143263/149‑1 (MQ=255)
                                                                 acGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGtt                                                                                     <  7:88587/149‑1 (MQ=255)
                                                                    atatTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa                                                                                  <  2:148302/149‑1 (MQ=255)
                                                                    atatTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa                                                                                  <  2:208734/149‑1 (MQ=255)
                                                                    atatTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGaa                                                                                  >  1:92768/1‑149 (MQ=255)
                                                                       ttttttATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATg                                                                                <  3:107337/148‑1 (MQ=255)
                                                                         ttttATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGat                                                                              <  3:209111/148‑1 (MQ=255)
                                                                            tATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATg                                                                           <  8:216890/148‑1 (MQ=255)
                                                                             aTTTTTATTAATTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGGCATTAATTAGCTCAACCAACTTCAGTTGAATGATATgg                                                                          <  8:112030/148‑1 (MQ=255)
                                                                                ttGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTac                                                                       <  7:95463/148‑1 (MQ=255)
                                                                                 tGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCTATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTaca                                                                      <  5:49295/148‑1 (MQ=255)
                                                                                   attattTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATAAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACtt                                                                   <  5:197864/149‑1 (MQ=255)
                                                                                           tctcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgcttt                                                            <  6:122266/148‑1 (MQ=255)
                                                                                           tctcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgcttt                                                            <  6:167550/148‑1 (MQ=255)
                                                                                           tctcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgcttt                                                            <  6:74859/148‑1 (MQ=255)
                                                                                             tcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTaca                                                                      >  7:43027/1‑136 (MQ=255)
                                                                                             tcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTaca                                                                      <  8:43027/136‑1 (MQ=255)
                                                                                             tcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAAAGATATGGTCCACTTCCgctttgct                                                         >  4:244/1‑149 (MQ=255)
                                                                                               aaaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgcttt                                                       >  8:35709/1‑149 (MQ=255)
                                                                                                aaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCtttt                                                      <  2:242434/149‑1 (MQ=255)
                                                                                                aaCTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCtttt                                                      <  7:120777/149‑1 (MQ=255)
                                                                                                   tctGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGcc                                                   >  8:85112/1‑149 (MQ=255)
                                                                                                         aTTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCctt                                             >  5:219805/1‑149 (MQ=255)
                                                                                                              aaTAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTc                                        >  3:76881/1‑149 (MQ=255)
                                                                                                                tAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCt                                      <  2:92768/149‑1 (MQ=255)
                                                                                                                 aaaaTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTg                                     >  5:94307/1‑149 (MQ=255)
                                                                                                                  aatttccATAGCGGTTCGATAAGTTGATTAACTGCTACAACCTTATTTTTCCTCCCATCTTTCGGTACAGAGTCATTAATTATCTCAACCTACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGc                                    <  8:215088/144‑1 (MQ=255)
                                                                                                                  aaaTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTg                                     >  7:239184/1‑148 (MQ=255)
                                                                                                                   aaTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCa                                   >  8:179527/1‑149 (MQ=255)
                                                                                                                          tAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTcc                                       <  1:160468/138‑1 (MQ=255)
                                                                                                                          tAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTcc                                       >  2:160468/1‑138 (MQ=255)
                                                                                                                               gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCctt                                             >  8:111281/1‑127 (MQ=255)
                                                                                                                               gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCctt                                             <  7:111281/127‑1 (MQ=255)
                                                                                                                               gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACt                       <  4:183977/149‑1 (MQ=255)
                                                                                                                                  cGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTc                     <  5:88422/148‑1 (MQ=255)
                                                                                                                                  cGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTc                     <  8:239184/148‑1 (MQ=255)
                                                                                                                                    aTAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGa                  <  4:231612/149‑1 (MQ=255)
                                                                                                                                     tAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGa                  >  7:33338/1‑148 (MQ=255)
                                                                                                                                         ttCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtt             >  1:98183/1‑149 (MQ=255)
                                                                                                                                           cATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttt           >  3:39901/1‑149 (MQ=255)
                                                                                                                                           cATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttt           >  3:138524/1‑149 (MQ=255)
                                                                                                                                            aTTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAttttt          <  7:35709/149‑1 (MQ=255)
                                                                                                                                                aCTGATACAACCTTATTTTTCCCCCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttttt       <  7:106665/148‑1 (MQ=255)
                                                                                                                                                aCTGATACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttttt       <  7:110064/148‑1 (MQ=255)
                                                                                                                                                 cTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGa     <  4:9007/149‑1 (MQ=255)
                                                                                                                                                 cTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGa     <  5:139904/149‑1 (MQ=255)
                                                                                                                                                  tGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttt         >  8:233657/1‑144 (MQ=255)
                                                                                                                                                  tGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttt         <  7:233657/144‑1 (MQ=255)
                                                                                                                                                  tGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGAt    >  2:65313/1‑149 (MQ=255)
                                                                                                                                                    cTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGAtt   <  4:138524/148‑1 (MQ=255)
                                                                                                                                                    cTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGAtt   <  6:94307/148‑1 (MQ=255)
                                                                                                                                                    cTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTa  <  7:179527/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTATCTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTA  >  NC_000913/2757699‑2757995

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