Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,760,589 C→A 5.2% R15R (CGG→AGG)  yfjJ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,760,5890CA 5.2% 558.5 / 15.2 155R15R (CGG→AGG) yfjJCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (73/74);  new base A (4/4);  total (77/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 155 total aligned reads displayed.

GGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCT  >  NC_000913/2760455‑2760735
                                                                                                                                      |                                                                                                                                                  
        gCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTa                                                                                                                               >  2:249081/1‑148 (MQ=255)
         cTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                              <  1:249081/148‑1 (MQ=255)
             tCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       <  4:989573/135‑1 (MQ=255)
              cATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                         >  4:984912/1‑148 (MQ=255)
                tCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGAc                                                                                                                       >  1:549506/1‑148 (MQ=255)
                           aaTATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                                                                                                                                   >  3:1023609/1‑109 (MQ=255)
                           aaTATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                                                                                                                                   <  4:1023609/109‑1 (MQ=255)
                                   ggTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTc                                                                                                                     <  2:10253/131‑1 (MQ=255)
                                       cATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                              <  4:650999/118‑1 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTTCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                            >  2:488851/1‑118 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTTCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                            <  1:488851/118‑1 (MQ=255)
                                         tAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGaa                                                                                                                                           <  3:99170/103‑1 (MQ=255)
                                          aGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAAt                                                                                             >  3:188116/1‑148 (MQ=255)
                                           gACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  2:172346/1‑117 (MQ=255)
                                            acacATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                              <  1:733028/145‑1 (MQ=255)
                                            acacATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                              >  2:733028/1‑145 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  3:649070/1‑110 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           <  4:649070/110‑1 (MQ=255)
                                                  cAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTAc                                                                                     >  1:625202/1‑148 (MQ=255)
                                                    gAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       <  3:1134946/96‑1 (MQ=255)
                                                    gAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       >  4:1134946/1‑96 (MQ=255)
                                                      aCGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTaa                                                                                                   >  1:674398/1‑130 (MQ=255)
                                                          tCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                             >  4:1126316/1‑148 (MQ=255)
                                                           cACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATc                                                                            <  4:188116/148‑1 (MQ=255)
                                                           cACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATc                                                                            <  1:1022357/148‑1 (MQ=255)
                                                             cTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                                    >  1:168605/1‑122 (MQ=255)
                                                              ttACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  3:989492/1‑98 (MQ=255)
                                                               tACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                         >  3:663717/1‑99 (MQ=255)
                                                               tACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATCTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                         <  4:663717/99‑1 (MQ=255)
                                                                 cTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATg                                                                      <  3:1126316/148‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACt                                                                                                 <  1:130701/119‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                             <  3:198713/139‑1 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGAt                                                                             >  4:198713/1‑139 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAc                                                                    >  3:306352/1‑148 (MQ=255)
                                                                   ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAc                                                                    >  4:108223/1‑148 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAtt                                                                                                                             <  4:794450/89‑1 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGa                                                                                                            >  3:394383/1‑106 (MQ=255)
                                                                     aaTATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGa                                                                                                            <  4:394383/106‑1 (MQ=255)
                                                                      atatCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                              >  2:592148/1‑119 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  1:504248/1‑88 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  1:1111873/1‑88 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           <  1:154634/88‑1 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           >  2:154634/1‑88 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                           <  2:1111873/88‑1 (MQ=255)
                                                                        atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                               <  4:859746/148‑1 (MQ=255)
                                                                            ttGGCCAGGATACTATATGTCATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGc                                                                         <  4:749684/134‑1 (MQ=255)
                                                                            ttGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       >  3:859547/1‑72 (MQ=255)
                                                                            ttGGCCAGGAGACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGAGTTGATCTGc                                                                         >  3:749684/1‑134 (MQ=255)
                                                                             tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCa                                                                                                                                        <  2:315066/70‑1 (MQ=255)
                                                                             tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATa                                                          <  3:108223/148‑1 (MQ=255)
                                                                               gCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaaa                                                                                                               <  4:977969/93‑1 (MQ=255)
                                                                                ccAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAg                                                       >  3:65312/1‑148 (MQ=255)
                                                                                  aGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaa                                                     >  1:312198/1‑148 (MQ=255)
                                                                                   ggATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaa                                                    <  4:306352/148‑1 (MQ=255)
                                                                                    gATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGa                                                                                                                                            <  4:298161/59‑1 (MQ=255)
                                                                                    gATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGa                                                                                                                                            >  3:298161/1‑59 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                            <  3:246840/138‑1 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                            >  4:246840/1‑138 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                            <  4:807536/138‑1 (MQ=255)
                                                                                     aTACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                            >  3:807536/1‑138 (MQ=255)
                                                                                      tACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAata                                                 <  2:625202/148‑1 (MQ=255)
                                                                                          atatGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCa                                                                 >  2:1092089/1‑128 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       <  1:767130/56‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATg                                                                      >  4:73815/1‑121 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                           <  4:891514/132‑1 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                           >  3:891514/1‑132 (MQ=255)
                                                                                            atGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAatat                                                <  3:911873/143‑1 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtg                                                                                                                            <  3:761571/60‑1 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTTAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                           >  1:1072709/1‑125 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCAtt                                    <  3:984912/148‑1 (MQ=255)
                                                                                                   aCGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAAGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGAt                                                           <  2:1072709/125‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                                <  2:701724/69‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCt                                                                                                      <  1:1007407/79‑1 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCt                                                                                                      >  2:1007407/1‑79 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                               >  3:414476/1‑118 (MQ=255)
                                                                                                      aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                               <  4:414476/118‑1 (MQ=255)
                                                                                                        aGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       <  3:627185/44‑1 (MQ=255)
                                                                                                        aGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaa                                                                                                                                       >  4:627185/1‑44 (MQ=255)
                                                                                                           gAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtg                                                                                  >  4:212757/1‑94 (MQ=255)
                                                                                                           gAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtg                                                                                  <  3:212757/94‑1 (MQ=255)
                                                                                                           gAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCttt                                <  2:50866/144‑1 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                                >  1:972184/1‑62 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaa                                                                                                                <  2:972184/62‑1 (MQ=255)
                                                                                                             aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTg                                                                          <  1:434201/100‑1 (MQ=255)
                                                                                                               tCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                                >  4:1017009/1‑108 (MQ=255)
                                                                                                               tCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                                <  3:1017009/108‑1 (MQ=255)
                                                                                                               tCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAg                        >  4:13069/1‑148 (MQ=255)
                                                                                                                       attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATaa                                          <  1:809202/118‑1 (MQ=255)
                                                                                                                        aTGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCa                                                                                               >  4:603190/1‑68 (MQ=255)
                                                                                                                        aTGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCa                                                                                               <  3:603190/68‑1 (MQ=255)
                                                                                                                         tctATCCCGATATAGGGTTAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                                >  1:534357/3‑98 (MQ=255)
                                                                                                                         tctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCaa                                                                <  2:534357/96‑1 (MQ=255)
                                                                                                                         tGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGc                                                                         >  2:970591/1‑89 (MQ=255)
                                                                                                                          ctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTg                                                                          >  2:131895/2‑87 (MQ=255)
                                                                                                                            aTCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCa                                                                        <  4:1081799/87‑1 (MQ=255)
                                                                                                                            aTCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCa                                                                        >  3:1081799/1‑87 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgt                                                                                   >  1:1074730/1‑74 (MQ=255)
                                                                                                                              cccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTg                   >  2:802841/1‑138 (MQ=255)
                                                                                                                               ccTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATa                                           >  2:45700/1‑113 (MQ=255)
                                                                                                                                     tCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATAtct  <  3:13069/148‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
GGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCT  >  NC_000913/2760455‑2760735

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