Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 739,616 C→A 6.5% A37D (GCT→GAT)  phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913739,6160CA 6.5% 548.0 / 20.6 153A37D (GCT→GAT) phrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (73/70);  new base A (5/5);  total (78/75)
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.85e-01

Only 100 of 159 total aligned reads displayed.

GGCCGCGTGTGATTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGC  >  NC_000913/739470‑739754
                                                                                                                                                  |                                                                                                                                          
ggCCGCGTGTGATTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTa                                                                                                                                          >  3:216482/1‑149 (MQ=255)
   cgcgTGTGATTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTacac                                                                                                                                       >  1:1210306/1‑149 (MQ=255)
      gtgtGATTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTAcaccac                                                                                                                                    <  3:1145258/149‑1 (MQ=255)
           aTTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAg                                                                                                                               >  4:462410/1‑149 (MQ=255)
             tAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTg                                                                                                                             <  3:462410/149‑1 (MQ=255)
                    cgcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACg                                                                                                                      <  3:578580/149‑1 (MQ=255)
                     gcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGc                                                                                                                     >  3:1048159/1‑149 (MQ=255)
                     gcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGc                                                                                                                     <  4:1048159/149‑1 (MQ=255)
                            aGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGc                                                                                                                                            <  2:1177600/119‑1 (MQ=255)
                             ggAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATaa                                                                                                                 >  3:1005872/1‑145 (MQ=255)
                             ggAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATg                                                                                                             >  1:542665/1‑149 (MQ=255)
                             ggAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGATACACCACGCCAGTGGGCGACGCATAACATg                                                                                                             >  3:734761/1‑149 (MQ=255)
                             ggAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGATACACCACGCCAGTGGGCGACGCATAACATg                                                                                                             >  1:1138030/1‑149 (MQ=255)
                                   ttATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTcgccg                                                                                                       >  2:98132/1‑149 (MQ=255)
                                       gACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTc                                                                                                   >  1:278502/1‑149 (MQ=255)
                                       gACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTc                                                                                                   >  1:653222/1‑149 (MQ=255)
                                        aCTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCa                                                                                                  >  1:1344212/1‑149 (MQ=255)
                                          tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGt                                                                                                                              >  2:1465310/1‑119 (MQ=255)
                                          tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGt                                                                                                                              <  1:1465310/119‑1 (MQ=255)
                                          tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTcg                                                                                                          >  1:645060/1‑139 (MQ=255)
                                          tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTcg                                                                                                          <  2:645060/139‑1 (MQ=255)
                                          tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAgg                                                                                                >  1:952576/1‑149 (MQ=255)
                                            cccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCt                                                                                              <  2:542665/149‑1 (MQ=255)
                                             ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTg                                                                                                                             >  3:984069/1‑117 (MQ=255)
                                               atgtggtctggtTTCGCCCGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCTAGTGGGCGACGCATAACa                                                                                                               >  4:483475/4‑129 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGc                                                                                                                                            >  1:930983/1‑100 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGc                                                                                                                                            <  2:930983/100‑1 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGc                                                                                                                                            <  4:637498/100‑1 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACa                                                                                                               <  3:483475/129‑1 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           <  1:613518/149‑1 (MQ=255)
                                               aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           <  2:278502/149‑1 (MQ=255)
                                                  tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGATACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACtt                                                                                        <  2:1138030/149‑1 (MQ=255)
                                                  tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGATACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACtt                                                                                        <  4:734761/149‑1 (MQ=255)
                                                    gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTc                                                                                                           <  3:1407210/128‑1 (MQ=255)
                                                    gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           >  2:16694/1‑144 (MQ=255)
                                                    gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           <  1:16694/144‑1 (MQ=255)
                                                    gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTAt                                                                                      <  2:863033/149‑1 (MQ=255)
                                                     tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATc                                                                                     >  1:1158424/1‑149 (MQ=255)
                                                     tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATc                                                                                     >  3:522544/1‑149 (MQ=255)
                                                       tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTAcaccac                                                                                                                                    >  2:457705/1‑100 (MQ=255)
                                                       tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTAcaccac                                                                                                                                    <  1:457705/100‑1 (MQ=255)
                                                               ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCAt                                                                                                                   >  2:678678/1‑109 (MQ=255)
                                                               ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAAc                                                                           <  2:1344212/149‑1 (MQ=255)
                                                               ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAAc                                                                           >  2:685537/1‑149 (MQ=255)
                                                                   gATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGa                                                                                            >  2:669613/1‑128 (MQ=255)
                                                                   gATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGa                                                                                            <  1:669613/128‑1 (MQ=255)
                                                                   gATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGaa                                                                       <  4:216482/149‑1 (MQ=255)
                                                                      ttACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATgg                                                                    >  1:600742/1‑149 (MQ=255)
                                                                           tCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATc                                                                                     <  4:115993/127‑1 (MQ=255)
                                                                             tGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACaa                                                             >  1:1002409/1‑149 (MQ=255)
                                                                             tGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTGGTGTAGCGATATACAACGCCa                                                                                                                                >  4:1269268/1‑82 (MQ=25)
                                                                             tGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTGGTGTAGCGATATACAACGCCa                                                                                                                                <  3:1269268/82‑1 (MQ=25)
                                                                               cACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTgg                                                                                                                            >  3:541955/1‑84 (MQ=255)
                                                                               cACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTgg                                                                                                                            <  4:541955/84‑1 (MQ=255)
                                                                                aCGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATa                                                          >  4:169156/1‑149 (MQ=255)
                                                                                aCGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATa                                                          <  3:169156/149‑1 (MQ=255)
                                                                                 cGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAg                                                         <  1:999544/149‑1 (MQ=255)
                                                                                  gATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGt                                                                                                                              >  3:1334093/1‑79 (MQ=255)
                                                                                  gATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGt                                                                                                                              <  4:1334093/79‑1 (MQ=255)
                                                                                         tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACg                                                                                                                      <  4:297808/80‑1 (MQ=255)
                                                                                         tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACg                                                                                                                      >  3:297808/1‑80 (MQ=255)
                                                                                         tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAAc                                                                                          <  3:585577/108‑1 (MQ=255)
                                                                                         tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCg                                                 >  3:448207/1‑149 (MQ=255)
                                                                                            cACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGaa                                              >  4:1121537/1‑149 (MQ=255)
                                                                                               tGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGaaaaa                                           <  2:653222/149‑1 (MQ=255)
                                                                                                ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAg                                          <  3:1192314/149‑1 (MQ=255)
                                                                                                ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAg                                          <  3:1121537/149‑1 (MQ=255)
                                                                                                  ctgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGt                                        <  1:685537/149‑1 (MQ=255)
                                                                                                   tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCGGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATc                                                                                     >  1:1310738/1‑103 (MQ=255)
                                                                                                   tgccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCGGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATc                                                                                     <  2:1310738/103‑1 (MQ=255)
                                                                                                    gccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTAt                                      >  1:527654/1‑149 (MQ=255)
                                                                                                     ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTAtt                                     <  4:1368819/149‑1 (MQ=255)
                                                                                                         ctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCtt                                 <  4:284197/149‑1 (MQ=255)
                                                                                                                      tCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGa                    >  2:440086/1‑149 (MQ=255)
                                                                                                                        tGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTa                                       <  1:1015165/128‑1 (MQ=255)
                                                                                                                        tGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAg                  >  1:347962/1‑149 (MQ=255)
                                                                                                                           aCGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAg                  >  1:67668/1‑146 (MQ=255)
                                                                                                                            cgcgTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGa              <  2:1002409/149‑1 (MQ=255)
                                                                                                                             gcgTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCaa                                                                            >  4:1235375/1‑86 (MQ=255)
                                                                                                                             gcgTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGAt             >  2:1093224/1‑149 (MQ=255)
                                                                                                                               gTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGa           <  4:448207/149‑1 (MQ=255)
                                                                                                                                   tGGCGTGGTGTAGCGATATACAACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTa                               >  4:197926/1‑125 (MQ=255)
                                                                                                                                   tGGCGTGGTGTAGAGATATACAACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTa                               <  3:197926/125‑1 (MQ=255)
                                                                                                                                       gTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGt                      <  4:867731/130‑1 (MQ=255)
                                                                                                                                       gTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGt                      >  3:867731/1‑130 (MQ=255)
                                                                                                                                        ttGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTAt                                                                                      <  3:285479/65‑1 (MQ=255)
                                                                                                                                        ttGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGc  >  2:1351045/1‑149 (MQ=255)
                                                                                                                                        ttGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGc  >  2:1035842/1‑149 (MQ=255)
                                                                                                                                        ttGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGc  >  1:372046/1‑149 (MQ=255)
                                                                                                                                        ttGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGc  <  1:1351045/149‑1 (MQ=255)
                                                                                                                                           tataTCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATg                                                                     <  4:797023/79‑1 (MQ=255)
                                                                                                                                           tataTCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATg                                                                     >  3:797023/1‑79 (MQ=255)
                                                                                                                                           tataTCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGa              <  4:488846/134‑1 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           <  3:706380/54‑1 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGaa                                                                                           >  4:706380/1‑54 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTCTCAATGCTCAACTGAATggg                                                                   >  2:1189239/1‑78 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATggg                                                                   >  2:546451/1‑78 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAg                                          >  4:1440859/1‑103 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAg                                          <  3:1440859/103‑1 (MQ=255)
                                                                                                                                              atCGCTACACCACGCCAGAGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTTAACTGAATggg                                                                   <  1:1189239/78‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GGCCGCGTGTGATTAACTTGCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTCGC  >  NC_000913/739470‑739754

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