Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 736,618 T→C Y98H (TAT→CAT)  ybfL → putative receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913736,6180TC99.7% 977.7 / ‑5.4 346Y98H (TAT→CAT) ybfLputative receptor protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (163/182);  total (163/183)
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.69e-01

Only 100 of 348 total aligned reads displayed.

TTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCAT  >  NC_000913/736538‑736707
                                                                                |                                                                                         
ttAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTg                                                                               <  2:8185907/93‑1 (MQ=17)
  aaaGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTa                                                                             <  2:3717877/93‑1 (MQ=18)
    aGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTaaa                                                                           >  1:9172357/1‑93 (MQ=18)
    aGGAAAAATCATCACAACTGATGCGATGGGGTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTaaa                                                                           >  2:1262883/1‑93 (MQ=16)
     ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTaaaa                                                                          >  1:4277454/1‑93 (MQ=18)
     ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTaaaa                                                                          <  1:841013/93‑1 (MQ=21)
     ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTaaaa                                                                          >  2:5590752/1‑93 (MQ=18)
        aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGa                                                                       <  1:5219803/93‑1 (MQ=18)
        aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGa                                                                       >  1:412887/1‑93 (MQ=18)
        aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGa                                                                       <  2:7172402/93‑1 (MQ=21)
        aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGAACATTTATTCGCTGTAAAAGGa                                                                       <  1:3928808/93‑1 (MQ=255)
         aaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaa                                                                      >  2:1767662/1‑93 (MQ=18)
         aaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaa                                                                      >  1:10010427/1‑93 (MQ=18)
          aaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaaa                                                                     <  1:9474490/93‑1 (MQ=18)
          aaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaaa                                                                     >  1:318725/1‑93 (MQ=18)
          aaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaaa                                                                     <  1:1178464/93‑1 (MQ=18)
          aaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATAGTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaaa                                                                     <  1:8094500/93‑1 (MQ=16)
          aaTCATCACAACTGATGCGATGGATTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGaaa                                                                     >  1:6697503/1‑93 (MQ=16)
           atcatcACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAAc                                                                    <  1:5463839/93‑1 (MQ=21)
            tcatcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAAcc                                                                   >  2:2261109/1‑93 (MQ=18)
             catcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCa                                                                  <  2:6126819/93‑1 (MQ=21)
              atcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAg                                                                 <  2:4084075/93‑1 (MQ=21)
              atcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAg                                                                 <  2:9638874/93‑1 (MQ=21)
                  caACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAggggg                                                             >  2:6696928/1‑92 (MQ=18)
                  caACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggc                                                             >  2:8757350/1‑93 (MQ=18)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTACAAGGAAACCAGGggcg                                                            >  2:3835612/1‑93 (MQ=16)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcg                                                            <  1:1312823/93‑1 (MQ=18)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcg                                                            >  2:2136246/1‑93 (MQ=18)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcg                                                            >  1:936442/1‑93 (MQ=18)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcg                                                            >  1:8044131/1‑93 (MQ=18)
                   aaCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcg                                                            >  1:6956357/1‑93 (MQ=18)
                    aCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGggcgg                                                           >  1:5867425/1‑93 (MQ=18)
                      tGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCt                                                         <  1:5825976/93‑1 (MQ=14)
                        aTGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTta                                                       >  1:2713405/1‑93 (MQ=17)
                         tGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaa                                                      >  1:9494111/1‑93 (MQ=18)
                         tGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaa                                                      <  2:7455472/93‑1 (MQ=18)
                           cGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaata                                                    <  2:7341587/93‑1 (MQ=18)
                            gATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaataa                                                   <  1:975230/93‑1 (MQ=21)
                            gATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaataa                                                   <  2:6378546/93‑1 (MQ=18)
                            gATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaataa                                                   >  2:2984896/1‑93 (MQ=18)
                            gATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaataa                                                   <  1:3568587/93‑1 (MQ=21)
                            gATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTtaataa                                                   <  1:1774285/93‑1 (MQ=21)
                             aTGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATaaa                                                  >  2:8038981/1‑93 (MQ=18)
                             aTGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATaaa                                                  <  1:2086068/93‑1 (MQ=21)
                              tGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAg                                                 <  2:7828963/93‑1 (MQ=18)
                              tGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAg                                                 <  2:1170366/93‑1 (MQ=21)
                              tGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAg                                                 >  1:2262029/1‑93 (MQ=18)
                               gggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGc                                                >  1:2561908/1‑93 (MQ=18)
                               gggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGc                                                >  1:7587608/1‑93 (MQ=18)
                               gggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGc                                                <  2:9496273/93‑1 (MQ=21)
                                ggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGcc                                               <  1:9042597/93‑1 (MQ=18)
                                 gTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCt                                              <  1:1600257/93‑1 (MQ=21)
                                   tGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCttt                                            <  2:4829784/93‑1 (MQ=25)
                                   tGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCttt                                            <  2:7611237/93‑1 (MQ=25)
                                   tGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCttt                                            <  1:599769/93‑1 (MQ=25)
                                    gCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTg                                           >  1:8202881/1‑93 (MQ=25)
                                    gCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGACTTTg                                           >  2:6921338/1‑93 (MQ=255)
                                     ccAGAAAGATATTGCAGAGAAGATACAAAGACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGa                                          <  1:9221718/93‑1 (MQ=255)
                                      cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAg                                         >  1:4356290/1‑93 (MQ=255)
                                      cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAg                                         <  1:5265047/93‑1 (MQ=255)
                                      cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAg                                         <  1:1557123/93‑1 (MQ=255)
                                       aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAgg                                        >  2:9179675/1‑93 (MQ=37)
                                       aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAgg                                        >  2:9026321/1‑93 (MQ=255)
                                       aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAgg                                        >  2:5125579/1‑93 (MQ=37)
                                       aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAgg                                        <  2:305307/93‑1 (MQ=255)
                                        gAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGa                                       >  1:4643140/1‑93 (MQ=255)
                                          aaGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGaaa                                     >  2:5706757/1‑93 (MQ=255)
                                            gATATTGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGaaaaa                                   >  1:6897569/1‑93 (MQ=255)
                                                ttGCAGAGAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTc                               <  1:2286240/93‑1 (MQ=255)
                                                    agagAAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCt                           >  1:7268003/1‑93 (MQ=255)
                                                     gagaAGATACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCta                          >  1:4865653/1‑92 (MQ=255)
                                                       gaagaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGaa                        >  1:5635323/1‑93 (MQ=255)
                                                        aagaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                       >  2:8916976/1‑93 (MQ=255)
                                                        aagaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                       <  2:7252310/93‑1 (MQ=37)
                                                         agaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAg                      <  2:2456699/93‑1 (MQ=255)
                                                         agaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAg                      >  2:1942804/1‑93 (MQ=255)
                                                         agaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAg                      <  2:2860038/93‑1 (MQ=255)
                                                          gaTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGa                     <  2:4534016/93‑1 (MQ=255)
                                                           aTACAAATACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGaa                    <  2:7689769/93‑1 (MQ=255)
                                                           aTACAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGaa                    >  2:592363/1‑93 (MQ=255)
                                                              cAAAAACAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTa                 >  1:3029687/1‑93 (MQ=255)
                                                                aaaaCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTaaa               >  1:235701/1‑93 (MQ=255)
                                                                aaaaCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTaaa               <  1:5564797/93‑1 (MQ=255)
                                                                 aaaCAGGGAGGTGATGATTTAGTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaat              <  1:7773161/93‑1 (MQ=255)
                                                                  aaCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaata             <  2:265394/93‑1 (MQ=37)
                                                                   aCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataa            <  2:847598/93‑1 (MQ=37)
                                                                   aCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataa            >  1:9107940/1‑93 (MQ=255)
                                                                   aCAGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataa            <  2:3472834/93‑1 (MQ=255)
                                                                     aGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATc          >  1:834170/1‑93 (MQ=255)
                                                                     aGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATc          <  1:7539997/93‑1 (MQ=255)
                                                                     aGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATc          >  1:3486237/1‑93 (MQ=255)
                                                                     aGGGAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATc          >  1:3750503/1‑93 (MQ=255)
                                                                       ggAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa        >  1:7039717/1‑93 (MQ=255)
                                                                       ggAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa        <  2:5813291/93‑1 (MQ=37)
                                                                       ggAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCa        <  1:6578340/93‑1 (MQ=255)
                                                                        gAGGTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCag       <  2:3268994/93‑1 (MQ=255)
                                                                           gTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGc    >  2:65615/1‑93 (MQ=255)
                                                                           gTGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGc    <  1:9968250/93‑1 (MQ=255)
                                                                            tGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCa   <  2:153279/93‑1 (MQ=255)
                                                                            tGATCATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCa   >  1:2425590/1‑93 (MQ=255)
                                                                                |                                                                                         
TTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCTGTAAAAGGAAACCAGGGGCGGCTTAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCAT  >  NC_000913/736538‑736707

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