Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,992 A→G intergenic (‑145/+13) ylaB ← / ← ylaC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,9920AG99.5% 616.1 / ‑5.2 192intergenic (‑145/+13)ylaB/ylaChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (79/112);  total (79/113)
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.72e-01

Only 100 of 201 total aligned reads displayed.

TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081
                                                                                          |                                                                                         
tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGtgtg                                                                                         >  1:905846/1‑93 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        <  1:1428063/93‑1 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        >  2:4228541/1‑93 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        <  2:3584321/93‑1 (MQ=255)
  gagaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTa                                                                                       <  1:2918111/93‑1 (MQ=255)
   agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTaa                                                                                      <  2:823330/93‑1 (MQ=255)
   agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTaa                                                                                      >  1:2871205/1‑93 (MQ=255)
   agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTaa                                                                                      <  2:1255852/93‑1 (MQ=255)
   agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTaa                                                                                      <  1:1329520/93‑1 (MQ=255)
        aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTa                                                                                 >  1:3510977/1‑93 (MQ=255)
         ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTaa                                                                                <  1:3473359/93‑1 (MQ=255)
          gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAAc                                                                               <  1:2234259/93‑1 (MQ=255)
            ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACtt                                                                             >  2:4377111/1‑93 (MQ=255)
              tAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTaa                                                                           <  2:390557/93‑1 (MQ=255)
                gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGa                                                                         <  2:4457057/93‑1 (MQ=255)
                  gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATg                                                                       >  1:3033425/1‑93 (MQ=255)
                      gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGa                                                                   <  2:3802738/93‑1 (MQ=255)
                      gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGa                                                                   <  2:3242705/93‑1 (MQ=255)
                       tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAt                                                                  <  1:671295/93‑1 (MQ=255)
                       tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAt                                                                  <  1:1154909/93‑1 (MQ=255)
                        aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAtt                                                                 >  1:556571/1‑93 (MQ=255)
                        aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAtt                                                                 >  2:1513814/1‑93 (MQ=255)
                         tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTc                                                                <  2:4557942/93‑1 (MQ=255)
                         tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTc                                                                <  2:2350374/93‑1 (MQ=255)
                         tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGGTGATTc                                                                >  2:1829549/1‑93 (MQ=255)
                          ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCg                                                               <  2:1733512/93‑1 (MQ=255)
                          ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCg                                                               <  2:38124/93‑1 (MQ=255)
                           tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCgg                                                              <  1:3537654/93‑1 (MQ=255)
                           tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCgg                                                              >  1:4059775/1‑93 (MQ=255)
                            cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgc                                                             >  2:4638849/1‑93 (MQ=255)
                              tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                           <  2:2901457/93‑1 (MQ=255)
                              tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                           <  2:3080456/93‑1 (MQ=255)
                              tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                           <  1:1521447/93‑1 (MQ=255)
                                aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcga                                                         <  1:4179540/93‑1 (MQ=255)
                                   tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                      >  1:4171420/1‑93 (MQ=255)
                                    taTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     <  2:1615785/92‑1 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     >  1:3649963/1‑93 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     >  1:2543627/1‑93 (MQ=255)
                                      tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   >  2:565808/1‑93 (MQ=255)
                                      tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   <  1:636219/93‑1 (MQ=255)
                                       gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                  <  2:2535527/93‑1 (MQ=255)
                                       gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                  >  1:2434915/1‑93 (MQ=255)
                                        ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                 <  2:1425020/93‑1 (MQ=255)
                                           ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                              <  2:4473788/93‑1 (MQ=255)
                                           ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                              <  1:1081279/93‑1 (MQ=255)
                                            gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                             >  2:4314541/1‑93 (MQ=255)
                                            gAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGGGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                             >  1:2462354/1‑93 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACTCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           <  2:2406815/93‑1 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           >  1:704041/1‑93 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           >  1:4382659/1‑93 (MQ=255)
                                               aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                          >  2:429813/1‑93 (MQ=255)
                                               aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                          <  1:2268882/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAACAATATc                                        >  2:2389030/1‑93 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        >  2:4144568/1‑93 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  2:2199109/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  2:2067172/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  1:4000726/93‑1 (MQ=255)
                                                     aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                    <  1:3335496/93‑1 (MQ=255)
                                                     aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                    <  2:3304226/93‑1 (MQ=255)
                                                     aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                    >  1:752028/1‑93 (MQ=255)
                                                      ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                   >  1:346789/1‑93 (MQ=255)
                                                      ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                   >  2:1192509/1‑93 (MQ=255)
                                                       tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                  <  1:2379223/93‑1 (MQ=255)
                                                        gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                 <  1:3197846/93‑1 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             <  2:137538/93‑1 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             <  2:3672016/93‑1 (MQ=255)
                                                              tGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                           >  2:2629488/1‑93 (MQ=255)
                                                               gACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                          <  1:1781837/93‑1 (MQ=255)
                                                                aCCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                         <  2:406735/93‑1 (MQ=255)
                                                                 ccAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCCCAAGCGAGGGTAAAAATATCGTAAAAAGCCAGTCCa                        >  2:4073729/1‑93 (MQ=255)
                                                                 acAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                        <  1:1703758/92‑1 (MQ=255)
                                                                  cAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       <  1:1606904/93‑1 (MQ=255)
                                                                  cAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       >  2:662159/1‑93 (MQ=255)
                                                                   aGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                      >  1:1396609/1‑93 (MQ=255)
                                                                   aGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                      >  1:85035/1‑93 (MQ=255)
                                                                    ggCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                     >  1:3927170/1‑93 (MQ=255)
                                                                    ggCGCAGAGCTTTACGCCTGATGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                     <  1:3127619/93‑1 (MQ=255)
                                                                      cgcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   <  2:3403147/93‑1 (MQ=255)
                                                                      cgcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   <  2:3463662/93‑1 (MQ=255)
                                                                         agagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                <  1:543132/93‑1 (MQ=255)
                                                                          gagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg               >  2:2754048/1‑93 (MQ=255)
                                                                           agCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa              <  1:3884236/93‑1 (MQ=255)
                                                                            gCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc             <  2:3745409/93‑1 (MQ=255)
                                                                            gCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc             <  1:2344732/93‑1 (MQ=255)
                                                                             cTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg            <  2:4091244/93‑1 (MQ=255)
                                                                              tttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa           <  2:2837248/93‑1 (MQ=255)
                                                                                tACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc         <  1:3463120/93‑1 (MQ=255)
                                                                                tACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc         <  2:1804164/93‑1 (MQ=255)
                                                                                 aCGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa        <  1:2966234/93‑1 (MQ=255)
                                                                                 aCGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa        <  1:3671760/93‑1 (MQ=255)
                                                                                   gCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt      >  1:379798/1‑93 (MQ=255)
                                                                                   gCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt      <  1:2576158/93‑1 (MQ=255)
                                                                                    ccTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAACAGACAGTTCACCTTTACGCACGATCAtct     >  1:2156410/1‑91 (MQ=255)
                                                                                     cTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt    <  2:3762627/93‑1 (MQ=255)
                                                                                     cTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt    <  1:3447616/93‑1 (MQ=255)
                                                                                     cTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt    >  1:4261316/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   <  2:3924042/93‑1 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  1:4521216/1‑93 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt  >  2:2980806/1‑93 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt  <  1:827854/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081

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