Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,286,677 A→G T75A (ACT→GCT)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,286,6770AG99.4% 1047.2 / ‑5.3 331T75A (ACT→GCT) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (152/177);  total (153/178)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

TCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCA  >  NC_000913/3286588‑3286766
                                                                                         |                                                                                         
 cGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTa                                                                                       >  1:4863504/1‑93 (MQ=255)
   aCCCCATAATGCTGTTACCCACAACCAGTAACATGCCAACAAACGAGCTGTAGCGCTGTCTGAAACAATATGCACAACTACCAGGAGCTCAtt                                                                                     >  2:3686318/1‑92 (MQ=255)
    ccccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCa                                                                                    >  1:764594/1‑93 (MQ=255)
    ccccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCa                                                                                    <  2:4206938/93‑1 (MQ=255)
     cccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAg                                                                                   <  1:6508663/93‑1 (MQ=255)
     cccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAg                                                                                   >  2:2222352/1‑93 (MQ=255)
     cccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAg                                                                                   >  1:5718488/1‑93 (MQ=255)
      ccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAgg                                                                                  <  2:5971736/93‑1 (MQ=255)
      ccATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGAACAACTACCAGGAGCTTATCAgg                                                                                  >  1:1905659/1‑93 (MQ=255)
       cATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGt                                                                                 <  2:5797447/93‑1 (MQ=255)
       cATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGt                                                                                 <  2:5740172/93‑1 (MQ=255)
           aTGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGat                                                                             >  2:1509918/1‑93 (MQ=255)
            tGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTCCCAGGAGCTTATCAGGTTGata                                                                            <  1:3994039/93‑1 (MQ=255)
              cTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATc                                                                          >  2:4056013/1‑93 (MQ=255)
              cTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATc                                                                          <  2:5200900/93‑1 (MQ=255)
               tGGAAACCACAAAAAGTGGTATGCAAACAGCCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCt                                                                         <  1:497573/93‑1 (MQ=255)
               tGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCt                                                                         >  2:6928981/1‑93 (MQ=255)
               tGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCt                                                                         <  1:148994/93‑1 (MQ=255)
               tGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCt                                                                         <  1:7280408/93‑1 (MQ=255)
                ggAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTg                                                                        <  1:5517935/93‑1 (MQ=255)
                ggAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTg                                                                        <  1:2514814/93‑1 (MQ=255)
                  aaaCCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGc                                                                      <  2:5017168/93‑1 (MQ=255)
                    aCCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTg                                                                    <  2:2692690/93‑1 (MQ=255)
                    aCCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTg                                                                    >  2:776623/1‑93 (MQ=255)
                      cacaAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGaa                                                                  <  2:1868207/93‑1 (MQ=255)
                      cacaAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGaa                                                                  <  1:1121383/93‑1 (MQ=255)
                      cacaAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGaa                                                                  <  2:3803424/93‑1 (MQ=255)
                       acaAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAAt                                                                 <  1:444566/93‑1 (MQ=255)
                        caAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATa                                                                >  1:4845119/1‑93 (MQ=255)
                        caAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATa                                                                >  1:5686637/1‑93 (MQ=255)
                        aaaaaaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATa                                                                <  1:5567882/92‑1 (MQ=255)
                         aaaaaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaa                                                               >  2:6707776/1‑93 (MQ=255)
                         aaaaaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaa                                                               <  2:5287794/93‑1 (MQ=255)
                          aaaaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaaa                                                              <  1:1052083/93‑1 (MQ=255)
                          aaaaGTGGTATGCAAACAACCGACCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaaa                                                              <  2:5641108/93‑1 (MQ=255)
                           aaaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaaaa                                                             <  1:1682856/93‑1 (MQ=255)
                            aaGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATaaaaa                                                            <  2:300155/93‑1 (MQ=255)
                             aGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAaag                                                           <  1:646673/93‑1 (MQ=255)
                              gTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAaaga                                                          <  1:5090403/93‑1 (MQ=255)
                                ggTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAaagaag                                                        >  1:2530581/1‑93 (MQ=255)
                                ggTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAaagaag                                                        <  2:3801294/93‑1 (MQ=255)
                                   aTGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGtt                                                     <  1:3072005/93‑1 (MQ=255)
                                   aTGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGtt                                                     <  1:3185139/93‑1 (MQ=255)
                                    tGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGttt                                                    >  1:3390819/1‑93 (MQ=255)
                                    tGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGttt                                                    <  1:1839057/93‑1 (MQ=255)
                                     gCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTCAATAAAAAGAAGGTTTc                                                   <  2:5465524/93‑1 (MQ=255)
                                      cAAACGACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTca                                                  <  1:4743550/93‑1 (MQ=255)
                                        aacaacCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTcaca                                                >  2:1412951/1‑93 (MQ=255)
                                         acaacCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAg                                               <  2:6577415/93‑1 (MQ=255)
                                         acaacCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAg                                               <  1:4211946/93‑1 (MQ=255)
                                           aacCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGaa                                             >  1:6235172/1‑93 (MQ=255)
                                            acCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGaaa                                            <  1:6003382/93‑1 (MQ=255)
                                            acCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGaaa                                            >  1:1269597/1‑93 (MQ=255)
                                               gATCTGTTGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGaaaaaa                                         >  1:2196805/1‑93 (MQ=255)
                                                 tCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAAt                                       <  2:7396619/93‑1 (MQ=255)
                                                  cTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAAtt                                      <  1:745724/93‑1 (MQ=255)
                                                   tGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTa                                     <  2:730776/93‑1 (MQ=255)
                                                    gTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTAc                                    <  1:4640727/93‑1 (MQ=255)
                                                     tCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACa                                   <  2:3381069/93‑1 (MQ=255)
                                                      cgcgTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACaa                                  <  2:3197634/93‑2 (MQ=255)
                                                       gcgTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACAtt                                 <  2:5306068/93‑1 (MQ=255)
                                                       gcgTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACAtt                                 <  2:3888668/93‑1 (MQ=255)
                                                       gcgTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACAtt                                 <  2:4577238/93‑1 (MQ=255)
                                                       gcgTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACAtt                                 <  1:3750792/93‑1 (MQ=255)
                                                         gTTTTTCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTa                               >  2:2763632/1‑93 (MQ=255)
                                                          tttttCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAc                              >  1:895688/1‑93 (MQ=255)
                                                          tttttCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAc                              >  1:4515965/1‑93 (MQ=255)
                                                          tttttCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAc                              <  1:6943447/93‑1 (MQ=255)
                                                           ttttCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAcc                             <  1:1746704/93‑1 (MQ=255)
                                                             ttCAAAAAAATATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGGATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTAACGc                           >  2:4497229/1‑93 (MQ=255)
                                                                aaaaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCaa                        <  1:2468495/93‑1 (MQ=255)
                                                                aaaaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCaa                        >  1:7757270/1‑93 (MQ=255)
                                                                 aaaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAAt                       >  1:5181747/1‑93 (MQ=255)
                                                                 aaaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAAt                       <  1:5969424/93‑1 (MQ=255)
                                                                 aaaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAAt                       >  1:537325/1‑93 (MQ=255)
                                                                  aaaaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATg                      >  2:4668869/1‑93 (MQ=255)
                                                                    aaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGcg                    >  2:1668918/1‑92 (MQ=255)
                                                                    aaaTATGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCa                    <  2:3093345/93‑1 (MQ=255)
                                                                      atatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCaga                  <  1:7227245/93‑1 (MQ=255)
                                                                      atatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCaga                  <  2:2002263/93‑1 (MQ=255)
                                                                      atatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCaga                  <  1:7097252/93‑1 (MQ=255)
                                                                      atatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCaga                  <  2:6302305/93‑1 (MQ=255)
                                                                      atatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCaga                  >  1:3540767/1‑93 (MQ=255)
                                                                       tatGCACAACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCagag                 >  2:1269879/1‑93 (MQ=255)
                                                                           cacaACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAAc             <  2:3513799/93‑1 (MQ=255)
                                                                            gcaACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACt            <  1:6059524/92‑1 (MQ=255)
                                                                            acaACTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACt            >  1:6793905/1‑93 (MQ=255)
                                                                              aaCTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTc          <  1:7907097/93‑1 (MQ=255)
                                                                                cTACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTg        >  1:5004954/1‑93 (MQ=255)
                                                                                 tACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGc       >  1:2088871/1‑93 (MQ=255)
                                                                                 tACCAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAAATACATTTACCGCCAATGCAGAGCAACTTCTGc       <  2:6235171/93‑1 (MQ=255)
                                                                                    cAGGAGCTTATCAGGTTGATATCTGGCTTAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGc    <  1:7555328/93‑1 (MQ=255)
                                                                                    cAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGc    >  1:2703391/1‑93 (MQ=255)
                                                                                    cAGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGAAGc    >  1:1950889/1‑93 (MQ=255)
                                                                                     aGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGcc   >  1:1402282/1‑93 (MQ=255)
                                                                                     aGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGcc   >  2:7571567/1‑93 (MQ=255)
                                                                                     aGGAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGcc   <  1:389858/93‑1 (MQ=255)
                                                                                      ggAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCca  >  1:20549/1‑93 (MQ=255)
                                                                                      ggAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCca  >  2:3259553/1‑93 (MQ=255)
                                                                                      ggAGCTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCca  >  2:7891677/1‑93 (MQ=255)
                                                                                         |                                                                                         
TCGACCCCATTATGCTGGAAACCACAAAAAGTGGTATGCAAACAACCGATCTGTCGCGTTTTTCAAAAAAATATGCACAACTACCAGGAACTTATCAGGTTGATATCTGGCTGAATAAAAAGAAGGTTTCACAGAAAAAAATTACATTTACCGCCAATGCAGAGCAACTTCTGCAGCCA  >  NC_000913/3286588‑3286766

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: ‑