Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,636 T→C intergenic (+282/+233) b0501 → / ← b0502 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,6360TC97.4% 721.6 / 4.8 232intergenic (+282/+233)b0501/b0502hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (87/139);  total (90/142)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 238 total aligned reads displayed.

TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725
                                                                                          |                                                                                         
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         <  2:1496732/93‑1 (MQ=255)
tCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACgcgt                                                                                         <  1:2908579/93‑1 (MQ=255)
 cAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGtt                                                                                        <  1:455456/93‑1 (MQ=255)
 cAGAGAGTTTACATAGAATGGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGtt                                                                                        <  1:2654877/93‑1 (MQ=255)
   gagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTa                                                                                      >  1:3503930/1‑93 (MQ=255)
   gagagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTa                                                                                      >  2:4333404/1‑93 (MQ=255)
     gagTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAt                                                                                    <  2:3195844/93‑1 (MQ=255)
      agTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAAtg                                                                                   <  2:3999887/93‑1 (MQ=255)
      agTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtg                                                                                   >  2:683613/1‑93 (MQ=255)
        tttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAAtgtg                                                                                 <  2:3796492/93‑1 (MQ=255)
         ttACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGa                                                                                <  2:1096533/93‑1 (MQ=255)
          tACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAg                                                                               <  1:1395620/93‑1 (MQ=255)
           aCATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGt                                                                              <  2:2392916/93‑1 (MQ=255)
            cATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTc                                                                             >  2:1363711/1‑93 (MQ=255)
             aTAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTcc                                                                            <  2:1935428/93‑1 (MQ=255)
              tAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTccc                                                                           <  1:2400510/93‑1 (MQ=255)
              tAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTccc                                                                           <  1:3688259/93‑1 (MQ=255)
              tAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTccc                                                                           <  2:2346182/93‑1 (MQ=255)
                 aaTTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTAt                                                                        <  1:2727817/93‑1 (MQ=255)
                  aTTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATc                                                                       <  1:1211651/93‑1 (MQ=255)
                   ttGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCt                                                                      <  2:4083382/93‑1 (MQ=255)
                     gCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGa                                                                    >  1:4311859/1‑93 (MQ=255)
                      cACTGGTCTTTTACGATATCTGACATTGTGTAATACGTAATCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGaa                                                                   <  2:4219280/93‑1 (MQ=255)
                             cTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaataaat                                                            >  1:2288491/1‑93 (MQ=255)
                             cTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaataaat                                                            <  1:3251486/93‑1 (MQ=255)
                             cTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAaaataaat                                                            <  2:4371150/93‑1 (MQ=255)
                                ttACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAaataat                                                         >  2:2640046/1‑93 (MQ=255)
                                 tACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATc                                                        >  1:564622/1‑93 (MQ=255)
                                  aCGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATcc                                                       <  2:76335/93‑1 (MQ=255)
                                     atatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCTTTc                                                    >  1:4223282/1‑93 (MQ=255)
                                     atatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTc                                                    <  1:4117769/93‑1 (MQ=255)
                                      tatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTcc                                                   >  2:2789331/1‑93 (MQ=255)
                                      tatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTcc                                                   >  1:1699120/1‑93 (MQ=255)
                                      tatCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTcc                                                   >  2:2556845/1‑93 (MQ=255)
                                       atCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCg                                                  >  2:327372/1‑93 (MQ=255)
                                       atCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCg                                                  <  1:206447/93‑1 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 >  2:3648523/1‑93 (MQ=255)
                                        tCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCgg                                                 >  1:3884430/1‑93 (MQ=255)
                                         cTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGa                                                <  1:1040783/93‑1 (MQ=255)
                                         cTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGa                                                >  1:1314763/1‑93 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               >  1:3815250/1‑93 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               <  2:3212519/93‑1 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               >  2:4565792/1‑93 (MQ=255)
                                          tGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAt                                               >  2:1799648/1‑93 (MQ=255)
                                           gACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAtt                                              <  1:3363865/93‑1 (MQ=255)
                                           gACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGAtt                                              <  2:3382015/93‑1 (MQ=255)
                                            aCATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTa                                             >  2:1697688/1‑93 (MQ=255)
                                               ttGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTaaaa                                          <  2:598073/93‑1 (MQ=255)
                                               ttGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTATAATGTGAGACCCTATCTGAAAATAAATAATCCTTCCGGATTaaaa                                          <  2:3140032/93‑1 (MQ=255)
                                                tgtgTAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaat                                         >  1:1900993/1‑93 (MQ=255)
                                                  tgtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataa                                       <  2:2461224/93‑1 (MQ=255)
                                                  tgtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataa                                       >  1:1558792/1‑93 (MQ=255)
                                                   gtAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaa                                      <  1:383492/93‑1 (MQ=255)
                                                    tAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaat                                     <  2:959212/93‑1 (MQ=255)
                                                    tAATACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAaaataaat                                     >  1:4446117/1‑93 (MQ=255)
                                                      aTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAGTCCTTCCGGATTAAAATAAATTc                                   >  2:3229744/1‑93 (MQ=255)
                                                      aTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTc                                   >  2:2455403/1‑93 (MQ=255)
                                                      aTACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTc                                   <  2:3719168/93‑1 (MQ=255)
                                                       tACATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCt                                  >  2:2867268/1‑93 (MQ=255)
                                                        aCATATTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCtt                                 <  1:2961713/93‑1 (MQ=255)
                                                          atatTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGc                               <  1:1788308/93‑1 (MQ=255)
                                                          atatTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGc                               >  1:910203/1‑93 (MQ=255)
                                                           tatTCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGcc                              <  2:2881212/93‑1 (MQ=255)
                                                              tCAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCggg                           >  1:1943282/1‑93 (MQ=255)
                                                               cAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGga                          <  2:3675974/93‑1 (MQ=255)
                                                               cAGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGga                          <  2:2824476/93‑1 (MQ=255)
                                                                aGCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaa                         >  2:4698554/1‑93 (MQ=255)
                                                                 gctaTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  1:3455137/90‑1 (MQ=255)
                                                                 gCCATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaa                        <  1:2200004/93‑1 (MQ=255)
                                                                  ccATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaag                       <  1:2005304/93‑1 (MQ=255)
                                                                   cATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaaga                      >  1:687573/1‑93 (MQ=255)
                                                                   cATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaaga                      >  1:15040/1‑93 (MQ=255)
                                                                   cATGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaaga                      >  1:1410901/1‑93 (MQ=255)
                                                                    aTGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaa                     >  1:4152655/1‑93 (MQ=255)
                                                                     tGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaaa                    <  2:159599/93‑1 (MQ=255)
                                                                     tGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaaa                    <  1:3519498/93‑1 (MQ=255)
                                                                     tGCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGgaaagaaa                    <  2:3004938/93‑1 (MQ=255)
                                                                      gCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGaaagaaag                   <  1:2945396/93‑1 (MQ=255)
                                                                      gCATTAATTAAAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGaaagaaag                   <  2:280611/93‑1 (MQ=255)
                                                                           aattaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaa              <  2:2184200/93‑1 (MQ=255)
                                                                            attaaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGaaa             <  2:3070887/93‑1 (MQ=255)
                                                                              taaAATGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATa           <  2:2603104/93‑1 (MQ=255)
                                                                               aaaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaa          >  2:3542371/1‑93 (MQ=255)
                                                                                aaaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATaaa         <  1:2819338/93‑1 (MQ=255)
                                                                                 aaTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAc        <  1:1800050/93‑1 (MQ=255)
                                                                                  aTGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAAcc       <  2:1809941/93‑1 (MQ=255)
                                                                                   tGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCa      >  2:890486/1‑93 (MQ=255)
                                                                                   tGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCa      <  2:2133895/93‑1 (MQ=255)
                                                                                   tGTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCa      <  1:2532546/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     <  1:1075344/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     <  2:492094/93‑1 (MQ=255)
                                                                                    gTTACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAt     >  2:2225874/1‑93 (MQ=255)
                                                                                     ttACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCAtt    >  1:3992613/1‑93 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   <  2:1187150/93‑1 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   <  1:2477922/93‑1 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   >  2:3952641/1‑93 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   <  1:2556845/93‑1 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   <  2:3503930/93‑1 (MQ=255)
                                                                                      tACGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTa   <  2:3069983/93‑1 (MQ=255)
                                                                                       aCGCGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAg  >  1:4723019/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCAGAGAGTTTACATAGAATTGCACTGGTCTTTTACGATATCTGACATTGTGTAATACATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAG  >  NC_000913/528546‑528725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: