Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 588,320 C→T A879T (GCC→ACC)  nfrA ← bacteriophage N4 receptor, outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913588,3200CT99.6% 976.7 / ‑4.6 243A879T (GCC→ACC) nfrAbacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (118/124);  minor base A (1/0);  total (119/124)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 243 total aligned reads displayed.

GTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGCAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATT  >  NC_000913/588222‑588416
                                                                                                  |                                                                                                
gTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTaa                                                                                                >  2:781783/1‑101 (MQ=255)
 ttAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTaaa                                                                                               >  1:781784/1‑101 (MQ=255)
 ttAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTaaa                                                                                               <  1:781466/101‑1 (MQ=255)
   aGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAAcc                                                                                             >  1:781787/1‑101 (MQ=255)
     cTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAg                                                                                           >  2:781789/1‑101 (MQ=255)
      tACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGc                                                                                          >  1:781794/1‑101 (MQ=255)
         cTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTg                                                                                       >  1:781798/1‑101 (MQ=255)
         cTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTg                                                                                       >  2:781797/1‑101 (MQ=255)
          ttCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCa                                                                                            <  1:781572/95‑1 (MQ=255)
          ttCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGa                                                                                      >  2:781801/1‑101 (MQ=255)
            cTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAAc                                                                                    >  1:781802/1‑101 (MQ=255)
            cTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAAc                                                                                    <  1:781454/101‑1 (MQ=255)
                tGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGtt                                                                                <  1:781427/101‑1 (MQ=255)
                tGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGtt                                                                                >  2:781806/1‑101 (MQ=255)
                 gCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTc                                                                               <  1:781173/101‑1 (MQ=255)
                 gCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTc                                                                               <  2:781509/101‑1 (MQ=255)
                 gCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTc                                                                               <  1:781377/101‑1 (MQ=255)
                 gCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTc                                                                               <  1:781342/101‑1 (MQ=255)
                  ccAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCg                                                                              >  1:781809/1‑101 (MQ=255)
                   cAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCgg                                                                             >  2:781811/1‑101 (MQ=255)
                   cAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCgg                                                                             <  2:781383/101‑1 (MQ=255)
                    aGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCggg                                                                            >  1:781812/1‑101 (MQ=255)
                     gCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGt                                                                           <  1:781435/101‑1 (MQ=255)
                     gCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGt                                                                           <  1:781314/101‑1 (MQ=255)
                      cTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCggg                                                                            >  2:781815/1‑99 (MQ=255)
                        gACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGcc                                                                        >  1:781817/1‑101 (MQ=255)
                          cgcgATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCAt                                                                      <  2:781601/101‑1 (MQ=255)
                           gcgATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCAtt                                                                     <  2:781305/101‑1 (MQ=255)
                           gcgATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCAtt                                                                     >  1:781821/1‑101 (MQ=255)
                            cgATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtc                                                                    <  1:781392/101‑1 (MQ=255)
                            cgATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGCACAGGTTTTGGGTAAAGCAGCCTGAACCGTTCGGGTGCCATtc                                                                    >  2:781823/1‑101 (MQ=255)
                              aTAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtcgt                                                                  <  2:781442/101‑1 (MQ=255)
                              aTAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtcgt                                                                  <  2:781458/101‑1 (MQ=255)
                              aTAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtcgt                                                                  <  2:781479/101‑1 (MQ=255)
                               tAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtcgtc                                                                 >  2:781826/1‑101 (MQ=255)
                               tAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATtcgtc                                                                 >  1:781827/1‑101 (MQ=255)
                                  tCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCt                                                              >  1:781831/1‑101 (MQ=255)
                                    tGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGt                                                            <  2:781559/101‑1 (MQ=255)
                                     gCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTa                                                           <  2:781553/101‑1 (MQ=255)
                                      ccGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTAt                                                          >  1:781836/1‑101 (MQ=255)
                                      ccGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTAt                                                          <  2:781421/101‑1 (MQ=255)
                                      ccGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTAt                                                          >  2:781835/1‑101 (MQ=255)
                                      ccGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTAt                                                          <  2:781651/101‑1 (MQ=255)
                                         tCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTg                                                       <  2:781406/101‑1 (MQ=255)
                                          ccACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgc                                                      >  2:781839/1‑101 (MQ=255)
                                          ccACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgc                                                      <  1:781495/101‑1 (MQ=255)
                                           cACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgcc                                                     >  2:781841/1‑101 (MQ=255)
                                            aCGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgccg                                                    >  1:781842/1‑101 (MQ=255)
                                            aCGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgccg                                                    <  1:781472/101‑1 (MQ=255)
                                            aCGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgccg                                                    <  2:781568/101‑1 (MQ=255)
                                              gCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgccgcc                                                  <  2:781476/101‑1 (MQ=255)
                                               ccTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCa                                                 >  2:781846/1‑101 (MQ=255)
                                                cTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAAc                                                                                    <  1:781492/65‑1 (MQ=255)
                                                 tGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCAtt                                               >  2:781848/1‑101 (MQ=255)
                                                    atatCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTaaa                                            >  2:781852/1‑101 (MQ=255)
                                                     tatCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAg                                           <  2:781358/101‑1 (MQ=255)
                                                      atCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGa                                          >  1:781853/1‑101 (MQ=255)
                                                      atCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGa                                          <  2:781397/101‑1 (MQ=255)
                                                      atCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGa                                          <  1:781481/101‑1 (MQ=255)
                                                         cTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATg                                       >  2:781856/1‑101 (MQ=255)
                                                           ggCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAg                                     <  1:781538/101‑1 (MQ=255)
                                                           ggCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAg                                     <  2:781515/101‑1 (MQ=255)
                                                             cGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgc                                   <  2:781540/101‑1 (MQ=255)
                                                              ggATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcg                                  >  1:781863/1‑101 (MQ=255)
                                                              ggATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcg                                  <  1:781699/101‑1 (MQ=255)
                                                              ggATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcg                                  >  1:781861/1‑101 (MQ=255)
                                                              ggATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcg                                  >  2:781860/1‑101 (MQ=255)
                                                               gATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcgc                                 <  2:781525/101‑1 (MQ=255)
                                                               gATATATTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGgcgc                                 <  1:781446/101‑1 (MQ=255)
                                                                atatatTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCt                                <  2:781410/101‑1 (MQ=255)
                                                                atatatTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCt                                >  1:781867/1‑101 (MQ=255)
                                                                 tatatTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTg                               <  1:781521/101‑1 (MQ=255)
                                                                  atatTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAg                                     >  2:781868/1‑94 (MQ=255)
                                                                   tatTGCGCCGCATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGgc                             >  1:781870/1‑101 (MQ=255)
                                                                        cgccgcATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCgca                        >  2:781873/1‑101 (MQ=255)
                                                                        cgccgcATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCgca                        <  2:781580/101‑1 (MQ=255)
                                                                            gcATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCAt                    <  2:781501/101‑1 (MQ=255)
                                                                            gcATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCAt                    <  1:781847/101‑1 (MQ=255)
                                                                            gcATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCAt                    >  2:781877/1‑101 (MQ=255)
                                                                             cATCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATg                   >  1:781878/1‑101 (MQ=255)
                                                                              aTCGAGGTACAGGTTTTGGGTAACCCAGCCTGAACCGTTCgg                                                                             <  2:781226/42‑1 (MQ=255)
                                                                               tCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGt                 <  1:781627/101‑1 (MQ=255)
                                                                               tCGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGt                 <  1:781432/101‑1 (MQ=255)
                                                                                cGAGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTa                >  1:781881/1‑101 (MQ=255)
                                                                                 gatgTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTAt               >  1:781883/4‑101 (MQ=255)
                                                                                  aGGTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATc              <  1:781297/101‑1 (MQ=255)
                                                                                   ggTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCg             <  1:781584/101‑1 (MQ=255)
                                                                                   ggTACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCg             >  2:781885/1‑101 (MQ=255)
                                                                                     tACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGa           <  2:781661/101‑1 (MQ=255)
                                                                                     tACAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGa           >  1:781889/1‑101 (MQ=255)
                                                                                      aCAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGAt          >  2:781891/1‑101 (MQ=255)
                                                                                       cAGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATg         <  2:781557/101‑1 (MQ=255)
                                                                                        aGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATgc        >  1:781893/1‑101 (MQ=255)
                                                                                        aGGTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATgc        >  2:781894/1‑101 (MQ=255)
                                                                                         ggTTTTGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTgccg                                                    >  2:781898/1‑56 (MQ=255)
                                                                                           ttttGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGcc     <  2:781504/101‑1 (MQ=255)
                                                                                           ttttGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGcc     <  2:781463/101‑1 (MQ=255)
                                                                                            tttGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCa    <  1:781363/101‑1 (MQ=255)
                                                                                            tttGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCa    >  2:781901/1‑101 (MQ=255)
                                                                                              tGGGTAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCAtt  <  2:781434/101‑1 (MQ=255)
                                                                                                  |                                                                                                
GTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAATATCCTGGCGGATATATTGCGCCGCATCGAGGTACAGGTTTTGGGCAAACCAGCCTGAACCGTTCGGGTGCCATTCGTCGCTGTATTTGCCGCCATTAAAGAATGAGGCGCTGGCGCGCAGCATGGTATCGGATGCGCCATT  >  NC_000913/588222‑588416

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