Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 10.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,8780TA10.8% 246.6 / 18.5 84intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base A (3/6);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGG  >  NC_000913/3477733‑3478024
                                                                                                                                                 |                                                                                                                                                  
aCCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcc                                                                                                                                                 <  4:488041/149‑1 (MQ=255)
 ccGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTccc                                                                                                                                                >  1:383765/1‑149 (MQ=255)
  cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcccc                                                                                                                                               <  2:593971/149‑1 (MQ=255)
    cTCGTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGCCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                                             <  1:231876/149‑1 (MQ=255)
      cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                           <  4:115401/149‑1 (MQ=255)
       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttt                                                                                                                                          >  3:444393/1‑149 (MQ=255)
       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttt                                                                                                                                          <  3:136654/149‑1 (MQ=255)
       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttt                                                                                                                                          >  2:626351/1‑149 (MQ=255)
       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttt                                                                                                                                          <  4:503715/149‑1 (MQ=255)
       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttt                                                                                                                                          <  4:259981/149‑1 (MQ=255)
         gAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                        <  3:35072/149‑1 (MQ=255)
         gAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCACCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                        >  3:239611/1‑149 (MQ=255)
          aaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGGGGCGTAAAAAAAGCGGGGAATCCCCGCtttttt                                                                                                                                       <  1:71540/149‑1 (MQ=255)
          aaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCGCGCtttttt                                                                                                                                       <  4:587196/149‑1 (MQ=255)
             tGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCACCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTc                                                                                                                                    <  1:454293/149‑1 (MQ=255)
                agaTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACt                                                                                                                                 >  1:39811/1‑149 (MQ=255)
                 gaTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACtt                                                                                                                                >  1:455281/1‑149 (MQ=255)
                       gacagtGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAATGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCgg                                                                                                                          >  2:300478/6‑149 (MQ=255)
                           cTGCGCGATCATCTGCCTCTGCTGAGCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAGACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAtt                                                                                                                      <  1:527642/149‑1 (MQ=255)
                            tGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTa                                                                                                                     <  1:300478/149‑1 (MQ=255)
                                    catcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaac                                                                                                             >  3:365018/1‑149 (MQ=255)
                                               gctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                           <  1:312055/108‑1 (MQ=255)
                                               gctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                           >  2:312055/1‑108 (MQ=255)
                                                   aCCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                              <  2:600403/149‑1 (MQ=255)
                                                    ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGcc                                                                                             <  4:361186/149‑1 (MQ=255)
                                                      gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAACGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCtt                                                                                           >  4:386240/1‑149 (MQ=255)
                                                      gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCCTCGCAGGCTGAAGAAACACCACCGCCACATGATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTCTTTTGTCATTTTTTCGGCATTAGTGGCAACAGCAACCGCCGTTGCCtt                                                                                           >  2:149007/1‑149 (MQ=255)
                                                         aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGTCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                                                                        <  4:490104/149‑1 (MQ=255)
                                                         aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                                                                        <  4:244396/149‑1 (MQ=255)
                                                           gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgccg                                                                                                  <  2:475551/137‑1 (MQ=255)
                                                           gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgccg                                                                                                  >  1:475551/1‑137 (MQ=255)
                                                           gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCg                                                                                      >  1:283757/1‑149 (MQ=255)
                                                            ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                     >  2:77539/1‑149 (MQ=255)
                                                            ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                     >  2:282285/1‑149 (MQ=255)
                                                            ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                     >  3:565610/1‑149 (MQ=255)
                                                            ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                     >  4:561482/1‑149 (MQ=255)
                                                              aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                   >  3:224418/1‑149 (MQ=255)
                                                               aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                                                                  >  1:73387/1‑149 (MQ=255)
                                                                 agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgc                                                                                                          >  4:667486/1‑123 (MQ=255)
                                                                 agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgc                                                                                                          <  3:667486/123‑1 (MQ=255)
                                                                    cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagcag                                                                             >  2:409666/1‑149 (MQ=255)
                                                                         agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCTCTACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTCTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCATAGCAGCCTTc                                                                        <  4:444393/149‑1 (MQ=255)
                                                                         agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTc                                                                        <  3:289275/149‑1 (MQ=255)
                                                                         agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTc                                                                        <  3:465512/149‑1 (MQ=255)
                                                                            cGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGcc                                                                     >  4:364728/1‑149 (MQ=255)
                                                                                 aaCATGGCGTCGCAGGCTGAAGAAACGACACCGCCACATTATTGAGGCGTAAAAAAAGCGGGAAGTCCCCGCTCTTGTGACACTTTTGCGGTATTAGTGGCAACCGCGACCGCCGTTGCCTGTACCGCCACAGCAGCCTTCACCACCgt                                                                >  2:503542/1‑149 (MQ=255)
                                                                                 aaCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                >  3:200146/1‑149 (MQ=255)
                                                                                 aaCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                <  3:89339/149‑1 (MQ=255)
                                                                                        cgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATg                                                         >  2:468893/1‑149 (MQ=255)
                                                                                           cGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGAcc                                                      <  3:34253/149‑1 (MQ=255)
                                                                                                  tGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                <  2:692849/84‑1 (MQ=255)
                                                                                                  tGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                >  1:692849/1‑84 (MQ=255)
                                                                                                   gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                                                                        <  3:486370/107‑1 (MQ=255)
                                                                                                   gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                                                                        >  4:486370/1‑107 (MQ=255)
                                                                                                    aagaaACGCCACCGCCACATTACTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTg                                             <  3:418945/149‑1 (MQ=255)
                                                                                                      gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCTCTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                                                                  <  1:442872/110‑1 (MQ=255)
                                                                                                      gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                                                                  >  2:442872/1‑110 (MQ=255)
                                                                                                      gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                           <  1:282285/149‑1 (MQ=255)
                                                                                                      gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                           <  1:77539/149‑1 (MQ=255)
                                                                                                         aCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCAt                                        >  2:226050/1‑149 (MQ=255)
                                                                                                           gccaccgccatATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTCGTCACTTCTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGa                                      <  1:503542/149‑1 (MQ=255)
                                                                                                            ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAc                                     <  1:468893/149‑1 (MQ=255)
                                                                                                                cgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGGCCACGATCACGGTCATGaacacg                                 <  4:60978/149‑6 (MQ=255)
                                                                                                                 gccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAccgccg                                >  1:283946/1‑149 (MQ=255)
                                                                                                                  ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACcgccgc                               >  4:462719/1‑149 (MQ=255)
                                                                                                                    acaTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCgcag                             >  3:224726/1‑149 (MQ=255)
                                                                                                                          ttGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTg                                             <  4:118085/126‑1 (MQ=255)
                                                                                                                          ttGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTg                                             >  3:118085/1‑126 (MQ=255)
                                                                                                                            gAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtc                     <  2:206611/149‑1 (MQ=255)
                                                                                                                               gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtg                  >  4:125254/1‑149 (MQ=255)
                                                                                                                               gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtg                  >  1:464053/1‑149 (MQ=255)
                                                                                                                                  tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtc               <  2:429429/149‑1 (MQ=255)
                                                                                                                                      aaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATTGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGa           <  1:654558/149‑1 (MQ=255)
                                                                                                                                       aaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGTCCGCCGCAGCAACCGTCGTGGTCGTGAt          >  1:427548/1‑149 (MQ=255)
                                                                                                                                       aaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTTGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGAt          >  2:610728/1‑149 (MQ=255)
                                                                                                                                        aaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc         <  1:626351/149‑1 (MQ=255)
                                                                                                                                        aaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTTCCGCCACAGCAGCCTTCGCCACCGTGTTCATGa                                                        <  2:561327/102‑1 (MQ=255)
                                                                                                                                        aaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGa                                                        >  1:561327/1‑102 (MQ=255)
                                                                                                                                          gCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgt       <  3:364728/149‑1 (MQ=255)
                                                                                                                                            ggggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtgg     <  1:303028/149‑1 (MQ=255)
                                                                                                                                            ggggATTCCCCGCTGTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtgg     <  3:378818/149‑1 (MQ=255)
                                                                                                                                              ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg   >  3:12407/1‑149 (MQ=255)
                                                                                                                                              ggATACCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg   <  1:313654/149‑1 (MQ=255)
                                                                                                                                               gATTTCCCGCTTTTTTGTCACTTTTTAGGTATTAGTGGCAACCGCAACTGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtgg  <  2:199131/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGG  >  NC_000913/3477733‑3478024

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