Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,920,591 T→C T135A (ACC→GCC)  b1838 ← protein phosphatase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,920,5910TC80.9% 440.3 / 103.0 199T135A (ACC→GCC) b1838protein phosphatase 1
Reads supporting (aligned to +/- strand):  ref base T (23/15);  new base C (85/76);  total (108/91)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 206 total aligned reads displayed.

GGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAAT  >  NC_000913/1920501‑1920681
                                                                                          |                                                                                          
ggCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGgcgc                                                                                          >  2:5116867/1‑93 (MQ=255)
ggCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGc                                                                                          <  1:2466958/93‑1 (MQ=255)
 gCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGcgcg                                                                                         >  2:3458957/1‑93 (MQ=255)
 gCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGcgcg                                                                                         >  2:3090030/1‑93 (MQ=255)
 gCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCg                                                                                         <  1:5037287/93‑1 (MQ=255)
  cTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGa                                                                                        >  2:3055119/1‑93 (MQ=255)
    ccACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACt                                                                                      >  2:2575736/1‑93 (MQ=255)
     cacaAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACtg                                                                                     <  1:3836687/93‑1 (MQ=255)
     cacaAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACtg                                                                                     <  1:2314622/93‑1 (MQ=255)
     cacaAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACtg                                                                                     <  1:1806991/93‑1 (MQ=255)
      acaAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCTGCGCGACtgt                                                                                    <  1:5821349/93‑1 (MQ=255)
       caAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACtgtg                                                                                   >  1:4175942/1‑93 (MQ=255)
        aaGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTgggt                                                                                  >  2:6645456/1‑93 (MQ=255)
        aaGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTgtgt                                                                                  <  2:4753183/93‑1 (MQ=255)
        aaGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTgtgt                                                                                  <  2:4026333/93‑1 (MQ=255)
         aGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTa                                                                                 >  1:6587272/1‑93 (MQ=255)
         aGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTa                                                                                 <  1:787336/93‑1 (MQ=255)
           aCCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACt                                                                               >  1:6320289/1‑93 (MQ=255)
             cTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTc                                                                             >  1:4231462/1‑93 (MQ=255)
             cTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTc                                                                             >  2:198456/1‑93 (MQ=255)
             cTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTc                                                                             >  1:1509428/1‑93 (MQ=255)
              tGATGCAAATCAACGTCCTTTTGCCATGCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCa                                                                            <  2:791866/93‑1 (MQ=255)
               gATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCaa                                                                           >  1:3859980/1‑93 (MQ=255)
               gATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCaa                                                                           >  1:4305617/1‑93 (MQ=255)
               gATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCaa                                                                           >  1:5075503/1‑93 (MQ=255)
                aTGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAg                                                                          >  2:1577812/1‑93 (MQ=255)
                   cAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAAt                                                                       >  2:4063589/1‑93 (MQ=255)
                   cAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAAt                                                                       <  1:4533028/93‑1 (MQ=255)
                    aaaTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATa                                                                      <  1:620425/93‑1 (MQ=255)
                     aaTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATaa                                                                     >  1:3954963/1‑93 (MQ=255)
                     aaTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATaa                                                                     >  2:4471879/1‑93 (MQ=255)
                     aaTCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATaa                                                                     <  2:5602507/93‑1 (MQ=255)
                       tCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAg                                                                   >  1:3208400/1‑93 (MQ=255)
                         aaCGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAggg                                                                 >  1:3180258/1‑93 (MQ=255)
                         aaCGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAggg                                                                 >  2:589509/1‑93 (MQ=255)
                          aCGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGc                                                                <  1:3782090/93‑1 (MQ=255)
                             tCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCaaa                                                             <  1:5591152/93‑1 (MQ=255)
                              ccTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAAt                                                            >  1:5392335/1‑93 (MQ=255)
                              ccTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAAt                                                            >  1:1448392/1‑93 (MQ=255)
                              ccTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAAt                                                            <  1:1384431/93‑1 (MQ=255)
                               cTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATg                                                           <  2:1279472/93‑1 (MQ=255)
                                ttttGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGCGCAAATGc                                                          >  2:1360907/1‑93 (MQ=255)
                                 tttGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCt                                                         >  2:1863375/1‑93 (MQ=255)
                                 tttGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCt                                                         <  2:3315451/93‑1 (MQ=255)
                                 tttGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCt                                                         <  1:3259581/93‑1 (MQ=255)
                                  ttGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTg                                                        >  1:1335548/1‑93 (MQ=255)
                                    gCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGAc                                                      >  1:4061618/1‑93 (MQ=255)
                                      cATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAt                                                    >  1:175994/1‑93 (MQ=255)
                                      cATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAt                                                    <  2:6438487/93‑1 (MQ=255)
                                      cATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCCGACAt                                                    >  2:3700463/1‑93 (MQ=255)
                                        ttCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAttt                                                  <  1:5913480/93‑1 (MQ=255)
                                        ttCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAttt                                                  >  2:1049542/1‑93 (MQ=255)
                                          cATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACAttttt                                                >  1:2375945/1‑93 (MQ=255)
                                             aaaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCa                                             <  2:2835333/93‑1 (MQ=255)
                                             aaaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCa                                             <  1:3117676/93‑1 (MQ=255)
                                              aaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAg                                            >  2:3536169/1‑93 (MQ=255)
                                              aaCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAg                                            <  1:956811/93‑1 (MQ=255)
                                               aCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgc                                           <  1:1943453/93‑1 (MQ=255)
                                               aCATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgc                                           <  2:4542879/93‑1 (MQ=255)
                                                catcatCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcg                                          <  2:2345448/93‑1 (MQ=255)
                                                catcatCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcg                                          <  1:1286312/93‑1 (MQ=255)
                                                catcatCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcg                                          >  1:5460376/1‑93 (MQ=255)
                                                 atcatcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcgc                                         >  2:210255/1‑93 (MQ=255)
                                                 atcatcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcgc                                         >  1:4680228/1‑93 (MQ=255)
                                                 atcatcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAgcgc                                         <  1:572140/93‑1 (MQ=255)
                                                    atcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGt                                      >  2:110896/1‑93 (MQ=255)
                                                    atcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGt                                      >  2:6535600/1‑93 (MQ=255)
                                                    atcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGt                                      >  1:6063737/1‑93 (MQ=255)
                                                     tcTGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGtt                                     >  2:4984967/1‑93 (MQ=255)
                                                       tGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGtttt                                   >  2:5742841/1‑93 (MQ=255)
                                                       tGGATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGtttt                                   >  1:6236234/1‑93 (MQ=255)
                                                        ggATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTc                                  <  2:265620/93‑1 (MQ=255)
                                                        ggATAATCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTc                                  <  1:1120963/93‑1 (MQ=255)
                                                          aTAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGc                                >  1:2022746/1‑93 (MQ=255)
                                                             aTCGGCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTg                             >  1:1650084/1‑93 (MQ=255)
                                                                ggCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGttt                          <  1:3424016/93‑1 (MQ=255)
                                                                ggCATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGttt                          <  1:4330299/93‑1 (MQ=255)
                                                                  cATGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCt                        <  2:3416884/93‑1 (MQ=255)
                                                                    tGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGt                      >  2:3600777/1‑93 (MQ=255)
                                                                    tGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGt                      >  1:5024045/1‑93 (MQ=255)
                                                                    tGAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGt                      <  2:4386419/93‑1 (MQ=255)
                                                                     gAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGtt                     >  2:5202687/1‑93 (MQ=255)
                                                                     gAGCAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGtt                     <  2:4212759/93‑1 (MQ=255)
                                                                        cAATAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGAt                  >  2:6603172/1‑93 (MQ=255)
                                                                          aTAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTa                <  2:4178645/93‑1 (MQ=255)
                                                                          aTAACATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTa                <  1:136722/93‑1 (MQ=255)
                                                                            aaCATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc              >  2:6000173/1‑93 (MQ=255)
                                                                            aaCATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc              >  1:2989013/1‑93 (MQ=255)
                                                                            aaCATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc              <  2:2273858/93‑1 (MQ=255)
                                                                            aaCATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc              <  1:656641/93‑1 (MQ=255)
                                                                            aaCATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATc              <  1:5198994/93‑1 (MQ=255)
                                                                             aCATGCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCt             >  1:157425/1‑93 (MQ=255)
                                                                               aTGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGc           >  2:3856720/1‑93 (MQ=255)
                                                                                 gCTTGCCGGCGCGACTGTGTGCTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTGTGTTGATTATCTGccc         >  2:6287712/1‑92 (MQ=255)
                                                                                 gCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCa         <  1:6148843/93‑1 (MQ=255)
                                                                                 gCTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCa         >  2:5775359/1‑93 (MQ=255)
                                                                                  cTTGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAg        <  1:4831200/93‑1 (MQ=255)
                                                                                   ttGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgc       <  1:2625315/93‑1 (MQ=255)
                                                                                   ttGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgc       >  2:2720529/1‑93 (MQ=255)
                                                                                    tGCCGGCGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAgcg      >  2:5679299/1‑93 (MQ=255)
                                                                                          |                                                                                          
GGCTCCACAAGACCTGATGCAAATCAACGTCCTTTTGCCATTCATAAACATCATCTGGATAATCGGCATGAGCAATAACATGCTTGCCGGTGCGACTGTGTACTTCAAGAATAAAGGGCAAATGCTGACATTTTTCCAGCGCCGTTTTCGCTTGTTTCTGTTGATTATCTGCCAGCGCAAT  >  NC_000913/1920501‑1920681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: