Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,782,327 T→C V41V (GTA→GTG b2649 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,782,3270TC86.2% 422.3 / 56.7 167V41V (GTA→GTGb2649hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/9);  new base C (64/80);  total (78/89)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-02

Only 100 of 179 total aligned reads displayed.

TATCTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATC  >  NC_000913/2782236‑2782416
                                                                                           |                                                                                         
tATCTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAataata                                                                                          >  1:6523014/1‑93 (MQ=255)
 aTCTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAacac                                                                                         <  1:3176218/93‑1 (MQ=255)
 aTCTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAacac                                                                                         <  2:5851311/93‑1 (MQ=255)
   cTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACAccc                                                                                       <  2:3340232/93‑1 (MQ=255)
    ttGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATAcccc                                                                                      >  2:6376774/1‑93 (MQ=255)
    ttGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACAcccc                                                                                      >  2:5436552/1‑93 (MQ=255)
    ttGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACAcccc                                                                                      <  1:5923640/93‑1 (MQ=255)
    ttGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACAcccc                                                                                      <  2:4054308/93‑1 (MQ=255)
      gAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCt                                                                                    <  2:5601464/93‑1 (MQ=255)
       aGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCtt                                                                                   <  2:4119211/93‑1 (MQ=255)
       aGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCtt                                                                                   <  2:1550001/93‑1 (MQ=255)
       aGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCtt                                                                                   >  1:3102886/1‑93 (MQ=255)
        ggCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTg                                                                                  >  1:2814891/1‑93 (MQ=255)
        ggCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTg                                                                                  >  1:5407213/1‑93 (MQ=255)
        ggCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTg                                                                                  >  2:1411856/1‑93 (MQ=255)
            gAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACt                                                                              >  1:5050005/1‑93 (MQ=255)
            gAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACt                                                                              <  1:1245270/93‑1 (MQ=255)
            gAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACt                                                                              <  1:300813/93‑1 (MQ=255)
              gCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCa                                                                            <  2:3749236/93‑1 (MQ=255)
               ctctTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCAt                                                                           <  1:2447177/93‑1 (MQ=255)
               ctctTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCAt                                                                           <  1:4356082/93‑1 (MQ=255)
                 ctTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTc                                                                         <  1:74863/93‑1 (MQ=255)
                  ttGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTca                                                                        >  1:1423246/1‑93 (MQ=255)
                  ttGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTca                                                                        >  2:6430167/1‑93 (MQ=255)
                   tGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTcac                                                                       <  1:5971831/93‑1 (MQ=255)
                     gTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACaa                                                                     >  1:5081366/1‑93 (MQ=255)
                       aaaaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAAc                                                                   <  2:4160940/93‑1 (MQ=255)
                       aaaaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAAc                                                                   <  1:4304047/93‑1 (MQ=255)
                       aaaaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAAc                                                                   <  1:2904486/93‑1 (MQ=255)
                         aaaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTTATTCACAAACGc                                                                 >  2:3885249/1‑93 (MQ=255)
                          aaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACACACGCa                                                                >  1:2175422/1‑93 (MQ=255)
                          aaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCa                                                                <  2:5085670/93‑1 (MQ=255)
                          aaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCa                                                                >  1:1756658/1‑93 (MQ=255)
                          aaCTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCa                                                                <  2:5763170/93‑1 (MQ=255)
                            cTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATc                                                              >  1:3782331/1‑93 (MQ=255)
                            cTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATc                                                              >  1:4691500/1‑93 (MQ=255)
                            cTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATc                                                              >  1:2644579/1‑93 (MQ=255)
                             ttATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCt                                                             <  1:3480857/93‑1 (MQ=255)
                             ctataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCt                                                             <  2:1150711/92‑1 (MQ=255)
                               ataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCttt                                                           <  2:4928763/93‑1 (MQ=255)
                               ataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCttt                                                           <  1:1120205/93‑1 (MQ=255)
                               ataGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCttt                                                           >  1:1731676/1‑93 (MQ=255)
                                taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCtttt                                                          >  1:6235313/1‑93 (MQ=255)
                                taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACTCCCTTGAACTCATTCACAAACGCATCtttt                                                          >  2:1241405/1‑93 (MQ=255)
                                taGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCtttt                                                          <  2:5978766/93‑1 (MQ=255)
                                  gCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCtttttt                                                        <  2:2663570/93‑1 (MQ=255)
                                  gCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCtttttt                                                        >  1:5538750/1‑93 (MQ=255)
                                  gCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCtttttt                                                        <  1:701408/93‑1 (MQ=255)
                                     acacTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTc                                                     <  2:2851561/93‑1 (MQ=255)
                                     acacTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTc                                                     >  1:4202729/1‑93 (MQ=255)
                                      cacTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCt                                                    <  2:1073964/93‑1 (MQ=255)
                                      cacTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCt                                                    <  1:1378216/93‑1 (MQ=255)
                                       acTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCtt                                                   <  1:3395708/93‑1 (MQ=255)
                                       acTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCtt                                                   <  1:1471495/93‑1 (MQ=255)
                                        cTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTg                                                  >  2:6575334/1‑93 (MQ=255)
                                         tCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGa                                                 >  1:4247414/1‑93 (MQ=255)
                                         tCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGa                                                 >  2:3645226/1‑93 (MQ=255)
                                          cATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGaa                                                <  1:989078/93‑1 (MQ=255)
                                            tAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAata                                              <  1:3638928/93‑1 (MQ=255)
                                             aaaaaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAatat                                             >  1:5034051/1‑93 (MQ=255)
                                               aaaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCt                                           >  1:5371304/1‑93 (MQ=255)
                                                aaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTa                                          >  1:4606255/1‑93 (MQ=255)
                                                aaTTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTa                                          <  2:5335038/93‑1 (MQ=255)
                                                  ttAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATg                                        >  1:3963563/1‑93 (MQ=255)
                                                  ttAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATg                                        <  2:394172/93‑1 (MQ=255)
                                                  ttAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATg                                        <  1:3870928/93‑1 (MQ=255)
                                                     aaCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATgggg                                     >  1:6273854/1‑93 (MQ=255)
                                                      aCAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGgt                                    <  1:5636273/93‑1 (MQ=255)
                                                        aGAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGgtgt                                  >  2:745513/1‑93 (MQ=255)
                                                         gAAGCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGtt                                 >  2:2131468/1‑93 (MQ=255)
                                                            gCTCTGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCt                              <  1:5581076/93‑1 (MQ=255)
                                                             ctctGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCtt                             <  2:5707040/93‑1 (MQ=255)
                                                                tGAAGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGa                          <  2:3988712/93‑1 (MQ=255)
                                                                  aaGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGAct                        >  1:4760612/1‑93 (MQ=255)
                                                                  aaGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGAct                        >  1:4956046/1‑93 (MQ=255)
                                                                   aGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGActc                       >  1:2044333/1‑93 (MQ=255)
                                                                   aGACCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGActc                       >  2:5986217/1‑93 (MQ=255)
                                                                    gACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGActct                      <  2:4331772/93‑1 (MQ=255)
                                                                    gACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGActct                      <  1:360050/93‑1 (MQ=255)
                                                                     aCCTTGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTg                     <  1:2717385/93‑1 (MQ=255)
                                                                        ttGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACa                  <  2:5764831/93‑1 (MQ=255)
                                                                        ttGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACa                  <  1:231578/93‑1 (MQ=255)
                                                                         tGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACaa                 >  2:4809228/1‑93 (MQ=255)
                                                                         tGTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACaa                 >  2:5217125/1‑93 (MQ=255)
                                                                          ttATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACaaa                <  2:4285178/92‑1 (MQ=255)
                                                                          gTATTTTTATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACaaa                >  1:5440882/1‑93 (MQ=255)
                                                                            aTTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGa              <  1:1068711/93‑1 (MQ=255)
                                                                             tttttATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGAt             <  2:45004/93‑1 (MQ=255)
                                                                             tttttATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGAt             <  2:2357143/93‑1 (MQ=255)
                                                                             tttttATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGAt             <  1:3998972/93‑1 (MQ=255)
                                                                                 tATCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGaaa         <  2:941025/93‑1 (MQ=255)
                                                                                   tCAAATAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtc       >  1:4432986/1‑93 (MQ=255)
                                                                                     aaaTAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtctc     >  1:4926133/1‑93 (MQ=255)
                                                                                     aaaTAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtctc     >  2:4024341/1‑93 (MQ=255)
                                                                                     aaaTAACACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAAtctc     <  1:1626855/93‑1 (MQ=255)
                                                                                      aataacACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCa    <  2:2581440/93‑1 (MQ=255)
                                                                                       ataacaCCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCAt   >  2:4582660/1‑93 (MQ=255)
                                                                                       ataacaCCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCAt   >  2:2965173/1‑93 (MQ=255)
                                                                                        taacaCCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATc  >  1:5630134/1‑93 (MQ=255)
                                                                                        taacaCCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATc  >  1:5492459/1‑93 (MQ=255)
                                                                                           |                                                                                         
TATCTTGAGGCTGAGCTCTTGGTAAAAACTTATAGCAACACTCATAAAAATTAAACAGAAGCTCTGAAGACCTTGTATTTTTATCAAATAATACCCCCTTGAACTCATTCACAAACGCATCTTTTTTTTCTTGAATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATC  >  NC_000913/2782236‑2782416

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