Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 13.8% intergenic (+20/+29) slyX → / ← slyD phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,8780TA13.8% 251.3 / 29.9 87intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (5/7);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGT  >  NC_000913/3477734‑3478025
                                                                                                                                                |                                                                                                                                                    
ccGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTccc                                                                                                                                                  <  2:836239/149‑1 (MQ=255)
   cTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGc                                                                                                                                               <  4:1812153/149‑1 (MQ=255)
     cATGAAATGGAGATGGCGATACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            <  2:1351564/149‑1 (MQ=255)
     cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tt                                                                                                                                            >  1:919830/1‑149 (MQ=255)
     cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            <  2:1717705/149‑1 (MQ=255)
                   ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACttttt                                                                                                                              >  4:1150389/1‑149 (MQ=255)
                    gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTc                                                                                                                             <  2:919830/149‑1 (MQ=255)
                        aaCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCg                                                                                                                            >  1:419768/1‑146 (MQ=255)
                        aaCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCg                                                                                                                            <  2:419768/146‑1 (MQ=255)
                        aaCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTa                                                                                                                         <  3:1150389/149‑1 (MQ=255)
                             agcgATCATCTGCGTCTGCTGACCGAGATGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTAGGTATTAGt                                                                                                                    <  4:1628472/148‑1 (MQ=255)
                                gATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                 >  2:1563057/1‑149 (MQ=255)
                                gATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                 >  2:174210/1‑149 (MQ=255)
                                     tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaaccg                                                                                                            >  4:1058059/1‑149 (MQ=255)
                                     tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            <  2:914494/117‑1 (MQ=255)
                                     tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            >  1:914494/1‑117 (MQ=255)
                                     tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            <  4:297261/117‑1 (MQ=255)
                                     tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                                                                            >  3:297261/1‑117 (MQ=255)
                                       tGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                  >  2:1277844/1‑141 (MQ=255)
                                       tGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                  <  1:1277844/141‑1 (MQ=255)
                                               ctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGcaaccgcaaccg                                                                                                      <  1:424157/145‑1 (MQ=255)
                                               ctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGcaaccgcaaccg                                                                                                      >  2:424157/1‑145 (MQ=255)
                                                    cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                               <  3:1228094/147‑1 (MQ=255)
                                                    cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                               >  4:1228094/1‑147 (MQ=255)
                                                    cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCt                                                                                             <  3:1058059/149‑1 (MQ=255)
                                                     gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTa                                                                                                                      <  1:1239825/123‑1 (MQ=255)
                                                     gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTa                                                                                                                      >  2:1239825/1‑123 (MQ=255)
                                                       gaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTa                                                                                          >  3:1258931/1‑149 (MQ=255)
                                                           ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACt                                                                                                                                  <  3:1404336/105‑1 (MQ=255)
                                                           ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACt                                                                                                                                  >  4:1404336/1‑105 (MQ=255)
                                                             aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                    >  3:837143/1‑149 (MQ=255)
                                                               aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGt                                                                                                                    <  2:213509/115‑1 (MQ=255)
                                                               aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGt                                                                                                                    >  1:213509/1‑115 (MQ=255)
                                                                agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑t                                                                                                                                             <  3:822808/89‑1 (MQ=255)
                                                                agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑t                                                                                                                                             >  4:822808/1‑89 (MQ=255)
                                                                agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                                                                 <  1:174210/149‑1 (MQ=255)
                                                                 gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagc                                                                                >  1:1010273/1‑149 (MQ=255)
                                                                   cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagcag                                                                              >  3:40136/1‑149 (MQ=255)
                                                                            gTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑ttt                                                                                                                                           >  4:684006/1‑79 (MQ=255)
                                                                            gTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑ttt                                                                                                                                           <  3:684006/79‑1 (MQ=255)
                                                                            gTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCa                                                                     >  1:1502968/1‑149 (MQ=255)
                                                                                   aTCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt                                                              <  2:1502968/149‑1 (MQ=255)
                                                                                       cgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATg                                                          >  2:1206740/1‑149 (MQ=255)
                                                                                        gTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGa                                                         <  1:1563057/149‑1 (MQ=255)
                                                                                              ggCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGa                                                   <  4:1258931/149‑1 (MQ=255)
                                                                                               gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcc                                                                                                                                                   <  2:901795/53‑1 (MQ=255)
                                                                                               gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcc                                                                                                                                                   >  1:901795/1‑53 (MQ=255)
                                                                                               gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGcc                                                                                              >  4:1257358/1‑105 (MQ=255)
                                                                                               gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTGTTCGGTATTAGTGGCAACCGCAACCGCCGTTTcc                                                                                              <  3:1257358/105‑1 (MQ=255)
                                                                                                    agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             <  4:40136/149‑1 (MQ=255)
                                                                                                    agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             <  1:1206740/149‑1 (MQ=255)
                                                                                                    agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATc                                                                                                                                                    <  3:315937/47‑1 (MQ=255)
                                                                                                     gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                 <  2:1260358/80‑1 (MQ=255)
                                                                                                     gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                 >  1:1260358/1‑80 (MQ=255)
                                                                                                         cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             >  4:198234/1‑144 (MQ=255)
                                                                                                         cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             >  2:1147249/1‑144 (MQ=255)
                                                                                                         cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             <  3:198234/144‑1 (MQ=255)
                                                                                                         cgccaccgacacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             <  1:1147249/144‑1 (MQ=255)
                                                                                                           ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt                                                              >  1:1325544/1‑125 (MQ=255)
                                                                                                           ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt                                                              <  2:1325544/125‑1 (MQ=255)
                                                                                                                 ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGcc                                                                            <  4:475499/105‑1 (MQ=255)
                                                                                                                 ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGcc                                                                            >  3:126191/1‑105 (MQ=255)
                                                                                                                 ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGcc                                                                            >  3:475499/1‑105 (MQ=255)
                                                                                                                 ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGcc                                                                            <  4:126191/105‑1 (MQ=255)
                                                                                                                     attattGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCgcagca                            >  1:731283/1‑149 (MQ=255)
                                                                                                                                gacaaaaaaGCGGGGAATCCCCGCTTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAcc                                     >  1:1516786/4‑130 (MQ=255)
                                                                                                                                gacaaaaaaGCGGGGAATCCCCGCTTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAcc                                     <  2:1516786/127‑1 (MQ=255)
                                                                                                                                gTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGa                                                   <  4:863207/115‑1 (MQ=255)
                                                                                                                                gTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGa                                                   >  3:863207/1‑115 (MQ=255)
                                                                                                                                gTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                            <  2:763656/122‑1 (MQ=255)
                                                                                                                                gTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                            >  1:763656/1‑122 (MQ=255)
                                                                                                                                 tAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtc                >  1:1579992/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGCGTGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGAAGCAACCgtcgtggtcgt              >  1:118811/1‑149 (MQ=255)
                                                                                                                                    aaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             <  1:1480913/117‑1 (MQ=255)
                                                                                                                                    aaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                             >  2:1480913/1‑117 (MQ=255)
                                                                                                                                     aaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGa            <  2:1579992/149‑1 (MQ=255)
                                                                                                                                      aaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGAt           <  2:1010273/149‑1 (MQ=255)
                                                                                                                                       aaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc          >  2:421148/1‑149 (MQ=255)
                                                                                                                                         gCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgt        >  3:860791/1‑149 (MQ=255)
                                                                                                                                             ggATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg    >  2:1570969/1‑149 (MQ=255)
                                                                                                                                             ggATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg    >  4:745818/1‑149 (MQ=255)
                                                                                                                                             ggAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg    >  4:1694797/1‑149 (MQ=255)
                                                                                                                                              gATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtgg                  >  1:961939/1‑134 (MQ=255)
                                                                                                                                              gATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtgg                  <  2:961939/134‑1 (MQ=255)
                                                                                                                                               aTTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggt  <  4:837143/149‑1 (MQ=255)
                                                                                                                                               aTTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggt     >  2:118838/1‑146 (MQ=255)
                                                                                                                                               aTTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggt     <  1:118838/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGT  >  NC_000913/3477734‑3478025

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