Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,961 C→T 5.6% A214T (GCA→ACA)  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,9610CT 5.6% 347.7 / 6.0 106A214T (GCA→ACA) yfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (52/48);  new base T (3/3);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 106 total aligned reads displayed.

AATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCA  >  NC_000913/2757814‑2758108
                                                                                                                                                   |                                                                                                                                                    
aaTCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCa                                                                                                                                                     >  1:555255/1‑149 (MQ=255)
  tCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACt                                                                                                                                                   >  2:286498/1‑149 (MQ=255)
    ccATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTcc                                                                                                                                                 >  1:472067/1‑149 (MQ=255)
      aTAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGc                                                                                                                                               >  1:162180/1‑149 (MQ=255)
      aTAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGc                                                                                                                                               >  2:321484/1‑149 (MQ=255)
          cGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGa                                                                                                                                           <  1:176216/149‑1 (MQ=255)
            gTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACt                                                                                                                                         >  4:60910/1‑149 (MQ=255)
              tctaTAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACCCTTCCTCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACGTc                                                                                                                                       <  4:90020/146‑1 (MQ=255)
               cGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCt                                                                                                                                      <  1:511926/149‑1 (MQ=255)
                    aGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATaa                                                                                                                                 <  2:91060/149‑1 (MQ=255)
                    aGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATaa                                                                                                                                 <  1:187357/149‑1 (MQ=255)
                    aGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTAGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACACCGCCTGACTTCTGATaa                                                                                                                                 <  2:38272/149‑1 (MQ=255)
                     gTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAt                                                                                                                               >  2:421510/1‑149 (MQ=255)
                      ttCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtt                                                                                                                              >  2:438161/1‑149 (MQ=255)
                      ttCATTAACTGCTACAACCTTATTTTTCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAttt                                                                                                                             >  4:358157/1‑149 (MQ=255)
                        cATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttt                                                                                                                            <  3:241896/149‑1 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttt                                                                                                                          >  3:285303/1‑149 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttt                                                                                                                          <  2:428622/149‑1 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttt                                                                                                                          >  1:99797/1‑149 (MQ=255)
                            aaCTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAAtttttttt                                                                                                                        <  1:183704/149‑1 (MQ=255)
                                    cAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaa                                                                                                                >  4:35904/1‑149 (MQ=255)
                                     aaCCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               >  2:117182/1‑149 (MQ=255)
                                        cTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACag                                                                                                            >  3:73427/1‑149 (MQ=255)
                                        cTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACag                                                                                                            >  1:381086/1‑149 (MQ=255)
                                          tATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACagag                                                                                                          >  2:302848/1‑149 (MQ=255)
                                                  cACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTaaa                                                                                                  <  4:226827/149‑1 (MQ=255)
                                                  cACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTaaa                                                                                                  <  2:110707/149‑1 (MQ=255)
                                                   aCCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAAc                                                                                                 >  2:35165/1‑149 (MQ=255)
                                                   aCCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCCTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAAc                                                                                                 >  3:229726/1‑149 (MQ=255)
                                                    cccATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAAcc                                                                                                <  4:388734/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  3:60910/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  1:35165/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  2:472067/149‑1 (MQ=255)
                                                       atcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAg                                                                                             <  1:348866/149‑1 (MQ=255)
                                                       atcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAg                                                                                             <  2:162180/149‑1 (MQ=255)
                                                          atTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTTGTTTACAAACAGAGATGTTAAACCAAGTc                                                                                           >  3:523558/1‑149 (MQ=255)
                                                          atTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTc                                                                                           >  3:230536/1‑149 (MQ=255)
                                                              ggAACAGAGTCATTAATTAGCCCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTAGCCACCCTTCTTTCCTGCACTCCGCCTGCCTTCTGATAATTTTTTTTGATTACAAACAGAAATGTTAAACCAAGACATAcc                                                                                      >  2:297433/1‑149 (MQ=255)
                                                                aaCAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTa                                                                                    >  2:521575/1‑149 (MQ=255)
                                                                aaCAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTa                                                                                    >  4:455156/1‑149 (MQ=255)
                                                                    gagTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCa                                                                                <  4:511828/149‑1 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAgaga                                                                                                         >  3:315985/1‑123 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAgaga                                                                                                         <  4:315985/123‑1 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCAc                                                                               >  4:70978/1‑149 (MQ=255)
                                                                       tCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACa                                                                                                             >  2:465061/1‑117 (MQ=255)
                                                                       tCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACa                                                                                                             <  1:465061/117‑1 (MQ=255)
                                                                         attaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTAt                                                                           >  3:474850/1‑149 (MQ=255)
                                                                          ttaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATc                                                                          <  1:216453/149‑1 (MQ=255)
                                                                          ttaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATc                                                                          <  3:358157/149‑1 (MQ=255)
                                                                            aattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGa                                                                                                                      <  1:178551/103‑1 (MQ=255)
                                                                            aattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTACCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGa                                                                                                                      >  2:178551/1‑103 (MQ=255)
                                                                             attaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  2:381086/149‑1 (MQ=255)
                                                                             attaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTCTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  2:495225/149‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               <  4:295429/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               >  1:344160/1‑108 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               <  4:513562/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               >  3:513562/1‑108 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               <  2:344160/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAAt                                                                      >  4:217075/1‑149 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTGTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAATTTTTTTTGATTACaaa                                                                                                               >  3:295429/1‑108 (MQ=255)
                                                                               taGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATa                                                                     <  1:117182/149‑1 (MQ=255)
                                                                                aGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATaa                                                                    <  4:73427/149‑1 (MQ=255)
                                                                                  cTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACt                                                                  <  1:29360/149‑1 (MQ=255)
                                                                                       ccaacTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCa                                                             >  2:51902/1‑149 (MQ=255)
                                                                                         aacTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGc                                                           >  1:504315/1‑149 (MQ=255)
                                                                                            ttCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCg                                                        >  2:569906/1‑149 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       >  4:379953/1‑142 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       >  4:145744/1‑142 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       <  3:379953/142‑1 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       <  3:145744/142‑1 (MQ=255)
                                                                                                       gATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             <  1:75939/149‑1 (MQ=255)
                                                                                                        atatGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCttt                                            >  2:165676/1‑149 (MQ=255)
                                                                                                            ggTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAg                                        <  1:50475/149‑1 (MQ=255)
                                                                                                            ggTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAg                                        <  1:514239/149‑1 (MQ=255)
                                                                                                              tACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGc                                      <  4:509400/149‑1 (MQ=255)
                                                                                                              tACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGc                                      <  2:555255/149‑1 (MQ=255)
                                                                                                                cacTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCaa                                    <  3:76974/149‑1 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:470912/1‑149 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:122858/1‑149 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:137798/1‑149 (MQ=255)
                                                                                                                  cTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             >  3:108066/1‑138 (MQ=255)
                                                                                                                  cTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             <  4:108066/138‑1 (MQ=255)
                                                                                                                    tCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGtt                                <  3:84727/149‑1 (MQ=255)
                                                                                                                       gctttgctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttc                             >  1:569424/1‑149 (MQ=255)
                                                                                                                           tgctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattt                         <  3:57865/149‑1 (MQ=255)
                                                                                                                             ctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattttc                       <  1:21218/149‑1 (MQ=255)
                                                                                                                              ttttGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattttca                      >  4:406513/1‑149 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGCCTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGa               >  2:563616/1‑149 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       >  2:387118/1‑93 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  1:387118/93‑1 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGa               <  1:417391/149‑1 (MQ=255)
                                                                                                                                      cccTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGAt              <  1:302848/149‑1 (MQ=255)
                                                                                                                                          tcttTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAt          <  2:532280/149‑1 (MQ=255)
                                                                                                                                          tcttTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAt          <  2:189554/149‑1 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAGGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  3:443623/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  4:219245/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  4:197661/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  1:471579/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCACGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  2:482538/1‑149 (MQ=255)
                                                                                                                                              tCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTc      <  1:424546/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCA  >  NC_000913/2757814‑2758108

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