Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 356,794 C→T 5.7% Y208Y (TAC→TAT codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913356,7940CT 5.7% 282.7 / 4.3 88Y208Y (TAC→TATcodAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base C (38/44);  new base T (3/2);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGC  >  NC_000913/356646‑356940
                                                                                                                                                    |                                                                                                                                                  
ggATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATAc                                                                                                                                                    <  5:47044/149‑1 (MQ=255)
ggATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATAc                                                                                                                                                    <  3:82188/149‑1 (MQ=255)
 gATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACg                                                                                                                                                   >  6:101256/1‑149 (MQ=255)
  aTTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGa                                                                                                                                                  <  3:350825/149‑1 (MQ=255)
    tttGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTTGAGTCGCTGCATACACCCTTCGCCCTGGCGCAAAAATActa                                                                                                                                                  >  6:87835/1‑145 (MQ=255)
    tttGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGCGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAATTCGACCg                                                                                                                                               >  5:248093/1‑149 (MQ=255)
       gTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGt                                                                                                                                              <  6:108392/148‑1 (MQ=255)
           tATTCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTAAt                                                                                                                                         >  1:311438/1‑149 (MQ=255)
               ccAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGa                                                                                                                                      <  5:123101/148‑1 (MQ=255)
               ccAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGAc                                                                                                                                     >  5:175239/1‑149 (MQ=255)
                cAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACg                                                                                                                                    >  8:153442/1‑149 (MQ=255)
                cAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACg                                                                                                                                    >  7:64273/1‑149 (MQ=255)
                   cGGTGAAGCGTTGCTGGAAGAGGCGTTTCGCTTAGGGGCAGATGGAGGGGGTCCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGGGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTc                                                                                                                                 <  6:52943/149‑1 (MQ=255)
                    ggTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCa                                                                                                                                <  3:227573/149‑1 (MQ=255)
                        aaGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACtgt                                                                                                                            >  4:142557/1‑149 (MQ=255)
                         aGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACtgtg                                                                                                                           <  7:116337/149‑1 (MQ=255)
                          gCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAATATACGACCGTCTCATCGACGTTCACTGtga                                                                                                                          <  5:318163/149‑1 (MQ=255)
                           cGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGtgat                                                                                                                         <  5:350649/149‑1 (MQ=255)
                             ttGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGtgat                                                                                                                         >  3:136737/1‑147 (MQ=255)
                               gCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATgag                                                                                                                      <  2:82712/148‑1 (MQ=255)
                                cTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATgaga                                                                                                                     <  1:38296/148‑1 (MQ=255)
                                 tGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGtgat                                                                                                                         >  8:326611/1‑143 (MQ=255)
                                 tGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGtgat                                                                                                                         >  3:352246/1‑143 (MQ=255)
                                  ggAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATc                                                                                                                   >  2:340220/1‑148 (MQ=255)
                                    aaGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATATGACCGTCTCATCGACGTTCACTGtga                                                                                                                          >  5:351803/1‑139 (MQ=255)
                                    aaGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATATGACCGTCTCATCGACGTTCACTGtga                                                                                                                          <  8:282062/139‑1 (MQ=255)
                                    aaGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATATGACCGTCTCATCGACGTTCACTGtga                                                                                                                          <  6:351803/139‑1 (MQ=255)
                                    aaGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATATGACCGTCTCATCGACGTTCACTGtga                                                                                                                          >  7:282062/1‑139 (MQ=255)
                                     agagGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGCCCGTCTCATCGACGTTCACTGTGATGAGATCGATg                                                                                                               <  3:142557/149‑1 (MQ=255)
                                      gagGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGa                                                                                                              <  1:82329/149‑1 (MQ=255)
                                       agGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGAc                                                                                                             <  7:153442/149‑1 (MQ=255)
                                        ggCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGAc                                                                                                             >  6:24491/1‑148 (MQ=255)
                                        ggCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACg                                                                                                            <  8:8090/149‑1 (MQ=255)
                                         gCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGa                                                                                                           <  7:82596/149‑1 (MQ=255)
                                                gCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGa                                                                                                              <  4:127239/139‑1 (MQ=255)
                                                gCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGa                                                                                                              >  3:127239/1‑139 (MQ=255)
                                                  ttAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGAt                                                                                                                <  8:156405/135‑1 (MQ=255)
                                                  ttAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGAt                                                                                                                >  7:156405/1‑135 (MQ=255)
                                                  ttAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTcgcg                                                                                                  >  8:343433/1‑149 (MQ=255)
                                                    aGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCTTAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTcgc                                                                                                   <  5:87835/146‑1 (MQ=255)
                                                     ggggCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCtt                                                                                               <  7:343433/149‑1 (MQ=255)
                                                        gCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGAcc                                                                                                                                                >  6:207959/1‑97 (MQ=255)
                                                        gCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGAcc                                                                                                                                                <  5:207959/97‑1 (MQ=255)
                                                        gCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGt                                                                                            <  6:175239/149‑1 (MQ=255)
                                                         cAGATGTTGTGTTGGCGATTGCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTTGCTGCATATAACCTTCGCCCTGGAGGAAAAATAGGAGGGTCTCATCGACGTTCACTGTGTTGAGATCGTTGTCGTGCAGTAGCGCTTTGTc                                                                                           <  6:150709/149‑1 (MQ=255)
                                                             tGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGaaa                                                                                       <  8:34196/149‑1 (MQ=255)
                                                               tAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAAcc                                                                                     >  6:201818/1‑149 (MQ=255)
                                                                  tGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGtt                                                                                  >  4:232179/1‑149 (MQ=255)
                                                                    ggggCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGa                                                                                                              >  2:235584/1‑119 (MQ=255)
                                                                    ggggCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGa                                                                                                              <  1:235584/119‑1 (MQ=255)
                                                                    ggggCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTtgc                                                                                >  2:15542/1‑149 (MQ=255)
                                                                      ggCGATTCCGCATTTTGAATTTACCCGTGAACACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTtgctg                                                                              >  5:255992/1‑149 (MQ=255)
                                                                       gCGATTCCGCATTTTGTATTTACCCGGTAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTc                                                                                           <  8:354417/135‑1 (MQ=255)
                                                                       gCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTc                                                                                           >  7:354417/1‑135 (MQ=255)
                                                                          aTTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCt                                                                          <  7:246673/149‑1 (MQ=255)
                                                                          aTTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCt                                                                          <  1:15542/149‑1 (MQ=255)
                                                                               gCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTtg                                                                                 >  3:347120/1‑137 (MQ=255)
                                                                               gCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTtg                                                                                 <  4:347120/137‑1 (MQ=255)
                                                                                 aTTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCAc                                                                   <  7:326611/149‑1 (MQ=255)
                                                                                  ttttGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCAc                                                                   >  8:325009/1‑148 (MQ=255)
                                                                                  ttttGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCAc                                                                   >  5:114067/1‑148 (MQ=255)
                                                                                       aaTTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCTCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGa                                                              <  2:311438/148‑1 (MQ=255)
                                                                                              ccGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATgg                                                      >  5:108765/1‑149 (MQ=255)
                                                                                              ccGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATgg                                                      >  1:157587/1‑149 (MQ=255)
                                                                                              ccGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATgg                                                      <  1:406293/149‑1 (MQ=255)
                                                                                                gTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGgc                                                    >  1:335547/1‑149 (MQ=255)
                                                                                                       cGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGa                                              <  1:340220/148‑1 (MQ=255)
                                                                                                       cGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGa                                              <  7:325009/148‑1 (MQ=255)
                                                                                                        ggCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGt                                            <  6:255992/149‑1 (MQ=255)
                                                                                                            tGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCAcc                                        <  2:335547/149‑1 (MQ=255)
                                                                                                                gTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCgcca                                    <  5:15447/149‑1 (MQ=255)
                                                                                                                gTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCgc                                      >  6:15447/1‑147 (MQ=255)
                                                                                                                       cATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCacac                             <  5:24491/149‑1 (MQ=255)
                                                                                                                          aaaaCCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCAcaccac                          >  7:67249/1‑149 (MQ=255)
                                                                                                                           aaaCCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACg                         >  4:128448/1‑149 (MQ=255)
                                                                                                                            aaCCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACgg                        >  3:373527/1‑149 (MQ=255)
                                                                                                                               cTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCaa                     >  1:153609/1‑149 (MQ=255)
                                                                                                                               cTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGGCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCaa                     <  5:101256/149‑1 (MQ=255)
                                                                                                                                    cccTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCAc                >  1:64566/1‑149 (MQ=255)
                                                                                                                                       tGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCacac                             <  3:58676/133‑1 (MQ=255)
                                                                                                                                       tGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCacac                             >  4:58676/1‑133 (MQ=255)
                                                                                                                                          cgcAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATg                   >  6:178386/1‑140 (MQ=255)
                                                                                                                                          cgcAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATg                   <  5:178386/140‑1 (MQ=255)
                                                                                                                                           gcAAAAATATGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCCTGGGCGCGCGAGTCACCGCCAGCCACACCCCGGCAATGCACTCCtat          >  1:248834/1‑148 (MQ=255)
                                                                                                                                           gcAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCtat          >  1:154218/1‑148 (MQ=255)
                                                                                                                                               aaaTACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACgg     <  8:67249/149‑1 (MQ=255)
                                                                                                                                                aaTACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGGAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACggg    <  6:114067/149‑1 (MQ=255)
                                                                                                                                                 aTACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAAt                    <  5:380977/132‑1 (MQ=255)
                                                                                                                                                 aTACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAAt                    >  6:380977/1‑132 (MQ=255)
                                                                                                                                                  tACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGc  <  3:128448/149‑1 (MQ=255)
                                                                                                                                                  tACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGc  <  2:64566/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGC  >  NC_000913/356646‑356940

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