Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,041,927 A→G Y50C (TAC→TGC)  b3861 → putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,9270AG99.5% 602.9 / ‑5.0 186Y50C (TAC→TGC) b3861putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (101/84);  total (102/84)
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.70e-01

Only 100 of 189 total aligned reads displayed.

ACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCG  >  NC_000913/4041839‑4042016
                                                                                        |                                                                                         
acaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGat                                                                                       <  2:432964/93‑1 (MQ=255)
 caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGata                                                                                      >  1:6145881/1‑93 (MQ=255)
 caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGata                                                                                      >  2:388164/1‑93 (MQ=255)
 caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGata                                                                                      >  2:5494793/1‑93 (MQ=255)
   tGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATg                                                                                    <  1:3523526/93‑1 (MQ=255)
   tGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATg                                                                                    >  1:4228665/1‑93 (MQ=255)
    ggTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGa                                                                                   <  1:975616/93‑1 (MQ=255)
     gTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGAc                                                                                  <  2:4055375/93‑1 (MQ=255)
      tACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACg                                                                                 >  1:2282565/1‑93 (MQ=255)
       acacAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGt                                                                                >  2:4830370/1‑93 (MQ=255)
        cacaGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGtt                                                                               <  1:335888/93‑1 (MQ=255)
          caGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTaa                                                                             <  2:273789/93‑1 (MQ=255)
           aGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAg                                                                            >  1:4718306/1‑93 (MQ=255)
           aGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAg                                                                            >  2:2889461/1‑93 (MQ=255)
            ggTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGt                                                                           <  2:81351/93‑1 (MQ=255)
            ggTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGt                                                                           >  1:3717806/1‑93 (MQ=255)
              tACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTa                                                                         >  2:4955044/1‑93 (MQ=255)
              tACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTa                                                                         >  2:5815299/1‑93 (MQ=255)
                cGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATc                                                                       >  2:5124359/1‑93 (MQ=255)
                cGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATc                                                                       >  2:3012857/1‑93 (MQ=255)
                cGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATc                                                                       <  1:1440220/93‑1 (MQ=255)
                 gCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCa                                                                      >  1:4479728/1‑93 (MQ=255)
                  cAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaa                                                                     >  1:4729514/1‑93 (MQ=255)
                  cAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaa                                                                     >  1:2070730/1‑93 (MQ=255)
                   aaaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaa                                                                    <  1:3589534/93‑1 (MQ=255)
                   aaaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaa                                                                    >  1:559888/1‑93 (MQ=255)
                    aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   >  1:839363/1‑93 (MQ=255)
                    aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   <  2:1772352/93‑1 (MQ=255)
                     aTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAAt                                                                  <  2:1705440/93‑1 (MQ=255)
                     aTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAAt                                                                  >  1:2677279/1‑93 (MQ=255)
                      ttCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAAtt                                                                 >  1:6242263/1‑93 (MQ=255)
                       tCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAAttt                                                                <  2:419544/93‑1 (MQ=255)
                             tGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACaa                                                          >  1:19129/1‑93 (MQ=255)
                              ggTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACaaa                                                         >  1:1077270/1‑93 (MQ=255)
                              ggTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACaaa                                                         <  1:5376321/93‑1 (MQ=255)
                               gtgtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAg                                                        >  1:3908401/1‑93 (MQ=255)
                               gtgtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAg                                                        <  2:69386/93‑1 (MQ=255)
                                 gtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGc                                                      <  2:4550893/93‑1 (MQ=255)
                                  tCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCt                                                     <  2:3482526/93‑1 (MQ=255)
                                    gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                   >  2:3632282/1‑93 (MQ=255)
                                     aaaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttt                                                  <  1:113711/93‑1 (MQ=255)
                                     aaaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttt                                                  <  1:1460834/93‑1 (MQ=255)
                                      aaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttttt                                                 >  1:6036624/1‑93 (MQ=255)
                                       aaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGGCAAAGGCtttttt                                                <  2:1579859/93‑1 (MQ=255)
                                       aaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttttt                                                >  2:3042645/1‑93 (MQ=255)
                                       aaaTTTATTAAAGATTTGAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttttt                                                <  2:5198584/93‑1 (MQ=255)
                                        aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTa                                               >  2:1658487/1‑93 (MQ=255)
                                        aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTa                                               >  2:1712156/1‑93 (MQ=255)
                                        aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGATATCAAAATTTTGACAAAGGCTTTTTTa                                               >  2:4229565/1‑93 (MQ=255)
                                         aTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGTCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTaa                                              >  1:4542884/1‑93 (MQ=255)
                                            tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTc                                           >  1:1964844/1‑93 (MQ=255)
                                            tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTc                                           >  2:221373/1‑93 (MQ=255)
                                            tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGACTGCGATATGCCGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTc                                           >  1:2008644/1‑93 (MQ=255)
                                             attaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGACTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcc                                          <  2:593431/93‑1 (MQ=255)
                                               taAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCg                                        <  2:2955404/93‑1 (MQ=255)
                                                aaaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcccat                                       <  2:1290336/93‑3 (MQ=255)
                                                aaaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGt                                       >  2:2272532/1‑93 (MQ=255)
                                                 aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGtt                                      <  2:1421103/93‑1 (MQ=255)
                                                 aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGtt                                      >  2:4655217/1‑93 (MQ=255)
                                                    aTTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTc                                   >  2:560315/1‑93 (MQ=255)
                                                     ttttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTca                                  >  2:2698877/1‑93 (MQ=255)
                                                      tttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaa                                 >  2:1167860/1‑93 (MQ=255)
                                                       ttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaa                                >  1:5777424/1‑93 (MQ=255)
                                                           cGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgca                            >  2:6253750/1‑93 (MQ=255)
                                                            gATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaa                           >  1:5775706/1‑93 (MQ=255)
                                                             aTGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaa                          >  2:3421499/1‑93 (MQ=255)
                                                             aTGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaa                          >  2:6171269/1‑93 (MQ=255)
                                                              tGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaat                         <  1:2534308/93‑1 (MQ=255)
                                                               gCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaatg                        <  2:2134023/93‑1 (MQ=255)
                                                                cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       >  1:3049978/1‑93 (MQ=255)
                                                                cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       <  1:170304/93‑1 (MQ=255)
                                                                cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       <  2:5694141/93‑1 (MQ=255)
                                                                 tAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGAc                      >  2:553008/1‑93 (MQ=255)
                                                                  aaaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACa                     <  2:110514/93‑1 (MQ=255)
                                                                   aaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAt                    >  1:4360741/1‑93 (MQ=255)
                                                                   aaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAt                    >  1:1798385/1‑93 (MQ=255)
                                                                    aaagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                   <  1:2335954/93‑1 (MQ=255)
                                                                    aaagaaagATGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                   <  2:3138715/93‑1 (MQ=255)
                                                                     aagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTc                  <  2:5064754/93‑1 (MQ=255)
                                                                     aagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTc                  <  2:522707/93‑1 (MQ=255)
                                                                     aagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTc                  <  2:1023799/93‑1 (MQ=255)
                                                                      agaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCg                 >  1:5158019/1‑93 (MQ=255)
                                                                       gaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                >  2:5905506/1‑93 (MQ=255)
                                                                       gaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                <  1:2002757/93‑1 (MQ=255)
                                                                        aaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGAt               >  1:1265462/1‑93 (MQ=255)
                                                                           ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            >  2:4784632/1‑93 (MQ=255)
                                                                           ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            >  2:3907152/1‑93 (MQ=255)
                                                                           ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            <  2:2794190/93‑1 (MQ=255)
                                                                           ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            <  1:5301631/93‑1 (MQ=255)
                                                                           ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            >  1:786393/1‑93 (MQ=255)
                                                                             tGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACgg          <  2:1089536/93‑1 (MQ=255)
                                                                              gAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGt         >  1:3465848/1‑93 (MQ=255)
                                                                               aaCATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTg        >  1:4930607/1‑93 (MQ=255)
                                                                                aCATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGc       <  1:5716678/93‑1 (MQ=255)
                                                                                 cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      >  2:4453573/1‑93 (MQ=255)
                                                                                 cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      >  2:6036705/1‑93 (MQ=255)
                                                                                 cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      >  1:5811801/1‑93 (MQ=255)
                                                                                 cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      >  1:729573/1‑93 (MQ=255)
                                                                                  aTGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAc     >  2:3699003/1‑93 (MQ=255)
                                                                                   tGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAcc    >  1:3182848/1‑93 (MQ=255)
                                                                                        |                                                                                         
ACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCG  >  NC_000913/4041839‑4042016

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