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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,760,5890CA 7.0% 579.6 / 20.4 172R15R (CGG→AGG) yfjJCP4‑57 prophage; protein YfjJ
Reads supporting (aligned to +/- strand):  ref base C (72/88);  new base A (9/3);  total (81/91)
Fisher's exact test for biased strand distribution p-value = 6.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 175 total aligned reads displayed.

TCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGC  >  NC_000913/2760441‑2760738
                                                                                                                                                    |                                                                                                                                                     
tctAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                                                                                                                                      >  2:82822/1‑150 (MQ=255)
 ctAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCgg                                                                                                                                                     >  1:341629/1‑150 (MQ=255)
      cccAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACAGATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAg                                                                                                                                                <  1:581671/150‑1 (MQ=255)
        cAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGaa                                                                                                                                              >  2:1127336/1‑150 (MQ=255)
             cGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaaa                                                                                                                                         >  2:1071166/1‑150 (MQ=255)
               gATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTACAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATACATATATCGGGATAGAATCAAACa                                                                                                                                       >  1:973691/1‑150 (MQ=255)
                 tAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAgg                                                                                                                                     <  1:1396992/150‑1 (MQ=255)
                     aGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAt                                                                                                                                 <  2:952210/150‑1 (MQ=255)
                      gCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAtt                                                                                                                                <  1:201657/150‑1 (MQ=255)
                        ttGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtga                                                                                                                              >  2:498477/1‑150 (MQ=255)
                          gTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                            <  1:591019/150‑1 (MQ=255)
                          gTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatg                                                                                                                            <  1:2726/150‑1 (MQ=255)
                           tCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatga                                                                                                                           >  1:16758/1‑150 (MQ=255)
                           tCATCTTGATAGCTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatga                                                                                                                           >  2:910386/1‑150 (MQ=255)
                            cATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGAc                                                                                                                          >  2:1218875/1‑150 (MQ=255)
                             aTCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACt                                                                                                                         >  2:1428806/1‑150 (MQ=255)
                             aTCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACt                                                                                                                         >  1:1502401/1‑150 (MQ=255)
                              tCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTc                                                                                                                        >  1:718797/1‑150 (MQ=255)
                              tCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTc                                                                                                                        >  1:480927/1‑150 (MQ=255)
                                ttGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTccc                                                                                                                      <  1:82822/150‑1 (MQ=255)
                                    tAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTaaa                                                                                                                  >  1:1052991/1‑150 (MQ=255)
                                     aGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAAc                                                                                                                 >  1:1515850/1‑150 (MQ=255)
                                      ggTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACg                                                                                                                <  2:1048083/150‑1 (MQ=255)
                                      ggTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACg                                                                                                                <  2:1039228/150‑1 (MQ=255)
                                       gTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGa                                                                                                               <  1:1402700/150‑1 (MQ=255)
                                         aaTATATAGGTACATAGACACATCAGAACGATCCCTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAAc                                                                                                             <  1:1173317/150‑1 (MQ=255)
                                           tatataGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACAt                                                                                                           <  2:629772/150‑1 (MQ=255)
                                               taGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCATGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                                       <  1:1395128/150‑1 (MQ=255)
                                               taGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                                       <  2:804707/150‑1 (MQ=255)
                                               taGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTa                                                                                                       >  2:1153314/1‑150 (MQ=255)
                                                  gTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACt                                                                                                    >  2:450496/1‑150 (MQ=255)
                                                   tACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTc                                                                                                   >  2:553805/1‑150 (MQ=255)
                                                     cATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCaa                                                                                                 >  2:220158/1‑150 (MQ=255)
                                                      aTAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAAGATCTTAACTCAAt                                                                                                <  2:57923/150‑1 (MQ=255)
                                                      aTAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAAt                                                                                                >  1:488167/1‑150 (MQ=255)
                                                       tAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATg                                                                                               >  2:457776/1‑150 (MQ=255)
                                                             caTCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTa                                                                                         >  2:1148172/1‑150 (MQ=255)
                                                             caTCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTa                                                                                         <  1:466054/150‑1 (MQ=255)
                                                               tCAGAACGATCACTGACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACg                                                                                       <  1:365980/150‑1 (MQ=255)
                                                                 aGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtg                                                                                     <  2:840226/150‑1 (MQ=255)
                                                                 aGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtg                                                                                     >  2:390145/1‑150 (MQ=255)
                                                                  gAACGATCACTTACTCCAATATCATTTGCCATGATACTATATTACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtgt                                                                                    <  2:1454606/150‑1 (MQ=255)
                                                                  gAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACgtgt                                                                                    <  2:1393654/150‑1 (MQ=255)
                                                                     cGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGa                                                                                 >  1:1285053/1‑150 (MQ=255)
                                                                         cACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTg                                                                             <  2:1114894/150‑1 (MQ=255)
                                                                          aCTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGc                                                                            >  1:667430/1‑150 (MQ=255)
                                                                           cTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCa                                                                           <  2:1349688/150‑1 (MQ=255)
                                                                            ttACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCAt                                                                          <  2:1323879/150‑1 (MQ=255)
                                                                               cTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAc                                                                       >  1:284402/1‑150 (MQ=255)
                                                                                 ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAccc                                                                     <  2:509753/150‑1 (MQ=255)
                                                                                 ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAccc                                                                     <  2:602631/150‑1 (MQ=255)
                                                                                 ccAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGAccc                                                                     <  2:617785/150‑1 (MQ=255)
                                                                                    atatCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAt                                                                  <  2:710817/150‑1 (MQ=255)
                                                                                     tatCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAAtt                                                                 <  1:910386/150‑1 (MQ=255)
                                                                                      atCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTg                                                                >  2:381551/1‑150 (MQ=255)
                                                                                       tCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                               >  2:794387/1‑150 (MQ=255)
                                                                                       tCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTAATTATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGa                                                               <  2:588149/150‑1 (MQ=255)
                                                                                         aTTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATa                                                             >  2:1273421/1‑150 (MQ=255)
                                                                                         aTTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATa                                                             >  1:1160299/1‑150 (MQ=255)
                                                                                         aTTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATa                                                             <  2:241110/150‑1 (MQ=255)
                                                                                           tGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACa                                                           <  1:739965/150‑1 (MQ=255)
                                                                                            ggCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAg                                                          <  1:1435761/150‑1 (MQ=255)
                                                                                                    tACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                  >  1:246056/1‑150 (MQ=255)
                                                                                                    tACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTAAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACTGAAAATATg                                                  <  2:786598/150‑1 (MQ=255)
                                                                                                    tACTATATGACATACGATAGTGAATTCGGAGCACATGTATCCCTATATCGGGATAGAATCAAACAGGGTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                  <  2:779450/150‑1 (MQ=255)
                                                                                                      cTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGa                                                <  2:329182/150‑1 (MQ=255)
                                                                                                        atatGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATa                                              >  2:478548/1‑150 (MQ=255)
                                                                                                         tatGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATaa                                             >  1:1088319/1‑150 (MQ=255)
                                                                                                          atGACATACGATAGTGAATTCGGATCACATGTATCACTATATCGTGATAGAATCAAACAGATTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAAc                                            >  2:366650/1‑150 (MQ=255)
                                                                                                              cATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                         <  2:259294/133‑1 (MQ=255)
                                                                                                              cATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                         >  1:259294/1‑133 (MQ=255)
                                                                                                               aTACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCAtt                                       >  1:642649/1‑150 (MQ=255)
                                                                                                                    aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTc                                  <  1:1465340/150‑1 (MQ=255)
                                                                                                                    aTAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTc                                  <  2:642649/150‑1 (MQ=255)
                                                                                                                       gTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAc                               <  2:466139/150‑1 (MQ=255)
                                                                                                                           aTTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAg                           >  1:1056247/1‑150 (MQ=255)
                                                                                                                              cGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTg                                                                >  1:43886/1‑110 (MQ=255)
                                                                                                                              cGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGt                        <  1:1328204/150‑1 (MQ=255)
                                                                                                                               ggATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                       <  2:146909/150‑1 (MQ=255)
                                                                                                                               ggATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                       <  2:1175117/150‑1 (MQ=255)
                                                                                                                                   tgattctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAct                   >  2:497975/7‑150 (MQ=255)
                                                                                                                                    gattctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                  >  2:431943/6‑150 (MQ=255)
                                                                                                                                    gattctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                  >  1:608230/6‑150 (MQ=255)
                                                                                                                                    acaTGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                  <  2:1415673/150‑1 (MQ=255)
                                                                                                                                    acaTGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACAACTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                  >  2:175475/1‑150 (MQ=255)
                                                                                                                                     attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAggg                         <  1:1279870/138‑1 (MQ=255)
                                                                                                                                     attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActct                 >  1:1300788/5‑150 (MQ=255)
                                                                                                                                     attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActct                 >  2:1248448/5‑150 (MQ=255)
                                                                                                                                     attctATCCCGATATAGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCAAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActct                 >  2:682633/5‑150 (MQ=255)
                                                                                                                                       tGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTgg               <  2:1056247/150‑1 (MQ=255)
                                                                                                                                        gTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGt              <  1:1218875/150‑1 (MQ=255)
                                                                                                                                         tATCCCTTTATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTg             <  2:718797/150‑1 (MQ=255)
                                                                                                                                         tATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTg             <  2:758/150‑1 (MQ=255)
                                                                                                                                         tATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTg             >  1:1353038/1‑150 (MQ=255)
                                                                                                                                         tATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTg             >  1:765983/1‑150 (MQ=255)
                                                                                                                                           tCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATAAAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCt           >  2:1465314/1‑150 (MQ=255)
                                                                                                                                           tCCCTATATCGGCATAGAATCAAACAGGTAATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTCCGAGTTGATCTGCATGACCCCAGTGATACAGAAAATATGGATAGGCCATTCTTTCCACCCAGGGGTGAGTCTGGAGcg           >  2:527572/1‑149 (MQ=255)
                                                                                                                                                    cGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGc  >  2:74410/1‑150 (MQ=255)
                                                                                                                                                    cGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGc  <  1:5311/150‑1 (MQ=255)
                                                                                                                                                    cGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGc  >  2:421072/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGC  >  NC_000913/2760441‑2760738

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