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.

AATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCA  >  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)
                     gTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑t                                                                                                                               >  2:421510/1‑149 (MQ=255)
                      ttCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tt                                                                                                                              >  2:438161/1‑149 (MQ=255)
                      ttCATTAACTGCTACAACCTTATTTTT‑CACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑ttt                                                                                                                             >  4:358157/1‑149 (MQ=255)
                        cATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tttt                                                                                                                            <  3:241896/149‑1 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tttttt                                                                                                                          >  3:285303/1‑149 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tttttt                                                                                                                          <  2:428622/149‑1 (MQ=255)
                          ttAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tttttt                                                                                                                          >  1:99797/1‑149 (MQ=255)
                            aaCTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑tttttttt                                                                                                                        <  1:183704/149‑1 (MQ=255)
                                    cAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaa                                                                                                                >  4:35904/1‑149 (MQ=255)
                                     aaCCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               >  2:117182/1‑149 (MQ=255)
                                        cTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACag                                                                                                            >  3:73427/1‑149 (MQ=255)
                                        cTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACag                                                                                                            >  1:381086/1‑149 (MQ=255)
                                          tATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACagag                                                                                                          >  2:302848/1‑149 (MQ=255)
                                                  cACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTaaa                                                                                                  <  4:226827/149‑1 (MQ=255)
                                                  cACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTaaa                                                                                                  <  2:110707/149‑1 (MQ=255)
                                                   aCCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAAc                                                                                                 >  2:35165/1‑149 (MQ=255)
                                                   aCCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCCTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAAc                                                                                                 >  3:229726/1‑149 (MQ=255)
                                                    cccATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAAcc                                                                                                <  4:388734/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  3:60910/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  1:35165/149‑1 (MQ=255)
                                                      catcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCaa                                                                                              <  2:472067/149‑1 (MQ=255)
                                                       atcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAg                                                                                             <  1:348866/149‑1 (MQ=255)
                                                       atcatTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAg                                                                                             <  2:162180/149‑1 (MQ=255)
                                                          atTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTTGTTTACAAACAGAGATGTTAAACCAAGTc                                                                                           >  3:523558/1‑149 (MQ=255)
                                                          atTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAATTTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTc                                                                                           >  3:230536/1‑149 (MQ=255)
                                                              ggAACAGAGTCATTAATTAGCCCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTAGCCACCCTTCTTTCCTGCACTCCGCCTGCCTTCTGATAA‑TTTTTTTTGATTACAAACAGAAATGTTAAACCAAGACATAcc                                                                                      >  2:297433/1‑149 (MQ=255)
                                                                aaCAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTa                                                                                    >  2:521575/1‑149 (MQ=255)
                                                                aaCAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTa                                                                                    >  4:455156/1‑149 (MQ=255)
                                                                    gagTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCa                                                                                <  4:511828/149‑1 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAgaga                                                                                                         >  3:315985/1‑123 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAgaga                                                                                                         <  4:315985/123‑1 (MQ=255)
                                                                     agTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCAc                                                                               >  4:70978/1‑149 (MQ=255)
                                                                       tCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACa                                                                                                             >  2:465061/1‑117 (MQ=255)
                                                                       tCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACa                                                                                                             <  1:465061/117‑1 (MQ=255)
                                                                         attaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTAt                                                                           >  3:474850/1‑149 (MQ=255)
                                                                          ttaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATc                                                                          <  1:216453/149‑1 (MQ=255)
                                                                          ttaattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATc                                                                          <  3:358157/149‑1 (MQ=255)
                                                                            aattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGa                                                                                                                      <  1:178551/103‑1 (MQ=255)
                                                                            aattaGCTCAACCAACTTCAGTTGAATGATATGGTACACTACCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGa                                                                                                                      >  2:178551/1‑103 (MQ=255)
                                                                             attaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  2:381086/149‑1 (MQ=255)
                                                                             attaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTCTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  2:495225/149‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               <  4:295429/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               >  1:344160/1‑108 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               <  4:513562/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               >  3:513562/1‑108 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               <  2:344160/108‑1 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAAt                                                                      >  4:217075/1‑149 (MQ=255)
                                                                              ttaGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTGTGCCACCCTTCTTTCCTGTACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACaaa                                                                                                               >  3:295429/1‑108 (MQ=255)
                                                                               taGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATa                                                                     <  1:117182/149‑1 (MQ=255)
                                                                                aGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATaa                                                                    <  4:73427/149‑1 (MQ=255)
                                                                                  cTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACt                                                                  <  1:29360/149‑1 (MQ=255)
                                                                                       ccaacTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCa                                                             >  2:51902/1‑149 (MQ=255)
                                                                                         aacTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGc                                                           >  1:504315/1‑149 (MQ=255)
                                                                                            ttCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCg                                                        >  2:569906/1‑149 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       >  4:379953/1‑142 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       >  4:145744/1‑142 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       <  3:379953/142‑1 (MQ=255)
                                                                                                    aaTGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGa                                                       <  3:145744/142‑1 (MQ=255)
                                                                                                       gATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             <  1:75939/149‑1 (MQ=255)
                                                                                                        atatGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCttt                                            >  2:165676/1‑149 (MQ=255)
                                                                                                            ggTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAg                                        <  1:50475/149‑1 (MQ=255)
                                                                                                            ggTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAg                                        <  1:514239/149‑1 (MQ=255)
                                                                                                              tACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGc                                      <  4:509400/149‑1 (MQ=255)
                                                                                                              tACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGc                                      <  2:555255/149‑1 (MQ=255)
                                                                                                                cacTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCaa                                    <  3:76974/149‑1 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:470912/1‑149 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:122858/1‑149 (MQ=255)
                                                                                                                 acTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAAc                                   >  4:137798/1‑149 (MQ=255)
                                                                                                                  cTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             >  3:108066/1‑138 (MQ=255)
                                                                                                                  cTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCtt                                             <  4:108066/138‑1 (MQ=255)
                                                                                                                    tCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGtt                                <  3:84727/149‑1 (MQ=255)
                                                                                                                       gctttgctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttc                             >  1:569424/1‑149 (MQ=255)
                                                                                                                           tgctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattt                         <  3:57865/149‑1 (MQ=255)
                                                                                                                             ctttTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattttc                       <  1:21218/149‑1 (MQ=255)
                                                                                                                              ttttGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGttttcattttca                      >  4:406513/1‑149 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGCCTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGa               >  2:563616/1‑149 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       >  2:387118/1‑93 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCaaa                                                                       <  1:387118/93‑1 (MQ=255)
                                                                                                                                     aCCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGa               <  1:417391/149‑1 (MQ=255)
                                                                                                                                      cccTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGAt              <  1:302848/149‑1 (MQ=255)
                                                                                                                                          tcttTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAt          <  2:532280/149‑1 (MQ=255)
                                                                                                                                          tcttTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAt          <  2:189554/149‑1 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAGGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  3:443623/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  4:219245/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  4:197661/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  1:471579/1‑149 (MQ=255)
                                                                                                                                            tttCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCACGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGAttt        >  2:482538/1‑149 (MQ=255)
                                                                                                                                              tCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTc      <  1:424546/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTCTTTCCTGCACTCCGCCTGACTTCTGATAA‑TTTTTTTTGATTACAAACAGAGATGTTAAACCAAGTCATACCTATGCACTTATCAAATAACTCCGCAGCCCGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCA  >  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: ‑