Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,017,437 T→C D29G (GAC→GGC)  b0955 ← putative ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,017,4370TC97.7% 543.4 / ‑0.3 175D29G (GAC→GGC) b0955putative ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base C (103/68);  total (107/68)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 190 total aligned reads displayed.

TAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTT  >  NC_000913/1017348‑1017526
                                                                                         |                                                                                         
tAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGccaa                                                                                        <  2:1678207/93‑1 (MQ=255)
  aaaGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGccaatc                                                                                      >  2:91177/1‑93 (MQ=255)
  aaaGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGccaatc                                                                                      <  2:1544815/93‑1 (MQ=255)
     gaggagGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCaaa                                                                                   >  1:3480691/1‑93 (MQ=255)
     gaggagGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCaaa                                                                                   <  1:2044744/93‑1 (MQ=255)
      aggaggATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAt                                                                                  >  2:1062459/1‑93 (MQ=255)
      aggaggATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAt                                                                                  <  1:3786924/93‑1 (MQ=255)
      aggaggATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAt                                                                                  <  1:3786923/93‑1 (MQ=255)
      aggaggATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTAGCACTGAATAAAGGATCGTTTTCGCCAATCAAAt                                                                                  >  1:4583562/1‑93 (MQ=255)
        gaggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAtgt                                                                                >  1:1014622/1‑93 (MQ=255)
        gaggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAtgt                                                                                <  2:1562887/93‑1 (MQ=255)
        gaggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAtgt                                                                                <  2:3169469/93‑1 (MQ=255)
        gaggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAtgt                                                                                <  1:2702794/93‑1 (MQ=255)
        gaggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAAtgt                                                                                >  1:1495882/1‑93 (MQ=255)
          ggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTgg                                                                              >  2:246665/1‑93 (MQ=255)
          ggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTgg                                                                              >  1:3954468/1‑93 (MQ=255)
          ggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTgg                                                                              >  1:350480/1‑93 (MQ=255)
          ggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTgg                                                                              <  2:3651288/93‑1 (MQ=255)
           gaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGc                                                                             <  1:3077577/93‑1 (MQ=255)
           gaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGc                                                                             >  1:4634561/1‑93 (MQ=255)
           gaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGc                                                                             >  1:1376139/1‑93 (MQ=255)
            aTGCTGTCGAATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCt                                                                            <  1:4456093/93‑1 (MQ=255)
             tGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTg                                                                           >  2:4105557/1‑93 (MQ=255)
             tGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTg                                                                           <  2:2503344/93‑1 (MQ=255)
              gCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGa                                                                          >  1:4749861/1‑93 (MQ=255)
                tCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGc                                                                        >  2:657546/1‑93 (MQ=255)
                 cgcgTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCa                                                                       >  1:1705195/1‑93 (MQ=255)
                 cgcgTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCa                                                                       >  1:2489045/1‑93 (MQ=255)
                   cgTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCaaa                                                                     >  2:3002294/1‑93 (MQ=255)
                     tATTCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAata                                                                   >  1:2096138/1‑93 (MQ=255)
                     tATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAata                                                                   <  1:2130634/93‑1 (MQ=255)
                     tATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAata                                                                   >  2:4138476/1‑93 (MQ=255)
                        gcagcaGCTGCAAAAGCGCAAATTGCAGAAGCAATTGAGTATCACTGAATAAAGAATCGATTTCGCAAATCAAATATGAAAAAGCAAATAttt                                                                <  2:3602206/93‑1 (MQ=255)
                        gcagcaACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATAttt                                                                <  2:2723260/93‑1 (MQ=255)
                           gcaACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCt                                                             <  2:1080877/93‑1 (MQ=255)
                              aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGAt                                                          >  2:3777606/1‑93 (MQ=255)
                              aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGAt                                                          >  1:4689421/1‑93 (MQ=255)
                              aCTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGAt                                                          <  1:4425456/93‑1 (MQ=255)
                               cTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATa                                                         <  2:3283047/93‑1 (MQ=255)
                                 gCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGc                                                       <  1:3570632/93‑1 (MQ=255)
                                  cTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCt                                                      <  2:3384805/93‑1 (MQ=255)
                                  cTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCt                                                      >  2:3974874/1‑93 (MQ=255)
                                   tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTa                                                     >  1:4532817/1‑93 (MQ=255)
                                   tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTa                                                     >  2:4731368/1‑93 (MQ=255)
                                    ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTAt                                                    <  2:3038947/93‑1 (MQ=255)
                                    ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTAt                                                    >  1:1923078/1‑93 (MQ=255)
                                     cAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATc                                                   <  2:571384/93‑1 (MQ=255)
                                       gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCgg                                                 >  2:4631534/1‑93 (MQ=255)
                                          cAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTAt                                              >  2:2110564/1‑93 (MQ=255)
                                           aaaTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                             <  2:1925342/93‑1 (MQ=255)
                                             aTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAg                                           >  2:4557466/1‑93 (MQ=255)
                                              ttGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAgg                                          <  2:4321794/93‑1 (MQ=255)
                                                gCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGaa                                        <  2:2774819/93‑1 (MQ=255)
                                                gCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGaa                                        <  1:3028257/93‑1 (MQ=255)
                                                  aGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAACTATCGGTATCAGGAAcc                                      >  1:995216/1‑93 (MQ=255)
                                                   gCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCa                                     >  1:1102714/1‑93 (MQ=255)
                                                    ccGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAg                                    <  1:3313322/93‑1 (MQ=255)
                                                       cgGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTc                                 >  2:3837500/1‑93 (MQ=255)
                                                       cgGTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTc                                 <  1:2594858/93‑1 (MQ=255)
                                                        ggTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCACATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCa                                >  2:3322582/1‑93 (MQ=255)
                                                        ggTTGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTACTGATAGCTATCGGTATCAGGAACAAGGTCa                                >  1:1285177/1‑93 (MQ=255)
                                                           tGAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGc                             <  1:3162648/93‑1 (MQ=255)
                                                            gAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGcc                            >  2:4044997/1‑93 (MQ=255)
                                                            gAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGcc                            >  2:504620/1‑93 (MQ=255)
                                                            gAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGcc                            <  2:2298677/93‑1 (MQ=255)
                                                            gAGTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGcc                            >  1:80371/1‑93 (MQ=255)
                                                              gTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATTTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCAt                          >  2:4064749/1‑93 (MQ=255)
                                                              gTATCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCAt                          <  1:3820999/93‑1 (MQ=255)
                                                                aTCACTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGc                        >  2:1409413/1‑93 (MQ=255)
                                                                   aCTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAg                     >  2:4746498/1‑93 (MQ=255)
                                                                    cTGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGt                    >  1:3406216/1‑93 (MQ=255)
                                                                     tGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGtt                   >  2:4347974/1‑93 (MQ=255)
                                                                     tGAATAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGtt                   >  1:2041734/1‑93 (MQ=255)
                                                                       aaTAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTc                 >  1:2784008/1‑93 (MQ=255)
                                                                         tAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGt               >  1:228654/1‑93 (MQ=255)
                                                                         tAAAGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGt               >  1:2796099/1‑93 (MQ=255)
                                                                          aaaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa              >  2:346749/1‑93 (MQ=255)
                                                                          aaaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa              >  2:264020/1‑93 (MQ=255)
                                                                          aaaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTa              >  2:2061432/1‑93 (MQ=255)
                                                                           aaGGATCGTTTTCGCTAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTaa             <  2:2935954/93‑1 (MQ=255)
                                                                           aaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTaa             <  1:4174929/93‑1 (MQ=255)
                                                                           aaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTaa             >  1:2281618/1‑93 (MQ=255)
                                                                           aaGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTaa             <  1:1739262/93‑1 (MQ=255)
                                                                            aGGATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAAt            <  1:609825/93‑1 (MQ=255)
                                                                             ggATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg           <  2:3552849/93‑1 (MQ=255)
                                                                             ggATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg           <  1:4032668/93‑1 (MQ=255)
                                                                             ggATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATg           >  2:1379736/1‑93 (MQ=255)
                                                                              gATCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATgg          >  1:3637026/1‑93 (MQ=255)
                                                                              gATCGTATTCGCCAAAAAAATATGGCAGAACAAATATTTCCTGATAGCTAACGGTATAAGAAACCAGGTCACGACATGAAAGTTACGTAAtag          <  2:4349644/93‑3 (MQ=255)
                                                                                tCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc        >  2:3316161/1‑93 (MQ=255)
                                                                                tCGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc        <  2:1347105/93‑1 (MQ=255)
                                                                                 cGTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCa       >  2:1198574/1‑93 (MQ=255)
                                                                                  gTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaa      >  2:1611079/1‑93 (MQ=255)
                                                                                  gTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaa      >  1:1835354/1‑93 (MQ=255)
                                                                                  gTTTTCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaa      >  2:2002479/1‑93 (MQ=255)
                                                                                   ttttCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaaa     <  2:3801784/93‑1 (MQ=255)
                                                                                   ttttCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaaa     <  1:4342071/93‑1 (MQ=255)
                                                                                   ttttCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaaa     >  2:2065196/1‑93 (MQ=255)
                                                                                   ttttCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaaa     >  1:1807703/1‑93 (MQ=255)
                                                                                   ttttCGCCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCaaa     >  2:1790543/1‑93 (MQ=255)
                                                                                         |                                                                                         
TAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTT  >  NC_000913/1017348‑1017526

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