Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,005,601 +C coding (237/282 nt) b4568 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,005,6001.C98.6% 660.7 / 1.4 210coding (236/282 nt)b4568hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (122/85);  total (124/86)
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.97e-01

Only 100 of 228 total aligned reads displayed.

GTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCT  >  NC_000913/4005509‑4005686
                                                                                              |                                                                                      
 tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcc                                                                                      >  1:2758048/1‑92 (MQ=255)
   aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcctg                                                                                    >  2:40177/1‑90 (MQ=255)
    tGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGcctgt                                                                                   <  1:3692418/93‑5 (MQ=255)
       tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtc                                                                                 >  2:2731667/1‑93 (MQ=255)
       tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtc                                                                                 <  2:4591410/93‑1 (MQ=255)
         tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                                               <  1:4416891/93‑1 (MQ=255)
         tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                                               <  2:3058129/93‑1 (MQ=255)
          gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGTCTGGCTGCATACGCTGATTATGGTGGTCGGCGGGCTGGATctct                                                                              >  2:4487259/1‑93 (MQ=255)
            gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGACCTGTATCTctgc                                                                            >  2:359421/1‑93 (MQ=255)
            gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgc                                                                            <  2:1942111/93‑1 (MQ=255)
              tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                          >  1:1707122/1‑93 (MQ=255)
              tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                          <  1:1002317/93‑1 (MQ=255)
               gggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                                                                         <  2:4632469/93‑1 (MQ=255)
                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                        >  1:1651102/1‑93 (MQ=255)
                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                        >  1:321128/1‑93 (MQ=255)
                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                        <  1:4529405/93‑1 (MQ=255)
                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGAATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                        <  1:4314411/93‑1 (MQ=255)
                 gATTGCGCTGCTTGTCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGAt                                                                       <  1:1586664/93‑1 (MQ=255)
                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGTTg                                                                      >  1:491787/1‑93 (MQ=255)
                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                      <  1:2696223/93‑1 (MQ=255)
                   ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                                     <  2:4113172/93‑1 (MQ=255)
                     gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                                   <  1:642497/93‑1 (MQ=255)
                       gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                                 >  2:3156086/1‑93 (MQ=255)
                       gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                                 >  2:3734096/1‑93 (MQ=255)
                       gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGGTa                                                                 >  2:3011818/1‑93 (MQ=255)
                        ctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc                                                                >  1:3036038/1‑93 (MQ=255)
                        ctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc                                                                <  1:133345/93‑1 (MQ=255)
                         tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcc                                                               >  2:369410/1‑93 (MQ=255)
                           ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAg                                                             >  1:1636543/1‑93 (MQ=255)
                            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                            >  1:2100722/1‑93 (MQ=255)
                            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                            <  2:110511/93‑1 (MQ=255)
                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGATGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAgctgct                                                        >  1:4351681/1‑93 (MQ=255)
                                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                       >  2:1669105/1‑93 (MQ=255)
                                 cTGCATTTGATTATCGAAAAAAAGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCATATGCTa                                                       >  2:698806/1‑93 (MQ=255)
                                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                      >  1:1118451/1‑93 (MQ=255)
                                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                      >  1:2649228/1‑93 (MQ=255)
                                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgg                                                    >  2:145885/1‑92 (MQ=255)
                                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  1:3394862/1‑93 (MQ=255)
                                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   <  1:2568917/93‑1 (MQ=255)
                                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  2:4287061/1‑93 (MQ=255)
                                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                  <  2:2931835/93‑1 (MQ=255)
                                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                  >  2:854012/1‑93 (MQ=255)
                                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                  >  1:3995372/1‑93 (MQ=255)
                                       ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                                 <  2:2078848/93‑1 (MQ=255)
                                       ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                                 <  1:2591597/93‑1 (MQ=255)
                                         gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                               <  1:2064118/93‑1 (MQ=255)
                                          aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                                              <  2:3260454/93‑1 (MQ=255)
                                          aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                                              >  1:4227350/1‑93 (MQ=255)
                                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                             >  1:2048387/1‑93 (MQ=255)
                                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                             >  1:3152473/1‑93 (MQ=255)
                                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                             >  1:638223/1‑93 (MQ=255)
                                            tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                            >  1:2311541/1‑93 (MQ=255)
                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                          >  1:3533701/1‑93 (MQ=255)
                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                          >  1:701957/1‑93 (MQ=255)
                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                          >  1:4069757/1‑93 (MQ=255)
                                               cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                                         >  2:3421150/1‑93 (MQ=255)
                                               cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                                         <  2:1262684/93‑1 (MQ=255)
                                                  aaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaa                                      <  1:3137195/93‑1 (MQ=255)
                                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                     >  2:1537641/1‑93 (MQ=255)
                                                    aaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAag                                    >  1:776222/1‑93 (MQ=255)
                                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAACaga                                   >  1:1487450/1‑93 (MQ=255)
                                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                   <  2:532021/93‑1 (MQ=255)
                                                      aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                  <  1:2852546/93‑1 (MQ=255)
                                                        ggCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                <  2:1493065/93‑1 (MQ=255)
                                                          ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg                              >  1:1266064/1‑93 (MQ=255)
                                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                             >  1:2971107/1‑93 (MQ=255)
                                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                             >  1:1499725/1‑93 (MQ=255)
                                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                             >  1:3492343/1‑93 (MQ=255)
                                                             ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt                           <  1:74532/93‑1 (MQ=255)
                                                               cTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt                         <  2:2501206/93‑1 (MQ=255)
                                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                       <  1:2400468/93‑1 (MQ=255)
                                                                   aTACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                     >  2:1814384/1‑93 (MQ=255)
                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                    <  1:1402481/93‑1 (MQ=255)
                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                    >  1:2369459/1‑93 (MQ=255)
                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTACTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                    >  1:2228199/1‑93 (MQ=255)
                                                                     aCGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt                   >  2:2112690/1‑93 (MQ=255)
                                                                      cGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTg                  >  1:975195/1‑93 (MQ=255)
                                                                       gCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGGGCACTGAAAAAAGAGGCGGTTTCTGCACCTgc                 >  2:3132276/1‑93 (MQ=255)
                                                                        cTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcg                >  2:3513575/1‑93 (MQ=255)
                                                                         tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc               >  2:2262612/1‑93 (MQ=255)
                                                                          gATTATGGTGGGCGGTGGCCTGTATTTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc              <  2:4150414/93‑1 (MQ=255)
                                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc              >  2:773909/1‑93 (MQ=255)
                                                                            ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            >  1:2876685/1‑93 (MQ=255)
                                                                            ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            >  1:3771032/1‑93 (MQ=255)
                                                                            ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            >  1:1531829/1‑93 (MQ=255)
                                                                            ttATGGTGGGCGATGTCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            <  2:321651/93‑1 (MQ=255)
                                                                             tATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca           <  1:1995210/93‑1 (MQ=255)
                                                                             tATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca           >  2:589254/1‑93 (MQ=255)
                                                                              aTGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAg          <  2:1718252/93‑1 (MQ=255)
                                                                              aTGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAg          <  1:3577316/93‑1 (MQ=255)
                                                                               tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  1:4278062/1‑93 (MQ=255)
                                                                               tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGGTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  1:4479676/1‑93 (MQ=255)
                                                                               tggtggGCGGGGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  1:3631573/1‑93 (MQ=255)
                                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAggc        >  2:2299866/1‑92 (MQ=255)
                                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt        <  1:3588623/93‑1 (MQ=255)
                                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt        <  2:4423123/93‑1 (MQ=255)
                                                                                ggtggGCGGTGGCCTGTATCTCGGCTGGAGGGGTTCCCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCCCCTGCGCCACAGGt        <  2:619615/93‑1 (MQ=255)
                                                                                 gtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTc       <  1:3813073/93‑1 (MQ=255)
                                                                                     gCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc   <  1:2454980/93‑1 (MQ=255)
                                                                                      cGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCt  >  1:219485/1‑93 (MQ=255)
                                                                                              |                                                                                      
GTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCT  >  NC_000913/4005509‑4005686

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